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HSRNOISE User’s and Theoretical Manual 


1.0 Introduction 

From 1989 to 1999, NASA and US industry partners conducted the High Speed Research 
(HSR) program to advance technology toward the development of a supersonic 
commercial transport. One element of the program was to ensure that the aircraft could be 
certified according to FAA FAR 36 Noise Certification regulations. During the HSR 
program, system noise studies were made using the NASA Fangley Aircraft Noise 
Prediction Program (ANOPP) and the Boeing Multiple Component Predictor (MCP) 
program. Since the source for both of these codes is proprietary to their respective 
agencies, neither was suitable to be archived at the end of the HSR program. The 
HSRNOISE computer code, described in this report, was developed to document the noise 
prediction effort of the HSR program and to provide a code for continued research on a 
supersonic commercial transport. 

Research was conducted on several engine concepts during Phase I of the HSR program. 
Among them were the turbine bypass engine, the variable cycle engine, and the mixed 
flow turbofan. The mixed flow turbofan was ultimately chosen as the primary engine for 
the Phase II studies. For this engine, jet noise was acknowledged as a noise source that 
would require substantial suppression. Suppressors were designed to promote rapid 
mixing of the supersonic jet with aspirated ambient flow by a multilobe mixer inside an 
ejector. The new propulsion concepts warranted the development of new prediction 
methods to study the noise impact of the High Speed Civil Transport (HSCT). Fan and jet 
noise algorithms were developed specifically for the mixed flow turbofan and mixer 
ejector nozzle. Combustor, turbine, and airframe noise prediction algorithms were 
obtained from published methods and correlated with model and/or full-scale data. 

This report provides instruction on executing the HSRNOISE code, documents the noise 
prediction algorithms, and provides comparisons of predicted noise levels from the MCP 
and HSRNOISE codes for the Technology Concept (TC) aircraft. A flow chart of the 
HSRNOISE code is provided in figure 1.0. Whenever possible, algorithms provided by 
Boeing from the MCP system were incorporated into the HSRNOISE code. Algorithms 
obtained from the MCP system include subroutines in the Aircraft Flight Definition 
section, the jet noise subroutines (JN8C4), and the atmospheric attenuation subroutines. 
Modern Technologies Corporation developed an alternative mixer ejector noise prediction 
method for General Electric Aircraft Engines, which has also been incorporated into the 
HSRNOISE code. The remainder of the HSRNOISE code was developed at NASA 
Langley. 
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2.0 User's Manual 


This section provides instructions on creating input files, and executing the HSRNOISE 
code. The HSRNOISE code was developed at NASA Langley on a DEC Alpha 
workstation using the UNIX OAS1 version 4.0 operating system. The code was compiled 
with the Digital FORTRAN F77/90 compiler. 

2.1 Creating the Executable 

All subroutines required to create an executable are provided in a tar file. Once the 
FORTRAN files have been extracted, the executable hsrnoise is created by entering 
the following command: 

> make -f makefile . hsrnoise 

HSRNOISE prompts the user for the names of the input and output files. To execute the 
program with an input file named sample . input enter: 

> hsrnoise (return) 

ENTER INPUT FILENAME: sample. input (return) 

ENTER OUTPUT FILENAME: sample . output (return) 

An output file named sample . output will be created which contains flight path 
information, ambient weather conditions as a function of altitude, and noise levels for each 
noise source. The name of the output file is arbitrary. 
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2.2 Creating an Input File 

Input to the HSRNOISE code is provided through NAMELIST statements. The following 
1 1 NAMELIST groups have been established and must be input in the order specified: 


GEOM 

Defines the geometric relationship between the aircraft and 
the observer 

WEATHER 

Defines ambient weather conditions at ground level 

FLIPATH 

Defines the flight path parameters 

SOURCE 

Defines noise source and propagation adjustment for 
atmospheric attenuation, ground reflections, and EGA 

DELTAS 

Defines spectral adjustment 
(Required if IDELTA = TRUE) 

FANIN 

Defines fan noise parameters 
(Required if IF AN = TRUE) 

COREIN 

Defines combustor noise parameters 
(Required if ICORE = TRUE) 

TURIN 

Defines turbine noise parameters 
(Required if ITUR = TRUE) 

JET 1 IN 

Defines jet noise parameters 

(Required if IJET = TRUE, JETMETHOD = 1) 

JET2IN 

Defines jet noise parameters 

(Required if IJET = TRUE, JETMETHOD = 2) 

AIRFIN 

Defines airframe noise parameters 
(Required if IAFRAM = TRUE) 


The parameters defined within each namelist are described in the following sections. Care 
should be taken to ensure that the parameters are in the proper units. 
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2.2.1 NAMELIST GEOM 


Variable 

Units 

Default 

Definition 

ANGLE 

degrees 

10,20,. ..,170 

Emission angle measured from engine inlet axis 

HMIC 

ft 

4. 

Height of microphone above ground level 

NANG 


17 

Number of elements in array ANGLE 

SL 

ft 

0 . 

Sideline distance 


2.2.2 NAMELIST WEATHER 


Variable 

Units 

Default Definition 

PSIAMB 

lb/in 2 

14.7 

Ambient pressure at ground level 

RH 

% 

70 

Relative humidity at ground level 

SIGMA 

lbm/sec-ft 3 

9.4E3 

Ground flow resistance 

TFAMB 

°F 

77 

Ambient temperature at ground level 

2.2.3 NAMELIST FLIPATH 


Variable 

Units 

Default 

Definition 

ALTEVO 

ft 

700. 

Aircraft altitude at overhead 

FPA 

degrees 

0 . 

Flight Path Angle. Angle between airplane 
instantaneous flight path and the horizontal 
direction. Positive for Climb. 

PAE 

degrees 

0 . 

Angle between engine axis and horizontal 
direction. Equals the sum of the airplane 
flight path angle (FPA), the airplane body 
angle-of-attack, and the angle between the 
airplane reference axis and the engine axis. 
Positive for inlet up. 

VAIR 

ft/sec 

0 . 

Aircraft velocity 
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2.2.4 NAMELIST SOURCE 


Variable 

Default 

Definition 

IAFRAM 

TRUE 

Airframe source flag, (TRUE / FALSE) 

IATM 

1 

Atmospheric absorption flag 
= 0 ; no atmospheric absorption 
= 1 ; SAE ARP866A 
= 2 ; Shields & Bass 

ICORE 

TRUE 

Core source flag, (TRUE / FALSE) 

IDELTA 

FALSE 

User supplied suppression flag, (TRUE / FALSE) 

IEGA 

TRUE 

Extra ground attenuation flag, (TRUE / FALSE) 

IFAN 

TRUE 

Fan source flag, (TRUE / FALSE) 

IGRND 

TRUE 

Ground reflection flag, (TRUE / FALSE) 

IJET 

TRUE 

Jet source flag, (TRUE / FALSE) 

IJETSHLD 

TRUE 

Jet shielding flag, (TRUE / FALSE) 

ITUR 

TRUE 

Turbine source flag, (TRUE / FALSE) 

JETMETHOD 

1 

Jet method flag 
= 1, Boeing JN8C4 Method 
= 2, Stone’s Method 

NENG 

1 

Number of engines 





2.2.5 NAMELIST DELTAS 


Variable 

Units 

Default 

Definition 

FANDIBB 

dB 

17*0.0 

Fan broadband directivity delta array 

FANSPBB 

dB 

24*0.0 

Fan broadband spectrum delta array 

FANBPF 

dB 

0.0 

Fan first harmonic delta 

FAN2BPF 

dB 

0.0 

Fan second harmonic delta 

FAN3BPF 

dB 

0.0 

Fan third harmonic delta 

FAN4BPF 

dB 

0.0 

Fan fourth harmonic delta 

CORDIBB 

dB 

17*0.0 

Combustor broadband directivity delta 

CORSPBB 

dB 

24*0.0 

Combustor broadband spectrum delta 

TURDIBB 

dB 

17*0.0 

Turbine broadband directivity delta 

TURSPBB 

dB 

24*0.0 

Turbine broadband spectrum delta 

TURDITN 

dB 

17*0.0 

Turbine tone directivity delta 

TURSPTN 

dB 

24*0.0 

Turbine tone spectrum delta 

AIRFDICW 

dB 

17*0.0 

Airframe directivity delta (Clean Wing) 

AIRFDILE 

dB 

17*0.0 

Airframe directivity delta (Leading Edge Slat) 

AIRFDITEI 

dB 

17*0.0 

Airframe directivity delta (Trailing Edge Inboard Slat) 

AIRFDITEO 

dB 

17*0.0 

Airframe directivity delta (Trailing Edge Outboard Slat) 

AIRFDIHT 

dB 

17*0.0 

Airframe directivity delta (Horizontal Tail) 

AIRFDIVT 

dB 

17*0.0 

Airframe directivity delta (Vertical Tail) 

AIRFDIMG 

dB 

17*0.0 

Airframe directivity delta (Main Gear) 

AIRFDING 

dB 

17*0.0 

Airframe directivity delta (Nose Gear) 

AIRFSPCW 

dB 

24*0.0 

Airframe spectrum delta (Clean Wing) 

AIRFSPLE 

dB 

24*0.0 

Airframe spectrum delta (Leading Edge Slat) 

AIRFSPTEI 

dB 

24*0.0 

Airframe spectrum delta (Trailing Edge Inboard Slat) 

AIRFSPTEO 

dB 

24*0.0 

Airframe spectrum delta (Trailing Edge Outboard Slat) 

AIRFSPHT 

dB 

24*0.0 

Airframe spectrum delta array (Horizontal Tail) 

AIRFSPVT 

dB 

24*0.0 

Airframe spectrum delta array (Vertical Tail) 

AIRFSPMG 

dB 

24*0.0 

Airframe spectrum delta array (Main Gear) 

AIRFSPNG 

dB 

24*0.0 

Airframe spectrum delta array (Nose Gear) 
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2.2.6 NAMELIST FANIN 


Variable 

Units 

Default 

Definition 

AFAN 

ft 2 

20 

Fan inlet annulus area 

FANRPM 

rev/min 

5000 

Fan RPM 

IDOOR 


0 

Flag for fan doors open or closed 
= 0 ; doors closed 
= 1 ; doors open 

ISHAPE 


0 

Flag for 2D Bifurcated Inlet Shape Correction 
= 0 ; No correction 
= 1 ; Correction 

NB1 


25 

Number of fan blades 

UT 

ft/sec 

1000 

First stage corrected tip speed 
(Equation 4 page 21) 

2.2.7 NAMELIST COREIN 


Variable 

Units 

Default 

Definition 

AB 

in 2 

100.0 

Burner cross sectional area 

BL 

ft 

1.0 

Burner length 

ICR 


1 

Burner type 

= 1; Can (JT3D, JT4, JT8D) 
= 2; Annular (JT9D, JT10D) 

NFN 


1 

Number of fuel nozzles 

PLD 


1.0 

Transmission loss factor 

P3 

lb/ in 2 

111.8 

Burner inlet pressure 

P4 

lb/ in 2 

104.2 

Burner exit pressure 

P7 

lb/ in 2 

19.4 

Low pressure turbine exit total pressure 

T3 

°R 

1007.5 

Burner inlet total temperature 

T4 

°R 

1883.5 

Burner exit total temperature 

T7 

°R 

1172.3 

Low pressure turbine exit total temperature 

W3 

lbm/sec 

450.0 

Burner air mass flow rate 

WFP4 

(lbm/sec)/ 
(lb/ ft 2 ) 

0.0155 

Fuel-flow-rate to burner-exit pressure ratio at takeoff 
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2.2.8 NAMELIST TURIN 


Variable 

Units 

Default 

Definition 

CSS 

% 

100 

Stator axial chord to stator rotor spacing ratio 

MDOT 

lbm/sec 

200 

Primary mass flow rate at turbine exit 

NBLADE 


60 

Number of turbine blades is last stage 

TS 

°R 

1300 

Turbine exit static temperature 

TURCAL 

dB 

10 

Turbine calibration constant 

TURRPM 

rev/min 

3600 

Last stage turbine rotational speed 

VTR 

ft/sec 

850 

Relative tip speed of the last rotor of the 
turbine. If unknown use 0.7 times the tip speed. 
(Equation 4 page 36) 


n 






2.2.9 NAMELIST JET1IN (Boeing JN8C4 Methodology) 


Variable 

Units 

Default 

Definition 

APRI 

ft 2 

8.6562 

Primary nozzle exit area 

ASEC 

ft 2 

16.832 

Secondary nozzle exit area 

ATHP 

ft 2 

8.5466 

Throat area of the primary nozzle 

CFG 


0.95 

Mixer-ejector thrust coefficient 

DELMIX 

dB 

0 

Mixed jet suppression 

DELPE 

dB 

0 

Primary external jet suppression 

DELPI 

dB 

0 

Primary internal jet suppression 

DIVANG 

degrees 

0 

Spoke divergence angle. Outer tangent line relative to 
ejector flow surface at nozzle downstream of mixer 
exit. 

EJASP 


1.5 

Mixer-ejector aspect ratio major/minor axis 

EJD 

ft 

5.274 

Equivalent inner diameter of the ejector at the mixer 
exit 

EJL 

ft 

13.3 

Ejector length 

EJLIN 

ft 

9.0 

Length of ejector lining 

FLIN 

Hz 

2000 

Liner design center frequency 

PEN 


0.925 

Mixer spoke penetration to nozzle full height ratio 

PLUGD 

ft 

0 

Plug diameter at nozzle exit station 

PSIO 

degrees 

0 

Rectangular exit plane major axis roll angle relative to 
wing plane 

SPOKE 


18.0 

Number of spokes or lobes of the mixer 

TPRI 

°R 

1553.0 

Primary jet total temperature 

TSEC 

°R 

529.0 

Secondary jet total temperature 

VPRI 

ft/sec 

2378.0 

Primary jet fully expanded velocity 

VSEC 

ft/sec 

449.0 

Secondary jet velocity 

WPRI 

lbm/sec 

811.3 

Primary jet mass flow rate 

WSEC 

lbm/sec 

566.9 

Secondary jet mass flow rate 

XMAR 


0.97 

Ejector exit area to mixing plane area ratio 
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2.2.10 NAMELIST JET2IN (Stone Methodology) 


Variable 

Units 

Default 

Definition 

APT 

ft 2 

8.5466 

Throat area of the primary nozzle 

AS 

ft 2 

16.832 

Secondary nozzle flow area 

CER 



Core expansion ratio 

DHP 

ft 

0.5 

Hydraulic diameter of primary flow 
DHP = 4 * Area / Perimeter 

DL 

dB/ft 

0.7 

Peak perforate suppression at FPK 

DM 


0 

Degree of mixing (calculated if not supplied) 

FPK 

Hz 

1995 

Peak frequency of broad-band suppression 

GAMMAC 


0 . 

Specific heat ratio of primary stream. Internally 
calculated if set to 0. 

HEX 

ft 

4.0 

Ejector height at nozzle exit 

IEX 


0 

Mixed jet parameter flag for TEX & VEX 
=0 calculate internally 
= 1 use input values 

ISUPPR 


0 

Suppression flag 

=0 negate suppression of internal noise sources 
= 1 suppress internal noise sources 

LBE 

ft 

10. 

Length to end of bulk treatment from primary exit 

LBS 

ft 

0 . 

Length to start of bulk treatment from primary exit 

LE 

ft 

10. 

Ejector length 

LPE 

ft 

0 . 

Length to end of perforate treatment from primary exit 
For perforate over bulk treatment, add bulk length to 
perforate 

LPS 

ft 

0 . 

Length to start of perforate treatment from primary exit 

MMC 


0 . 

Molecular mass (weight) of primary stream. 
Calculated internally if set to zero 

MPD 


1.2 

Design Mach number of primary nozzle 

PC 

lb/ft 2 


Total pressure of primary stream upstream of nozzle 

PEN 


0.925 

Mixer spoke penetration to nozzle full height ratio 

SAR 



Ratio of total mixing area to primary nozzle throat 

SUPPK 

dB/ft 

1.6628 

Peak bulk suppression 

TC 

°R 


Total temperature of primary stream upstream of 
nozzle 

TEX 

°R 


Mixed exit total temperature 

VEX 

ft/sec 


Mixed exit velocity 

WEX 

ft 

6.0 

Ejector width at nozzle exit 

WSWP 



Pumping ratio 
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2.2.11 NAMELIST AIRFIN 


Variable 

Units 

Defaults 

Definition 

ICWING 


2 

Flag for clean wing component 
= 0 ; Do not compute 
= 1 ; Clean Wing 
= 2 ; Super Clean Wing 

IVTAIL 


1 

Flag for vertical tail component 
= 0 ; Do not compute 
= 1 ; Compute 

IHTAIL 


1 

Flag for horizontal tail component 
= 0 ; Do not compute 
= 1 ; Compute 

ILEDGE 


1 

Flag for leading edge slat component 
= 0 ; Do not compute 
= 1 ; Compute 

ITEIN 


1 

Flag for inboard trailing edge flap component 
= 0 ; Do not compute 
= 1 ; Compute 

ITEOUT 


1 

Flag for outboard trailing edge flap component 
= 0 ; Do not compute 
= 1 ; Compute 

INGEAR 


1 

Flag for nose gear component 
= 0 ; Do not compute 
= 1 ; Compute 

IMGEAR 


1 

Flag for main gear component 
= 0 ; Do not compute 
= 1 ; Compute 

AW 

ft 2 

8878.0 

Wing area 

BW 

ft 

154.3 

Wing span 

AV 

ft 2 

435.0 

Vertical tail area 

BV 

ft 

16.8 

Vertical tail span 

AH 

ft 2 

726.0 

Horizontal tail area 

BH 

ft 

71.0 

Horizontal tail span 

AFIN 

ft 2 

100.0 

Wing inboard trailing edge flap area 

BFIN 

ft 

10.0 

Inboard flap span 

FLAIN 

degrees 

10.0 

Inboard flap deflection angle 
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NFSIN 


1 

Number of inboard trailing edge slots 

AFOUT 

ft 2 

200.0 

Wing outboard trailing edge flap area 

BFOUT 

ft 

20.0 

Outboard flap span 

FLAOUT 

degrees 

10.0 

Outboard flap deflection angle 

NFS OUT 


1 

Number of outboard trailing edge slots 

NNG 


1 

Number of nose landing gear 

NWNG 


2 

Number of wheels per nose landing gear 

DNG 

ft 

3.3 

Nose gear wheel diameter 

RLNG 

ft 

13.7 

Length of nose landing gear strut 

NMG 


2 

Number of main landing gear 

NWMG 


8 

Number of wheels per main landing gear 

DMG 

ft 

3.7 

Main gear wheel diameter 

RLMG 

ft 

16.0 

Length of main landing gear strut 
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3.0 Theoretical Manual 


The HSRNOISE program predicts one-third-octave band sound pressure levels for a 
supersonic commercial transport powered by mixed flow turbofan engines with 
rectangular mixer-ejectors. Total aircraft noise is modeled as a combination of fan, core, 
turbine, jet, and airframe noise sources. This section provides a description of the noise 
prediction methods used in the HSRNOISE program. A brief description of the methods 
used to determine the geometry, weather, and flight path is also provided. 

3.1 Geometry (GEOM) 

The relationship between the aircraft and the observer is established by the parameters 
defined in the GEOM namelist. Points on the flight path where noise predictions are to be 
made and subsequently propagated to the observer are defined by an array of emission 
angles. The variable ANGLE contains the array of desired emission angles. Each angle is 
checked to determine if it corresponds to a point that is above ground. Any angle not 
passing the test is removed from the array and a revised emission angle array is defined 
for subsequent noise calculation. The observer location and height are defined by the 
parameters SL and HMIC. 

3.2 Atmospheric Properties (WEATHER) 

Ambient weather conditions are required for the calculation of atmospheric absorption. 
The ambient pressure, ambient temperature, relative humidity, and ground flow resistance 
are input at ground level. Changes in atmospheric properties as a function of altitude are 
computed from the conditions at ground level and a temperature lapse rate of 3.566 
degrees per 1000 feet. The relative humidity does not change with altitude. 

Ambient sound speed: c a (z)= 49.01-^(r fl + 459.7)-.003566z (1) 

Ambient pressure: P a (z) = P a \z = Oj 

144P (A 

Ambient density: P„(z) = — 7 c (3) 

j r a \ > R ^ T + 4597 _ qo 3566 z) 

The ground flow resistivity is an important parameter in the calculation of the ground 
reflection attenuation. Table 3.2.1 provides a range of values for ground flow resistivity 
typical of aeroacoustic predictions. 


[T a +459.7) -,003566z 


T + 459.7 


( 2 ) 
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Table 3.2.1 Range of Flow Resistivity for Various Ground Surfaces 


Surface Type 

Pa- sec/nr 

lbm/sec-ft 3 

Grass 

Asphalt 

Concrete 

1.5 - 3 x 10 3 
1.5 - 3x 10 7 
3 - 10 x 10 7 

9.4- 19.0 x 10 J 
9.4- 19.0 x 10 6 
1.9- 6.2 xlO 8 


from NASA RP 1258 vol. 2 Table 1 page 66 


3.3 Flight Path (FLIPATH) 

When predicting noise for a single observer, a complete description of the flight path is 
not necessary. Only the portion of the flight path that is significant to that observer is 
required. The flight path is defined as a straight line with the flight path angle, pitch 
angle, and aircraft velocity remaining constant. The altitude is established by defining the 
altitude when the aircraft is directly above the observer. 
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3.4 Aircraft Noise Sources (SOURCES) 


The parameters in the SOURCE namelist allow the user to select the noise sources that 
will model the total aircraft noise. Noise is computed for the fan, core, turbine, jet, and 
airframe when the appropriate source noise flag is set to TRUE. The input also allows the 
user to control certain propagation effects, such as atmospheric attenuation, ground 
reflection attenuation, and extra ground attenuation. These propagation effects are 
computed when the proper parameters are set to TRUE. Atmospheric attenuation is 
controlled by the IATM parameter, which has the following options: 

• no atmospheric attenuation 

• atmospheric attenuation determined using the SAE ARP866A method 

• atmospheric attenuation determined using the Shields & Bass method 

Two criteria are used to determine the attenuation frequency for the one-third-octave band 
analysis. For the SAE ARP866A method, the attenuation coefficient is determined using 
the following criteria: 

• one-third-octave band center frequency for frequencies below 4000 Hz. 

• one-third-octave band lower band edge frequency for frequencies 5000 Hz. and 
above 

When the Shields & Bass method is selected, the atmospheric attenuation coefficients are 
determined using the one-third-octave band center frequency. 
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3.5 Suppression (DELTAS) 


The DELTAS namelist allows the spectra of a noise source to be altered by a suppression 
function. Suppression is applied as 

^^^suppressed Ainsuppressed + A (f,0) (1) 

where A is a function of frequency and emission angle . The suppression function is 
calculated as a product as indicated by equation 2. 

A = A(0)B(f) (2) 


The function A is a multiplier and is usually positive. The function B provides the 
spectral suppression and is usually negative. Note that the sign of the suppression 
function can be positive or negative depending on whether noise is to be added or 
subtracted from the original spectra. 

There are several important points to note concerning the suppression function. 

• The default values for the fan harmonics (see parameters FANBPF, FAN2BPF, 
FAN3BPF, and FAN4BPF in section 2.2.5) are not arrays but single valued 
parameters. Instead of the suppression function, each fan harmonic is suppressed by a 
constant value that is not a function of the emission angle. 

• Jet noise suppression has been incorporated into the jet noise methods. See JET1IN 
and JET2IN for details. 


• Airframe noise is comprised of eight subcomponents. A unique suppression function 
can be applied to each subcomponent. When predicting airframe noise for a high 
speed civil transport other than the TC aircraft, an adjustment to the subsonic airframe 
prediction method will be required. The suppression function derived for the TC 
aircraft is defined in section 3.11. 
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3.6 Fan Noise 


The fan noise prediction method is based on the empirical formulation of XF-120 fan data 
derived by Gliebe [1]. The method assumes that inlet fan noise radiates from the first 
stage only and that noise from the 2nd and 3rd stages do not contribute. The method 
employs empirical functions to define the broadband and discrete tone noise components. 
The broadband and tone noise peak sound pressure levels are correlated as a function of 
fan corrected tip speed and polar directivity angle. The broadband spectral function is 
predicted using the spectral shape function from the Heidmann fan noise prediction 
method [2]. The total noise is obtained by summing the broadband and tone noise mean- 
square acoustic pressures for each frequency and polar directivity angle. 

Symbols 

a broadband coefficient for 2-D inlet correction 

A fan first stage rotor inlet annulus area, (ft ) ( input parameter - AFAN) 

Ajj broadband correlation coefficients 

c tone coefficient for 2-D inlet correction 

C I; tone correlation coefficients 

d fan rotor diameter, ft 

/ frequency, Hz 

f b blade passing frequency, Hz 

n harmonic number 

N h number of first stage fan blades ( input parameter - NB1) 

N 2 fan rotational speed, (rpm) ( input parameter - FANRPM) 

p ref reference pressure, (AMI xKY 1 lb I ft 1 ) 

r s distance from source to observer, (ft) 

SI spectrum index function, dB 

T a ambient temperature, °R 

U nondimensional corrected tip speed, re U ref 

U ref reference tip speed, (1000 ft/ sec) 

U t first stage corrected tip speed, (ft/sec) ( input parameter - UT) 

t] frequency parameter 

6 polar directivity angle, degrees 

o spectral shape function constant ( o = 2.2 ) 

Superscript: 

* dimensionless quantity 
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Subscript: 


a ambient 

bb broad band 

i tip speed index 

j directivity angle index 

/ one-third-octave band lower limit 

n harmonic index 

ref reference 

tone discrete tone 

u one-third-octave band upper limit 

Inlet Broadband Noise 


The prediction of inlet broadband noise is based on the correlation of fan data from the GE 
XF-120 engine. The sound pressure level (SPL) at the observer is 

SPL tb (r„B) - SPL p „ t + 10L O {r 10 (A / „,/5.584)+ S/(>/) (1) 


The broadband peak sound pressure level is determined by the following polynomial 
function of polar angle and corrected tip speed 


SPL peak - ^ 


/=! 


2 »/ 90 ) 


j~ 1 


j-i 


(u-r 


( 2 ) 


The broadband correlation coefficients A arc given in table 3.6.1, and the normalized 
corrected tip speed is defined as 


U = 


U, 


u. 


ref 


( 3 ) 


The corrected tip speed is calculated by 

TJ _ mIN 2 / 60 
JTJ 51 8.67 


(4) 


Table 3.6.1 - Broadband Noise Correlation Coefficients 



Polar Angle Index (j) 

Tip Speed Index (i) 

1 

2 

3 

1 

61.85 

- 183.92 

143.42 

2 

21.667 

705.25 

- 653.16 

3 

2.0071 

- 790.26 

731.70 

4 

- 6.2733 

272.17 

- 242.45 
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Figure 3.6.1 shows a contour plot of the peak SPL function defined by equation 2. Figure 
3.6.2 shows the peak SPL as a function of directivity angle for corrected tip speeds of 0.5, 
1.0, and 1.5. 

Scaling for airflow was not possible due to the classified nature of the XF-120 noise data. 
However, the broadband noise level is scaled for size based on rotor inlet area by the 
expression 10Log 10 (A fan / 5.584). 

Broadband Spectrum Index Function 

The broadband noise spectrum is computed by a modified version of the Heidmann 
spectral shape function [2]. The spectrum function Sir/) (figure 3.6.3) is given by 


SKrj) = 10 Log l0 



( 5 ) 


where the nondimensional frequency parameter >/ is 
?/ = (l- M fl costf)y- 

J peak 


(6) 


The peak frequency is Doppler shifted to account for the effect of source motion. The 
peak frequency is correlated as a function of tip speed and blade passing frequency. The 
equation for the peak frequency is 


/, 


peak 


2.265 + 0.07298(f7 * ) - 0.2039(f/ * ) 2 


f b 


( 7 ) 


where the blade passing frequency is 


f b = N h 



( 8 ) 


The constant o is given the value 2.2 as in the Heidmann method. The XF-120 data 
showed no significant dependence of the spectral shape function on polar angle. A plot of 
the peak frequency function is provided in figure 3.6.4. 

Inlet Tone Noise 

The correlation of the discrete tone noise levels is similar to the correlation of the 
broadband levels defined earlier. The fundamental tone and 2nd, 3rd, and 4th harmonics 
are computed as a function of polar angle and tip speed by 
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SPL. - Y. 


i = 1 


£C>/90) 


j~ 1 


7-1 


(u-r 


n = 1, 2, 3, 4 


(9) 


where n represents the harmonic number. The coefficients C fj for the first four harmonics 
are defined in tables 3.6.2 through 3.6.5. 


Table 3.6.2 - Fundamental Tone Correlation Coefficients C ij (n = 1) 



Polar Angle Index (j) 

Tip Speed Index (i) 

1 

2 

3 

1 

125.00 

-112.91 

14.02 

2 

-245.82 

688.66 

-392.05 

3 

330.86 

-883.86 

534.18 

4 

-117.90 

315.53 

-187.56 


Table 3.6.3 - 2nd Harmonic Correlation Coefficients C„ (n = 2) 

U 



Polar Angle Index (j) 

Tip Speed Index (i) 

1 

2 

3 

1 

92.291 

-59.79 

38.335 

2 

-127.82 

601.89 

-576.26 

3 

195.66 

-821.53 

751.96 

4 

-74.199 

298.90 

-258.03 

Table 3.6.4 - 3rd Harmonic Correlation Coefficients C tj 

(n = 3) 


Polar Angle Index (j) 

Tip Speed Index (i) 

1 

2 

3 

1 

135.93 

-247.51 

102.80 

2 

-185.54 

856.79 

-501.16 

3 

210.26 

-942.98 

597.45 

4 

-74.01 

322.24 

-208.53 

Table 3.6.5 - 4th Harmonic Correlation Coefficients C„ 

y 

O = 4) 


Polar Angle Index (j) 

Tip Speed Index (i) 

1 

2 

3 

1 

232.20 

-793.54 

742.31 

2 

-571.51 

2836.35 

-2612.96 

3 

621.24 

-3001.07 

2627.40 

4 

-205.75 

962.39 

-797.31 
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Equation 9 is plotted in figure 3.6.5 for the first four harmonics. Note that the 4th 
harmonic is not calculated when the tip speed is greater than 1300 ft/sec and the emission 
angle is greater than 90 degrees. The discrete tone noise levels are corrected for engine 
size, distance from source to observer, and source motion in the same manner as the 
broadband noise level. 

spk,M„8) - sPL,(u\e ) + io lo S Ja^ /5.584> (io) 

The frequency associated with each harmonic is obtained by Doppler shifting the blade 
passing frequency and multiplying by the appropriate harmonic number. 


fn = 


n f h 

1 - M „ cos# 


n = 1, 2, 3, 4 


( 11 ) 


2D Bifurcated Inlet Correction 

Inlet shape effects are applied using correlation equations based on HSCT inlet data taken 
in Boeing's low speed aeroacoustic facility (LSAF) in 1995. These equations detail 
suppression differences between the 2D bifurcated inlet and typical subsonic inlets. 

Deltas from these correlations are applied spectrally and include directivity effects. The 
broadband correction is expressed as a polynomial function of frequency and emission 
angle. 


A SPL hb = a 0 + a x f + a 2 f 2 + a 3 d + a 4 6 2 + a 5 fd (12) 

Table 3.6.6 defines the coefficients for approach, sideline, and centerline microphone 
positions. 


Table 3.6.6 Broadband Coefficients for 2D Bifurcated Inlet Correction 



a 0 

ai 

a 2 

a3 

a 4 

as 

Sideline 

4.26E+00 

3.86E-04 

-1.10E-07 

-7.46E-02 

2.30E-04 

1.57E-05 

Approach 

1.99E+00 

-5.60E-04 

-6.65E-08 

7.63E-02 

-9.76E-04 

-4.99E-06 

Centerline 

7.34E+00 

-5.35E-04 

4.75E-08 

-1.05E-01 

1.50E-04 

-7.25E-06 

(Doors Closed) 







Centerline 

3.76E+00 

1.28E-03 

-1.44E-07 

-1.08E-01 

5.49E-04 

1.67E-06 

(Doors Open) 








The tone correction is expressed as a polynomial function of the emission angle. 

A SPL tone =c 0+Cl 6 + c 2 6 2 (13) 

Tables 3.6.7 through 3.6.10 define the coefficients for the first four fan harmonics. 
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Table 3.6.7 1st Harmonic Coefficient for 2D Bifurcated Inlet Correction 



Co 

Cl 

C2 

Sideline 

-5.34E+00 

1.52E-01 

-6.06E-04 

Approach 

-1.40E+01 

4.25E-01 

-2.82E-03 

Centerline (Doors Closed) 

0.0 

0.0 

0.0 

Centerline (Doors Open) 

0.0 

0.0 

0.0 


Table 3.6.8 2nd Harmonic Coefficient for 2D Bifurcated Inlet Correction 



Co 

Cl 

c 2 

Sideline 

-8.56E+00 

3.48E-01 

-3.36E-03 

Approach 

-3.90E+01 

1.21E+00 

-1.01E-02 

Centerline (Doors Closed) 

0.0 

0.0 

0.0 

Centerline (Doors Open) 

0.0 

0.0 

0.0 


Table 3.6.9 3rd Harmonic Coefficient for 2D Bifurcated Inlet Correction 



Co 

Cl 

C2 

Sideline 

6.17E+00 

-1.26E-01 

1.06E-05 

Approach 

-5.40E-01 

-1.89E-01 

0.0 

Centerline (Doors Closed) 

0.0 

0.0 

0.0 

Centerline (Doors Open) 

0.0 

0.0 

0.0 


Table 3.6.10 4th Harmonic Coefficient for 2D Bifurcated Inlet Correction 



c 0 

Cl 

c 2 

Sideline 

6.09E+00 

-1.31E-01 9.11E-04 

Approach 

1.61E+00 

-8.57E-02 6.25E-04 

Centerline (Doors 

0.0 

0.0 

0.0 

Closed) 

Centerline (Doors Open) 

0.0 

0.0 

0.0 


Total Noise 

The total mean-square acoustic pressure is the sum of the broadband component and the 
discrete tones 


(p 2 ) = (p ) +(p ) 

/ total SI / bb M' /tone 


where 


[p) hh =Plf m 


(SPL bb + ASPL bb )llO 


(14) 

(15) 
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and 


= Pref 


10 


( SPL lone +ASPL lone ) / 10 


(16) 


Since broadband noise is computed at the one-third-octave band center frequency and tone 
noise occurs at discrete frequencies, it is necessary to identify the one-third-octave band 
which contains the tone. For a given center frequency, / , the lower limit of the band is 


f = 2~ 1/6 / 

and the upper limit is 

f u =T ll6 f. 


(17) 


(18) 


The tone mean-square acoustic pressure at each harmonic number, n , is then added to the 
band in which the relation f ,< f < f u is true. 

The free-field lossless total inlet radiated fan sound pressure level is 

SPL,„Jr,.0)-lOU,g lo {p\ m , (19) 

Propagation effects are defined in section 3.13. 
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SPL , ,dB SPL ,,dB 








3.7 Combustor Noise 


The combustion noise prediction method adopted for the HSR program is based on the 
method developed by Mathews and Rekos [3]. This method was developed from several 
JT8D burner configurations and full-scale turbofan engine data. The full-scale engines 
included the JT8D-109, JT9D-7A, JT9D-70, and the JT10D. Comparisons of predictions 
with data from both the turbine rigs and the full-scale engines showed that the equation for 
the OAPWL fits the data with a standard deviation of 1.7 dB. Reference 3 provided no 
comparisons with flight data. 


Symbols 


B, 


Di(e ) 

f 

F 


F b 


f, 


peak 


ICR 

K f 

(. L/jrD ) 
OAPWL 

Pref 

2 

< P > 

P,3 


R 

Slip) 

SPL 

TL 


t4 

T 
1 a 


2 

burner cross-sectional area, (ft") (input parameter - AB) 
burner length, (ft) (input parameter - BL) 
specific heat at constant pressure = 0.28 BTU/lbm °R 
directivity index function 
frequency, Hz 

ratio of acoustic impedances across turbine 
burner fuel to air ratio (input parameter - FB) 
combustion noise peak frequency, Hz 
Stochiometric fuel to air ratio = 0.068 
fuel heating value 18,500 BTU/lbm 
burner type (input parameter - ICR) 
coefficient in peak frequency equation 

transmission loss through burner exit duct (input parameter - PLD) 
number of fuel nozzles (input parameter - NFN) 
overall power level, dB 

reference pressure, 4.177 x 10" 7 (lb/ft 2 ) 

2 2 

mean-square acoustic pressure, (lb/ft") 

burner inlet total pressure, (lb/ft 2 ) (input parameter - P3) 

burner exit total pressure, (lb/ft") (input parameter - P4) 

low “pressure” turbine exit total pressure, (lb/ft") (input parameter -PI) 

distance from source to observer, (ft) 

gas constant = 53.35 ft-lbf/ lbm °R 

spectrum function 

sound pressure level, dB 

combustion noise transmission loss, dB 

burner inlet total temperature, (°R) (input parameter - T3) 

burner exit total temperature, (°R) (input parameter - T4) 

low “pressure” turbine exit total temperature, (°R) (input parameter - T7) 

burner mass air flow rate, (lbm/sec) (input parameter - W3) 
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(w f / P t4 fuel-flow-rate to bumer-exit-pressure ratio at takeoff, 

(lbm/sec / lb/ft 2 ) ( input parameter — WFP 4 ) 
i] nondimensional frequency parameter 

6 polar directivity angle, degrees 

Method 

The combustion noise prediction method adopted for the HSR program is based on the 
method developed by Mathews and Rekos [3]. The free-field lossless sound pressure 
level is 

SPL(r s ,6) = OAPWL -TL + DI(0) + SI(rj\ ) (1) 

where the overall acoustic power is 

. f t' U T? \ 2 

1 2 w, ,//’ H ft 

OAPWL = 10 Log 10 A: PA ^ ,4 1 + ^^ F; +131.3 (2) 

MPA r T 

iy f y r t3 A b ) y c P 1 t3 

Propagation effects are defined in section 3.13. 

Equation 2 was derived from Strahle’s equation for sound produced by open flames in 
terms of burner geometry and performance parameters. The acoustic power is directly 
proportional to the burner air flow rate, the fuel-to-air ratio, and the total temperature at 
the burner exit. Increasing either the air flow rate or the fuel to air ratio will increase the 
acoustic power. The acoustic power is inversely proportional to the number of fuel 
nozzles, the cross sectional area, and the burner inlet total pressure and temperature. 
Increasing either the number of fuel nozzles or the burner cross sectional area will reduce 
the acoustic power. 

Transmission Losses 

Transmission losses are caused by the loss of energy through the duct leading into the 
turbine and across the turbine rotor and stator blades. The transmission loss factor is 
defined by two expressions. The equation for the total transmission loss is 

TL = 10Log 10 [(l + f) 2 I AF(LI jtD) (3) 

The transmission loss associated with the turbomachinery is given by 

F-(PJPn)FiFF (+ 

Equation 4 was derived by modeling the turbine as a discontinuity in characteristic 
impedance. The ratio of the characteristic impedance across the turbine is a function of 
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the total pressures and total temperatures at the combustor and turbine exit locations. An 
increase in combustor exit total pressure or turbine exit total temperature increases the 
transmission losses and reduces the combustor noise. An increase in combustor exit total 
temperature or turbine exit total pressure reduces the transmission losses and increases the 
combustor noise. The expression Lt jzD is associated with losses in the duct and has a 
value of 0.15 for the mixed flow turbofan. Note that the transmission loss through the 
duct can be ignored by setting Lt jiD = 1 . 


Directivity Function 

The directivity function, defined in table 3.7.1, was established for observers located at a 
distance of 150 feet from the source. The directivity function is plotted in figure 3.7.1. 


Table 3.7.1 - Directivity Function 


6 

Dl{0) 

0 

DI^) 

0.0 

-60.0 

95.0 

-42.9 

10.0 

-58.3 

100.0 

-42.0 

20.0 

-56.5 

105.0 

-41.2 

30.0 

-54.8 

110.0 

-40.4 

40.0 

-53.0 

115.0 

-39.8 

50.0 

-51.1 

120.0 

-39.2 

60.0 

-49.3 

130.0 

-40.9 

70.0 

-47.5 

135.0 

-42.6 

80.0 

-45.7 

140.0 

-44.3 

90.0 

-43.8 

150.0 

-48.0 



180.0 

-59.0 


Spectrum Function 


The spectrum function was derived from a combination of combustor rig and full scale 
engine data. The peak frequency is a function of the burner geometry and burner type. 
The expression for the peak frequency is 


f peak 


RH , 


'V 


V / re j 


A,B 


L b"i 


(5) 


where K f = 8 for can type burners and K f = 3 for annular type burners. Equation 5 

indicates that can type burners have higher peak frequencies than annular burners of 
comparable size. The spectrum function in table 3.7.2 is tabulated as a function of the 
nondimensional parameter 


V = 10 Log l0 


\ f peak J 


(6) 
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The spectrum function is plotted in figure 3.7.2. Note that the spectrum is not Doppler 
shifted. 


Table 3.7.2 - Spectrum Function 


V 


r-~ 

i 

V 

Sr- 

S/( //)= -24.7 + 2.5(?/ + 7) 

-7 

-24.7 

-6 

-22.3 

-5 

-19.5 

-4 

-16.8 

-3 

-14.1 

-2 

-11.4 

-1 

-8.8 

0 

-6.8 

1 

-8.1 

2 

-9.9 

3 

-11.5 

4 

-13.3 

5 

-15.0 

6 

-16.8 

7 

-18.6 

8 

-20.5 

9 

-22.4 

10 

-24.3 

11 

-26.3 

?]> 11 

Sl{ v) = -26.3 - 2 . 0(77 - 1 1) 



3.8 Turbine Noise 


The turbine noise prediction method adopted for the HSR program is based on the method 
recommended by Krejsa and Valerino [4]. The empirical correlation was developed from 
turbine rig and full scale engine data. The uncertainty in predicted levels can be as high as 
+/- 5 dB for the fundamental tone and +/- 9 dB for the broadband noise when compared 
against data from commercial turbofan engines. The large prediction uncertainty is due to 
the wide range in turbine designs, geometry, and the complex interaction with noise from 
upstream turbine stages. An evaluation of this method was made with data from the GE 
XF120 engine [5]. These comparisons indicated that the trend with increasing turbine 
speed was not well predicted and that the model would over predict HSCT turbine noise 
by as much as 10 dB. Consequently, it is suggested that the turbine calibration constant 
TURCAL be set to - 10 dB . 

Symbols 

c ref reference sound speed, c re j = 11 16.45 ft/ sec 

c tex sound speed at the turbine exit, ft/sec 

CIS stator axial chord / stator rotor spacing ( input parameter - CSS) 

d last stage turbine rotor diameter, ft 

D1 directivity index, dB 

/ one-third-octave band center frequencies, Hz 

f hpf blade passing frequency, Hz 

2 

g c gravitational constant, 32.174ft/sec 

m primary mass flow rate at turbine exit, lbm/sec ( input parameter - MDOT) 

m ref reference mass flow rate at the turbine exit, m re j = lbm/sec 

N b number of turbine blades in the last stage ( input parameter - NBLADE) 

N 2 rotational speed, rpm ( input parameter - TURRPM) 

r s distance from source to observer, r s = 150 ft 

R gas constant, 53.35 ft-lbf / lbm °R 

SI spectrum index, dB 

T s turbine exit static temperature, °R ( input parameter - TS) 

TURCAL turbine calibration constant, dB ( input parameter - TURCAL) 

V ref reference tip speed, V re j = ft/sec 

V t tip wheel speed, ft/sec 

V tr relative tip speed, ft/sec ( input parameter - VTR) 

y tex ratio of specific heats at the turbine exit 

j] nondimensional frequency parameter 

6 emission angle, degrees 
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Method 


The free-field sound pressure levels for the broadband and discrete tone noise at a radius 
of 150 feet are given by equations 1 and 2. 


SPL hh {r s =150, f,8) = lOLog 10 


( v c ^ 

V tr C ref 1 

3 

( . \ 


m 

v r 

y ref tex J 


l <f) 


+ DI bb (d) + SI bb (rj)~ 10 


( 1 ) 


SPL tone (r s =150, f bpf ,d)=lOLog lo 


Yjl 
V V ref J 


'ref 


V C tex J 


\ 


m 


V m ref J 




DI,„M 


( 2 ) 


+ W«(/J+ 56 


The free-field lossless total turbine noise is obtained by combining the broadband and 
discrete tone noise components and removing the atmospheric absorption effects. 

SPfo,al(r„f,0) - 10Los 10 (l0° ls "» + 10° lspL >one ) + TURCAL (3) 


The function F abs is corrects the total noise to a lossless condition by the removal of the 
standard day (i.e. 59 °F, 70% relative humidity) atmospheric attenuation at a distance of 
150 feet as defined in reference 6. Note that turbine calibration constant TURCAL allows 
the user to choose the amount of noise reduction to be applied to the total turbine noise. 
TURCAL was determined to be -10 dB for the HSR program based on data from the GE 
XF120 engine. Propagation effects are defined in section 3.13. 

It is assumed that the design of the HSR turbine rotor blades will remove most or all of the 
swirl. To the accuracy of the method, this allows the relative tip flow velocity to be 
calculated (by the user) as 

v; = v 2 + v 2 (4) 


where the tip speed of the last stage rotor is 

T/ JtdN 2 
f 60 


(5) 


The axial flow velocity is given by 
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(6) 


v: = 


r, e , - 1 


~,T, . 


where the ratio of specific heat is 


0.4 

1 + 0.4 

_2 T 

r e 



( 7 ) 


and r represents the temperature ratio defined by 

r = 5000/r s (8) 

The speed of sound at the turbine exit is 

= a / Y,e X Sc RT s ( 9 ) 

Stator Rotor Spacing 

The stator rotor spacing is defined in figure 3.8.1 and usually has a value greater than 1. 
Directivity Function 

The broadband and discrete tone directivity functions are tabulated in table 3.8.1. 


Table 3.8.1. Broadband and Discrete Tone Directivity Functions 


0 

DI bb ( e ), dB 

Dltone(0), dB 

0 

- 37.0 

- 56.0 

40 

- 21.0 

- 27.5 

60 

- 12.9 

- 18.0 

80 

- 7.0 

- 10.0 

100 

- 1.3 

- 2.4 

110 

0.0 

0.0 

120 

- 1.3 

- 2.2 

140 

- 8.8 

- 14.0 

160 

- 19.0 

- 26.0 

180 

- 31.5 

- 40.0 
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Spectrum Index Function 

The variation in broadband noise as a function of frequency is given by the spectrum 
index function. The broadband spectra peaks at the blade passing frequency as indicated 
by equation 10. The discrete tone spectrum decays 10 dB for each harmonic greater than 
the blade passing frequency. 


siM- 


\QLog H) tj i] < 1 
20A16Log w ?] ?/ > 1 


( 10 ) 


SI".— 10(>7-l) t, -1,2,3... 


( 11 ) 


The spectrum index function is defined in terms of the nondimensional frequency 
parameter >j where rj is the ratio of the observed one-third-octave band center frequency 
to the Doppler shifted blade passing frequency. 

rj—f- ( 12 ) 

Jbpf 


f hpf = N b 


AT, 


60(l - M a cos d) 


(13) 


Table 3.8.2 defines the atmospheric absorption in decibels that must be added to correct 
the predicted levels to a lossless condition. The atmospheric absorption levels are 
calculated for a distance of 150 feet using reference 6 for standard day conditions defined 
as 59 °F and 70% relative humidity. 
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Table 3.8.2 Standard Day Atmospheric Absorption 


One-third-octave Band 
Center Frequency, Hz 

Fabs(f), 

dB 

50 

0.02 

63 

0.02 

80 

0.02 

100 

0.02 

125 

0.03 

160 

0.03 

200 

0.05 

250 

0.06 

315 

0.08 

400 

0.09 

500 

0.11 

630 

0.14 

800 

0.18 

1000 

0.23 

1250 

0.29 

1600 

0.36 

2000 

0.45 

2500 

0.60 

3150 

0.81 

4000 

1.14 

5000 

1.37 

6300 

1.91 

8000 

2.78 

10000 

4.11 
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Stator Vanes Rotor Blades 


Axial Flow 
Direction 



C Stator Axial Chord 

css - — 

S Stator Rotor Spacing 

Figure 3.8.1 Definition of stator chord to stator rotor spacing ratio 
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SPL, dB SPL, dB 



Figure 3.8.2 Turbine broadband and discrete tone directivity functions 



Figure 3.8.3 Turbine broadband spectrum function 
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3.9 Jet Noise (Boeing JN8C4 Methodology) 


Adoption of the mixer-ejector as a means of reducing jet noise necessitated the 
development of new jet noise prediction methods. Research was conducted on a number 
of mixer-ejector configurations during the HSR program for both circular and rectangular 
configuration. As data became available, empirical jet noise prediction methods were 
developed at Boeing with the designation JN8. The method described in this section is 
the JN8C4 version developed by Lu and Viswanathan [7,8,9] for a supersonic jet issuing 
from a multilobe mixer into a rectangular ejector. 

The jet is conceptually divided into three regions, the primary internal jet, the primary 
external jet, and the mixed jet. High frequency noise, termed “primary internal noise”, is 
generated by the turbulent mixing of the supersonic primary flow with the entrained air. 
Shocks may be present in this region as a result of non-ideal expansion of the supersonic 
primary flow. Mid-frequency noise, termed “primary external or pre-merged noise”, is 
produced by a combination of two mechanisms; (a) the continued mixing of the primary 
and the entrained flows inside the ejector and immediately downstream of the ejector 
exit, and (b) the initial mixing of the engine flow with the freestream. If the mixed jet is 
supersonic, the shock-associated noise of the mixed jet will also appear in the mid- 
frequency region. Low frequency noise, termed “mixed or merged noise”, is generated 
by the mixing of the “mixed” jet (the initially rectangular jet having decayed to a near 
circular shape) with the freestream. The peak noise source is several equivalent 
diameters downstream of the ejector exit, where the length scale is quite large, usually on 
the order of the ejector equivalent diameter. Figure 3.9.1 shows a graphical 
representation of the three mixing noise components in the physical and spectral 
domains. The procedure predicts the one-third-octave band sound pressure level (SPL) 
for each component. The total jet noise SPL is the logarithmic sum of the five 
components: 

1. primary internal noise 

2. primary external or pre-merged noise 

3. mixed or merged noise 

4. broadband shock noise from the primary jet 

5. broadband shock noise from the mixed jet 

For certain operating conditions, the contribution of broadband shock noise to the total 
noise produced by the mixed ejector could be significant. The shock noise prediction 
procedure incorporated into this module is based on the instability wave/shock cell 
interaction source model developed by Tam [10,11], The prediction of broadband shock 
noise in the internal jet region follows Tam’s model for circular jets. Empirical 
adjustment is made for the rectangular geometry, liner effects, and frequency. Shocks in 
the primary external jet are predicted by Tam model for rectangular nozzles [12]. 
Empirical adjustments are made for flight effects. 
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Symbols 

Aibs 

A 

A 

A T 


A 

b 

K 

c a 

C m 

C j 

C 

C fg 

D ej 


D 


D, 


plug 


D 2 

d 

Ip 

^ 'asp 
^ asp 1 

F 

FV t 

f 

flin 

L 

A 

G 

A 

A 

A 

K 

L 


lin 
l div 

M a 

Ml 


atmospheric attenuation coefficient, dB/ft 

2 

primary nozzle exit area, ft” ( input parameter - APRI) 

secondary nozzle exit area, ft' (input parameter - ASEC) 

throat area of the primary nozzle, ft' (input parameter- ATHP j 

area of a single lobe, ft' 

shock cell strength function 

source length, ft 

width of nozzle at exit plane, ft 

ambient speed of sound, ft/sec 
wave number parameter 
speed of sound in secondary jet, ft/sec 
broadband shock constant 

mixer-ejector thrust coefficient (default =0.95) (input parameter - CFG) 
equivalent inner diameter of the ejector at the mixer exit, ft 

( input parameter - E JD) 

Mixed jet nozzle exit diameter, ft 

primary nozzle exit diameter, ft (input parameter - PLUGD) 
diameter of a single lobe, ft 
d,j, =(*P -ip 0 )/(0.5jr) 

mixer-ejector aspect ratio major/minor axis (input parameter - E JASP) 

mixer-ejector aspect ratio minus one ( E - 1) 

spectral distribution function 

Source strength function 

one-third-octave band center frequency, Hz 

liner design center frequency, Hz (input parameter - FLIN) 

Doppler shifted one-third-octave band center frequency, Hz 

/ L — 

p 1 - M a cos 6 p 

ground effects, dB 

height of nozzle at exit plane, ft 

fully expanded nozzle height, ft 

first order Bessel function 

wave number 

shock cell spacing 

ejector length, ft (input parameter - E JL) 

length of ejector lining, ft (input parameter - E JLIN) 

exit divergence angle, degrees (input parameter - DIVANG) 

aircraft Mach number ( input parameter - MA) 

special flight effect aircraft Mach number 
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M c 

M d 

Mj 

N en g 

N, 

n , 

P a 

P en 



P 

r s 

rh 

R. 


R, 


iin 


S i 

^ lin 

T. 


T s 

K 

V c 

V m 

V P 

y P e 

K 

Wp 

W s 

X, 

yobs 


Pi 

X 

<P 

Pa 

Pi 


convection Mach number 
design Mach number 
fully expanded jet Mach number 
number of engines 

number of spokes or lobes of the mixer (input parameter - SPOKE) 
velocity exponent 

ambient pressure, lb/ft" (input parameter - PAMB) 
mixer spoke penetration to nozzle full height ratio, 

(input parameter - PEN) 
standard day pressure, lb/ft' 

mean-square pressure, (lb/ft 2 ) 2 

distance from apparent source (in the mixed jet ) to observer, ft 
distance from apparent source (in the primary jet ) to observer, ft 
distance from nozzle exit to observer at the time of emission, ft 
relative humidity, % 

ratio of ejector length to minimum inner diameter of the ejector 
ratio of ejector lining length to minimum inner diameter of the ejector 
Strouhal number 
liner Strouhal number 

ambient temperature at microphone location, (°R) (input parameter - TAMB) 

primary jet total temperature, (°R) (input parameter - TPRI) 

secondary jet total temperature, (°R) (input parameter - TSEC) 

aircraft velocity, ft/sec (input parameter - VAIR) 

convection velocity, ft/sec 

mixed jet velocity, ft/sec 

primary jet fully expanded velocity, ft/sec (input parameter - VPRI) 
external velocity, ft/sec 

secondary jet velocity, ft/sec (input parameter - VS EC) 

primary jet mass flow rate, lbm/sec (input parameter - WPRI ) 

secondary jet mass flow rate, lbm/sec (input parameter - WSEC) 

noise source location down stream of the primary of common exit, ft 

observer sideline distance from flight path, ft (input parameter - YOBS) 

observer altitude, ft 

source altitude, ft 

pressure ratio parameter 

Mach number ratio 

azimuth angle, degrees 

ambient density, lbm/ft 

jet density, lbm/ft 
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r modal summation parameter 

y ratio of specific heats 

y> azimuth angle measured from the ejector major axis, radians 

ip {) rectangular exit plan major axis roll angle relative to wing plane, degrees 

( input parameter - PS 10) 

II acoustic power 

d ohs observer angle measured from inlet, degrees 

<9, wave normal angle measured from inlet, radians 

o m zeros of the zero-order Bessel function 

£ 1 -Mlsm(ip) 

Subscripts 


i general index (i = p,pi,pe,orm) 

p primary 

pe primary external 

pi primary internal 

m mixed 

s secondary 
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Component SPL Formulation 

The general formula for calculating the one-third-octave band SPL for each jet mixing 
noise component is given by 

SPL,(f,e)-[{zi)lM g (FV,) + (Z2 l )\lM S S l -(Z3 l )lMg(FVy(Z4 l )\ 

+ (Z5, )Log(FV j ) + (Z6, ) 

Formulas for the source strength FV i , the Strouhal number Sj. , and the coefficients Zl ( . 
through Z6, are provided in following sections. The subscript i is a general index, which 
indicates the jet component. Here i = pi, pe, or m to indicate the primary internal, 
primary external, or the mixed jet component. When the primary internal and the 
primary external components share calculations, the index i = p is used. 

The source strength function FV t is correlated as a function of nondimensional shear 
layer velocity difference, eddy convection velocity, and for the primary jet, the Mach 
number and the ambient flow effects. The velocity exponent provides the scaling of SPL, 
with the jet and ambient flow acoustic Mach numbers (nondimensional velocities). The 
velocity exponent has a strong directivity dependence and a weak Strouhal number 
dependence. The Strouhal number is defined for each component, using the shear layer 
velocity difference, the component jet diameter and source frequency. 

Figure 3.9.2 shows graphically how each coefficient affects the SPL. Coefficients Z1 1 
and Z2 ( determine the spectrum distribution. Coefficients Z3 ( . and Z4 ( determine the 
location of the peak frequency. Coefficients Z5, and Z6 ( determine the SPL amplitude. 
All the coefficients are more directivity dependent than Strouhal number dependent. 

Coefficients Zl ( through Z6 ( are functions of the source radiation angle and other 
parameters. The mixing noise components are modeled as a distribution of point sources. 
Each source has a unique radiation angle. Calculation of the radiation angle from the 
distributed source location is explained in the following sections. 

Jet Noise Source Distributions 

Noise sources are distributed over a wide range of the jet mixing region, which extends 
from the nozzle exit to a distance of 15 mixed jet diameters downstream. Each source has 
a radiation angle that differs from the geometric angle associated with the position of the 
aircraft. The location of the source within the mixing region, the radiation distance, the 
source radiation angle, and the source frequency are uniquely determined for each source 
and observer position. This section describes the procedure for locating the sources in 
the primary internal jet, the primary external jet, and the mixed jet. 
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Source Distribution for the Primary Internal and Primary External Jets 
The source distribution for the primary internal and external jets is given by 





(\W n ^ 

A 

4 + 4tan 1 

p 9 

+ 2X 


K JT ) 



D 


V 1 + A ^ 


( 2 ) 


Sources radiating in the forward arc are located close to the nozzle exit, while sources 
radiating in the aft arc are located further downstream of the nozzle exit. Note that the 
source distribution formula for the primary jets is independent of frequency. Frequency 
dependent forms were tested but did not improve agreement with measured data. 

Source Distribution for the Mixed Jets 

The mixed (merged) jet source location normalized by mixed jet diameter is given 
empirically by 


X 


m 


3 + exp(- S m )- 


2 + l.ltan 


-i 


\ Jt J 


1 + 0.55. 


f V ) 

0 5 + 05'" 

/2 

f V. 1 

1 D m 

U • 1 U • 


{v m -vj 

+ 

O 

o 

+ 


(3) 


Although the mixed jet does not have a potential core, the characteristics of source 
distributions are similar to that of the primary internal and primary external jets. The 
source distribution function for the mixed jet is a function of frequency. 

Determination of Radiation Angle and Source Frequency 

Since noise sources are distributed over a wide range of the jet mixing region, the 
geometric angle, which identifies the observer location with respect to the nozzle exit, 
differs from the source radiation angle. The source location and radiation angle 
associated with a given observer is determined by an iterative procedure described below. 

Static or Wind Tunnel Coordinates 

When there is no relative motion between the observer and the nozzle, 

(figure 3.9.3a) the source for a given component is located at X and radiates sound at an 
angle d i (wave normal angle). The subscript i = p or m and the angles (J are in radians. 
The observer (microphone) angle is 6 obs at sideline distance y ohs as shown in 
figure 3.9.3b. The source location and radiation angle are determined by the following 
iterative procedure. 
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Step 1. For a given one-third-octave band frequency (Strouhal number Sj), initialize 

9 t = jt / 2 and compute the source location X i (dimensional) for the desired region 
from equations 2 or 3. 

Step 2. Compute a new radiation angle 6 j new using 


V,,ne K = sin 


-1 


- M B + JB 2 - M: + 1 


B 2 + 1 


where 
B = (X i /SL)+ cot G t 


obs 


M„ 


< 1 


X. is the source location from step 1 , SL is the sideline distance, and 9 ohs is the 
angle for the nozzle exit to the observer measured from the inlet axis. Compute a 
new source location X ( . using 6 i = 6 i new . 


Step 3a. Let 9, = 9 inew . 

Step 3b. Compute a new radiation angle 9 inew using X ; .. 

Step 3c. Let d 2 = 9 i new obtained from step 3b. 

Step 4a. Let 6 j = {d x + d 2 )/ 2 . 

Step 4b. Compute a new source location X. using 6 j = (0 ] + 0 2 )/ 2 . 
Step 4c. Let X i new = X ; . . 

Step 5. Repeat steps (3) and (4) until |X ( - X ( new | < 0.05 Z>. 


Transformation from Wind Tunnel to Ground-Fixed Coordinates 

The iteration procedure for ground-fixed coordinates is the same as the wind tunnel case 
except, the source radiation angle (wave normal) is given by 


9: 


sin 


1 

B 2 + 1 


(4) 


Calculation of the Mixed Jet Parameters for the Aspirated Mixer-Ejector 

Calculation of the mixed jet noise component level requires mixed jet parameters as well 
as some primary, secondary, and the ambient flow parameters. The mixed jet parameters 
are not independent parameters, since they are calculated from the primary and the 
secondary jet parameters under the fully-mixed assumption. For most mixer-ejectors, the 
“secondary” jet is the aspirated (induced or entrained) ambient flow. To calculate the 
mixed jet velocity, it is assumed that the total momentum of the primary stream and the 
aspirated stream, before and after mixing, is conserved (i.e., the net thrust is independent 
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of aspiration). Since the aspirated flow has the ambient stream momentum, calculation of 
the mixed jet velocity V m for the aspirated mixer-ejector is given by: 


V,., 


V W +VW 

p p a s 

W + W 

p s 


( 5 ) 


where the ambient (wind tunnel or airplane) velocity V a is used. Other conditions of the 
mixed jet are calculated by using conservation of mass flow and total enthalpy. 


T _ T _T p W p+ T,W, 


' m t.m 


w p +w. 


w p +w s 

A - p 

m T j 

P„y,n 


(6) 

( 7 ) 


D 


m 



( 8 ) 


The mixed jet density in addition to other jet parameters is determined by a two step procedure. 
The ratios of specific heats for the primary, secondary, and mixed jets are calculated initially 
using the total temperatures as follows 


y = 1.357 + 0.058 tan" 


y s = 1.368 + 0.050 tan" 


y m = 1.363 + 0.054 tan" 


1200 -T, 


t,S 


v 870 J 


1200 -T, 


t,s 


870 

1200 -T. 


t,m 


830 


The static temperature are computed by 


Y - 1 

r r i y ' P 

1 s,p ~ t,p ~ ^ 


T =T - 

S,S t,S 


Ys "I 


7 y 2 3 

p 

0', R J 

7 y 2 3 


\Ys R j 


T = T - 

s,m t,m 


Y - 1 

/ m 


7 y 2 3 

m 

\Y m R j 


( 9 ) 

( 10 ) 
( 11 ) 


( 12 ) 

(13) 

(14) 
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Finally, the static temperature for the primary, secondary, and mixed jets is determined by 



( 21 ) 

( 22 ) 

( 23 ) 


From the quantities computed above, the fully expanded jet Mach number, jet densities, and 
pressure ratio parameters can be computed. 
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Fully expanded jet Mach number parameters 


M J , = 


y P - 1 


— ^ — i 
T 

\ s ’P J 


M us = 


Ys - 1 


T 

-^-1 

T 

v ^ y 


M i,m = 


y - 1 

/ m 


T 
v J 


-1 


Jet density parameters 


Pp RT 


s,p 


Ps = 




Pm 


RT 


(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


Pressure ratio parameters 


Pp = 


(30) 

fis 


(31) 

fh = 

iK-' 1 

(32) 
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Primary Internal Component of Jet Mixing Noise 

The high frequency noise, termed “primary internal noise”, is generated by the turbulent 
mixing of the supersonic primary flow with the entrained air. The length scales in this 
regime are small with the macro scales being on the order of the lobe width. The sound 
pressure level for the primary internal jet mixing noise is 

SPL pl (f,e^,r, ) - [(Zl pl )L B gjFV r ,) + (z2„ - {z3 jLog m (FV p , )- (z4„. J 

+ (Z5„ )Log m (FvJ+ (Z6„ )+ A SPL„ 

The last term in the equation above represents propagation effects which are discussed 
later in this section. Formulas for the source strength FV pi , the Strouhal number S pi , and 

the coefficients Zl pi through Z6 pi are provided below. 


Source Strength Function 


Primary Internal Jet 


FV = 

pi 


V -V 

p s 


,0.6 


'a J 


V + V 

p s 


, 0.4 


'a J 


(34) 


Strouhal Number 


Primary Internal Jet 


p> 



(35) 


D 


p 



(36) 


Coefficient Z1 

Primary Internal Jet 


( 1.8(9 V 

Z1 " __18 t f- 0 - 6 ) 




(37) 

Coefficient Z2 

Primary Internal Jet 




(1.80,, 3 

2 / \\ 

f 1 - W p 1 

23 - 20 

p 0.6 

+ 5(l + Log l0 S pj K 1 + 

-- 0 ’ 6 


l 71 ) 

l n J) 


+ C2l(R ej ,0 p ) (38) 


Adjustment for ejector length to inner diameter ratio 
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Z4 ; , = O.GLogJl + iV / ){l - O.2exp[-5(0, -l.s) 2 }-0.6 

W l l-86>„ 'j 

-0.2 — - +0.l(/? +l) — - 0.6 +0.5 

TIT- P ' 


V -V (1.80 
— - 0.6 


+ C4\(P en ) + C42(E asp ,0 p ) + C43(R ej , 6 p ) 

Adjustment for lobe penetration 
C41 = 0.2 P en 

Adjustment for aspect ratio 


f 

C42 = +0.1E 1 — - 0.6 

JT 


Adjustment for ejector length to inner diameter ratio 

(\ 80 A 4 

C43 = -0.4Log 10 (l + R ej ) 1- — — — -0.6 

TC , 


a /Y ~ 1 for E aspl > 0 
E asp\ = i / X mar ^ 

0 for E asp t < 0 
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Coefficient Z5 


Primary Internal Jet 


Z5 pi = 50 + 20exp[- (ff p - 2.6) 2 _ 


Coefficient Z6 


Primary Internal Jet 


Z6 pi = 142 - 56^^ - 0.6j + 46exp[- (o p - 2.5) 2 ]- 8Log 10 (l + N, )- 5^- 

+ 2exp - lo[d p - 2.s) 2 ] + 5exp[- 5^6 p - 1.9) 2 (49) 

+ ^6 1 (E asp , E aspl ,</,) + C62(C fg ) + C63(M; , d ¥ ) + C6A{R ej , R lin ,S p , S lin ) 

Adjustment for aspect ratio and azimuth angle 


C61 - I 1 + 6exp[- lofc, - 2.6) 2 }- 6 E a 


+ 5 MC - 8 (l - OMl } . 0.6 


+ 2 E t „ l , ex s>r w ( e p ~ 2 - 5 ) 2 . +5 V E OSJ ; 1 exp-5(e j; -1.4) 2 _i-o. 5 rfy) 

Adjustment for thrust loss 
C62 = -40(l - C /g ) 

Adjustment for flight effects and azimuth angle 

C63-+3VMy{exp[-10(tf„ -l) 2 ]-exp[-lo(tf, -2.4> = } 

+ 15 1 /MT<,exp(-5 J ff p )(xp[-10(^ -23) 2 ]-0.8exp[-2o((? p -1.5) 2 } 


Adjustment for ejector length to inner diameter ratio, ejector liner, and Strouhal number 
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C 64 = — lo[l - exp(- R ej )] | - 0.5exp[- 5^ - 1.6) 2 ] + 0.4exp[- 5 ( 0 p - 2.2 ) 2 ] (53) 

- 0-3/? 2 exp[- 5{d p - 1) 2 ] 

2.5 - - °- 6 j Run [2 - exp (R, in - R ej )]exp[- 0.5(Log H) S pi - Log l0 S lin ) 2 _ 

- 1 - 3 {l + 0.4exp[- 5 ( 0 p - 1.6) 2 }/?„,, [2 - cxpf/?,., - R ej )]exp[- 0.5 (lo <?10 5 ; „. - Log m S Un 


v 0.5 ji 


(54) 


, \M for M < 0.3 
M'=\ 

a [0.3 for M > 0.3 


(55) 


/? — ^lin 

^ in ~ 


(56) 


(/», +0.0l)D 


S »" (V -V,) 


(57) 


Primary External Component of Jet Mixing Noise 

The mid-frequency noise, termed “primary external or pre-merged noise”, is produced by 
a combination of two mechanisms; (a) the continued mixing of the primary and the 
entrained flows inside the ejector and immediately downstream of the ejector exit, and (b) 
the initial mixing of the engine flow with the freestream. The sound pressure level for the 
primary external jet mixing noise is 

SPLjf,e,0,r, ) - [(zi Jlog(FV„)+ (z2„)IlogS„ - (z3„ )log(FvJ- (z4„ )f 

+ (Z5,„ )log(FvJ+ (z6 pp )+ A SPL p , 

The last term in the equation above represents propagation effects which are discussed 
later in this section. Formulas for the source strength FV , the Strouhal number S , and 
the coefficients Zi pe through Z6 pe are provided below. 
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Coefficient Z4 Primary External Jet 

/ a r / v-vfi.80 


Z4 /;e = 0.6Log 10 (l +iV,){l - 0.2exp[-5(^ - 1 .s) 2 }- 0.6 - 


Z L — - 0.6 

c„ n 


W f 1 -8 <9 ) ( 1.8# 

-0.2^+0Aj3„ — - 0.6 - 0.6 + 0.1 ^-0.6 

Wp \ * J l * 

-0.2exp[-lo(# p -l) 2 _ 

+ C4l(P en ,l div ) + C42(E aspl , E aspl , R eJ ) + C43(d^,M' a ff p ) 
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Adjustment for lobe penetration and divergence angle 


C41 = 0.2P en - 0.01l div 


(67) 


Adjustment for aspect ratio and ejector length to inner diameter ratio 
C42 = -0.3 E„ 


J aspl 


+ 0.1 E 


asp 



r 1.80 p 

v 

1- 

P 

- 0.6 


Y JT 

)_ 


fi. 80 „ 

y 

1- 

p 

- 0.6 


{ * 

)_ 


Adjustment for azimuth angle 

C43 = 0.2 exp [5 - (0 p - 1 . 9) 2 - 2(1 - d: ) 
- 2 V^T exp(- 5(3 p ) 


1.80 „ 


pe 


It 


w. 


w„ 


( 68 ) 


(69) 


Coefficient Z5 

Primary External Jet 


Z5 pe = 50 + 20exp[- [d p - 2.6) 2 _ 

(70) 

Coefficient Z6 

Primary External Jet 



Z6 „„ = 137 - 40| 22 - 0.6 + 46exp[- (0 p - 2.5 ) 2 ]- 8Log 10 (l + N, ) - 5 Yf (71) 


+ 2 exp[- 10(0, - 2.5) 2 ] + (2 + 3/3 p ^ - 0.6 

- 3exp[- 10(0 p - 2.3 ) 2 ]+ 2exp[- io(0 p - 1 . 9) 2 

+ C61(l div ) + C62(M' a ffi p ) + C63(d ip , E aspt ) + C64(R lin ,R ej ) + C65(C fg ,M' a ) 


Adjustment for divergence angle 




f 1.80 ) 


C61 = 0.4/ d;v 

0.4- 

P 

-0.6 


l x J 



(72) 


Adjustment for flight effects 

C62 = 15 jMl exp(- 5/3 p )exp[- io(0 p - 2.4) 2 _ 


Adjustment for aspect ratio and azimuth angle 


(73) 
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Adjustment for thrust loss and flight effects 

C65 --40(1-Cjl-M:) (76) 

Mixed Jet Component of Jet Mixing Noise 

The low frequency noise, termed “mixed or merged noise”, is generated by the mixing of 
the so-called “mixed” jet (the initially rectangular jet having decayed to a near circular 
shape) with the freestream. The peak noise source is several equivalent diameters 
downstream of the ejector exit, where the length scale is quite large, usually on the order 
of the ejector equivalent diameter. The characteristics of the merged noise are similar to 
those of subsonic round jets. The sound pressure level for the mixed jet is 

SPL„ ) - [(Zl„ )Log m (FV„ ) + (Z2„ )ILo* i0 S. - (z 3. )Log „ (. FV „ ) - (Z4„, )] 2 

+ (Z5 . )LogjFV„)+ (Z6„ ) + ASM,,, 

The last term in the equation above represents propagation effects which are discussed 
later in this section. Formulas for the source strength FV m , the Strouhal number S m , and 
the coefficients Z\ m through Z6 m are provided below. 


Source Strength Function 


Mixed Jet Component 


FV„, = 


V -V V + V 

m a \ I m a 


( V 'V^ 2 0 2 ( S V 

n = \ 0.6 + exp - 0.3 8 m + — — - 2.7 

c „ 0.2 + 5,,, 1 + 5,„ 
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Coefficient Z4 


Mixed Jet Component 


f 4.50 ■ ) 0.2V - 0.7V A 

Z 4 m = 0.44 - 0.5 exp ^ - 4 + - 0.2 Log w 1 + — 

V 71 ) C a \ \ 

1 + Log 10 (1 + N i ) 

f V "V" 18^? 5 2 

+ 0.1 — ^ + 5/3 m — — - 0.6 + C4l(E ,C f ) 

c„ A n ) 


where 0 mj = 0 m 


A a a 8 

and - < 0 mj <-ji 


59 







Adjustment for thrust loss 


C4 1 = 0. 1 5E asp - (l - C js ) 


( 85 ) 


Coefficient Z5 

Mixed Jet Component 


8 10. ( 1.80 V 

Z5 = 34 + ^ - 20 a - 0.6 

TC X jt J 


(86) 



Coefficient Z6 

Mixed Jet Component 



Z6,„ =123 + 12 


1.80. 




- 0.6 


J 


y Jt 

+ C61(C f ,M') + C62(E.,d u ,,M') 


7exp[-8(«„ -2.4) 2 ]+5 J 8„exp[-15(9 m -2.3) 2 ] (87) 


Adjustment for thrust loss and flight effects 

C61 = -40(l-Cjl-Af:) (88) 


Adjustment for aspect ratio, azimuth angle, and flight effects 
C62 - 2£„,, l cxp[- Kl(0„ - 1 ,8) 2 ] (l + ) 


1.80, 


- 0.6 


(89) 


Broadband Shock Noise 

For a supersonic jet inside an ejector, formation of broadband shocks in the internal flow 
is unavoidable due to the presence of pressure gradients and non-ideal expansion at off- 
design operating conditions. For very high nozzle pressure ratios, the mixed jet velocity 
may be supersonic at the ejector exit, which could lead to shock formations in the 
primary external plume. For certain operating conditions, the contribution of shock noise 
to the total noise produced by the mixer-ejector could be significant. The shock noise 
prediction procedure incorporated into this module is based on the instability wave/shock 
cell interaction source model developed by Tam [10,11]. The prediction of broadband 
shock noise in the internal jet region follows Tam’s model for circular jets. Empirical 
adjustment is made for the rectangular geometry, liner effects, and frequency. Shocks in 
the primary external jet are predicted by Tam’s model for rectangular nozzles [12]. 
Empirical adjustments are made for flight effects. 
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Primary Internal Broadband Shock Noise 

The mean-square acoustic pressure for the broadband shock noise associated with the 
primary internal jet is 


l (r s ,0,<p,f) 


shock, pi 


ha 2 

r; 




(90) 


where II is the overall acoustic power, A 2 is the area of a single lobe given by, 


\ = A p / N, 


(91) 


and F is the spectral distribution function. The right side of equation 90 represents the 
mean-square pressure for a stationary observer in nozzle fixed coordinates. Figure 3.9.4 
shows the relationship between the nozzle fixed coordinate system and the ground based 
observer coordinate system. Since the noise at the observer is independent of the 
coordinate system, the mean-square pressure can be transformed from nozzle fixed 
coordinates to the ground fixed coordinates by the following geometric relations 


ijj = tan 


sin# 

cos# - M a 


(92) 


and 


R 


sin# 
sin ij) 


(93) 


Primary Internal Shock Acoustic Power 

The broadband shock acoustic power adjusted for the mixer configuration is 


n = 


_ v 

cl} A 2 — | 
c; 


r 


where 


and 


c =0.0001 


i 2 =l-M fl 2 sin 2 M 


(94) 

(95) 

(96) 
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L is the shock cell spacing parameter and is assigned a constant value 


L = 3.3 


( 97 ) 


The shock cell strength function ^ is defined for over and underexpanded jets as 


A =< 


X 


1 + 3 ^ A, 


J>P 


X 


l + 6x 


M >M d Underexpanded Jet 
M < M d Overexpanded Jet 


(98) 


where 


X- 


M) 

1 -M 2 

,p 1Y1 d 

\ 

1 + 1 

W-K 


(99) 


A 2 is the area of a single lobe and A. is the fully expanded jet area given by 


A i,p - \ 



l( Yp+ 1 N 

M “ ) 



M ,r y 


( 100 ) 


Spectral Distribution Function 

The one-third-octave band spectrum function is obtained by integrating 



(101) 

The upper and lower limits of integration are 


i 

(102) 

and 

l 


f, = 2~ 6 f 

(103) 


where / is the one-third-octave band center frequency. The argument of the integral is 
evaluated at 20 points within a one-third-octave band and integrated numerically using 
Simpson's rule. In equation 101, f re is 
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(104) 


f = 

J re 


AD 


2 


D -D . 

ej plug 


f 


The parameter f H represents the sum of the modal contributions for a given frequency 


35 

-i 


3 o. 


W(> 


a. 


exp 



( f V 

J m j 

1 4- M 

( MJ + c os(^)Y 

2 (V 

2 2 ] 

L 2 


yf re J 

1 + m c 

l ^ -K) l 

UcJ 

2 In 2 


(105) 


Summation takes place over the first 35 modes or until the quantity (f re D, / V ) becomes 
greater than 5.5, at which point the summation is terminated. The factors o m are the 
zeros of the zero-order Bessel function and /, (cr m ) is the first order Bessel function. 
Values for o m and /, (o m ) are given in table 3.9.1. The modal frequency parameter f m 
is defined in terms of the wave number k m and the convection velocity V c by 


k V \ 
m C \\ + M r 

M J + cos{i/j)]} 

— 1 r 

2jv 

(1 ~M 2 a )z JJ 


The convection velocity is 
V c -0.1V p+ 0.3V a 

and 

M c =V c /c a 

The wave number is given by 

, 2 2 c “ g “ , 1 . 2....35 

The parameters c m are given by 
q = 1.1596 - 0.1113M p 

c 2 =1.1 

and wave numbers 3 through 35 are given by 


(106) 

(107) 

(108) 

(109) 

( 110 ) 
( 111 ) 

( 112 ) 
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Primary Internal Broadband SPL 

The broadband shock noise associated with the primary internal jet is determined by 

SPL,^, - 10 Log m ((p 2 )/ pl f )+ R„ (l + SM J - 3^ + 0.5 exp[- 5(s„ - 1 4 } ( 1 13) 

x 2 - ex pip... - P., )]cxp[- 0.s(Log,„S p - Log m S m } J+ ASPL M 

Lu and Viswanathan developed empirical adjustment for geometry, liner effects, and 
frequency. The last term in the equation above represents propagation effects, which are 
discussed later in this section. 
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Table 3.9.1 - Values of o m and 7, (or,, ) 


Mode number, m 

°m 


1 

2.404826 

0.519147 

2 

5.520078 

- 0.340265 

3 

8.653728 

0.271452 

4 

11.791534 

- 0.232460 

5 

14.930918 

0.206546 

6 

18.071064 

-0.187729 

7 

21.211637 

0.173266 

8 

24.352472 

-0.161702 

9 

27.493479 

0.152181 

10 

30.634606 

-0.144166 

11 

33.775820 

0.137297 

12 

36.917100 

-0.131325 

13 

40.058426 

0.126069 

14 

43.199792 

-0.121399 

15 

46.341188 

0.117211 

16 

49.482613 

-0.113429 

17 

52.624054 

0.109991 

18 

55.765514 

-0.106848 

19 

58.906986 

0.103960 

20 

62.048470 

-0.101293 

21 

65.189964 

0.098823 

22 

68.331474 

- 0.096524 

23 

71.472984 

0.094379 

24 

74.614510 

-0.092371 

25 

77.756027 

0.090485 

26 

80.897560 

-0.088711 

27 

84.039093 

0.087037 

28 

87.180634 

- 0.085454 

29 

90.322174 

0.083955 

30 

93.463722 

- 0.082532 

31 

96.605270 

0.081179 

32 

99.746819 

- 0.079890 

33 

102.888374 

0.078661 

34 

106.029938 

- 0.077487 

35 

109.171494 

0.076364 
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Mixed Jet Broadband Shock Noise 

The mean-square acoustic pressure for the broadband shock noise associated with the 
mixed jet is 


(p 2 (r, 6, (/>,/)) 


shock , m 


R 2 


f(v, </>,/) 


(114) 


where II is the acoustic power and A. m is the fully expanded jet mach number. The 
acoustic power of a jet with a rectangular nozzle is 


IT = 


cA;M\ m 


r 


\ 




where 


c = 2.886 x 10' 4 

and 


(115) 


(116) 

(117) 


The shock cell strength function for over and under expanded jet is 


A?=j 


A 2 M j m > M d (Under expanded jet) 


A 


2 J,m 


Mj,m ^M d {0 ver expanded jet) 


(118) 


— 2 

where A is 


A 2 = 


1 + 


V 

2 

E asp 2 

1 

' V 

M l m ~ l 

2 

1 

2 J 

1 + Ej 

i + \ir m -i) 

v 2 J 

i +E as ; 

( 

\ 2 

E a sp 

3 + 

2 

1 + 

( 

M lm - 1 

\ 2 

1 

V 2 

1 + Ej 

i + ^Ou, -i) 

V 2 

1 + V 

/ 


(119) 


The fully expanded jet area is calculated using the following equation 
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f 


(120) 






1 ( y m + 1 


M d 

v M jJ 


where M d is the design jet Mach number and M . m is the fully expanded jet Mach 
number. The rectangular nozzle exit area is 


A p =b ejKj 


( 121 ) 


The parameter represents the aspect ratio of the rectangular nozzle defined by 
b„ 


E ~- SL 

a sp h . 

ej 


( 122 ) 


Rectangular Nozzle Spectrum Function 

The one-third-octave band spectrum function in the fly-over plane (i.e. 0= 0 ) is obtained 
by integrating 


f(^ = o,/)= \ F ^ F ^df 

J re 


(123) 


The argument of the integral is evaluated at 20 points within a one-third-octave band and 
integrated numerically using Simpson's rule. 

The functions F l and F 2 are expressions for the odd and even modes, respectively. 

The function F 1 is 


30 30 


exp 


p.-ZE- 



K , M c (cosy + gM a ) ' 

Y v o 

2 7 1 

L 2 




2 In 2 


m = 1 n = 1 


(2m - 1) 2 (2/t - ljl + 2 {n - l)]exp{- (E asp - l)} 


(124) 


The convection velocity, V c , of the large turbulent structure of the jet is given by 
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v,- 0.5-0.2exp 4 ^— 1 -0.06[^--l] (V.-Vj+V, 

'l ^ / * r. 


(125) 


and the half- width of the similarity noise source is 


/ r p ' 

L = 3.0^- 1.0+ 1.114-0.36-^- M, 


(126) 


with a core length to jet diameter ratio defined as 


4-3 + 1.2 M 2 . T tJ JT a >1.0 


J ' rj-t \ 

%c = |4.3 + 1.2M7,„+1.2 1.0--^ T.JT <1.0 


(127) 


The spacing of the odd shock cells is defined as 


L = 


2-\j M j, m ~ 1 L0+ 0.812-0.254^ M c 


\ 2 f \ 2 

2 n - 1 2m - 1 
b- h : 


m = 1,2,3, ...30 
n = 1,2,3,. ..30 


(128) 


where C m is an empirical correction factor, which refines the spacing for the /nth mode 


c = k + [0.756-0.195 M jm - (a - l)]exp{- \()A(E asp - 1) 2 } m = 1 
| a, + (l . 1 - a, )exp{- 0.4(e - 1) 2 } m > 2 


(129) 


and where 


a, = 1.0 + 0.578 exp{- [4.75 (m Jm - 1 


(130) 


The fully expanded jet height, h f , is defined with respect to the fully expanded jet area as 


/ / A i- m 

h , = Kj -T— 

b ej 


(131) 
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The function F, is defined as 


exp- 


j 

-I- 



K x M c {cosip +£M a ) 


2 ^ i 

Lr 


_L m (^-M a cos^f 4(l 

Uj 

2 In 2 


4m 2 


( 132 ) 


with the spacing of the even mode shock cells given by 


2 -1 


= 


1.0 + 


0.812-0.254 


1 t,m 


M, 


l a J 


c „ 


m = 1,2,3 


f . 3 


\ h ej J 


r 2 n? 

V h J J 


(133) 


Mixed Jet Broadband SPL 

The broadband shock noise associated with the mixed jet is determined by 
SPL^ l .M.e,'p,r,).10Lo gm ((p , -)lp^)-HE mr J) 


- 10 1 exp 




k TC 


J 


- 1 + 3 M' ^ + A SPL 


'shock 


(134) 


The prediction is extended to azimuthal angles other than zero, ^^ asp ,(p) which was 
developed by Ponton, Manning and Seiner [13]. The azimuthal correction factor is given 

by 




0.3992 -0.5228£ fl?p +0.1 192£ fl , p j 


\jzI2 j 


+ 

- 1.48 + 1.616AE - 0.192 E n J 

f $ \ 


cisp asp 

/ 2y 


(135) 


E 

and depends on the aspect ratio of the jet as P and the azimuthal angle <!> in radians. The 
last term in the equation for SPL is an empirical flight effects adjustment developed by 
Lu and Viswanathan. The last term in the SPL equation above represents propagation 
effects, which are discussed next. 
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Propagation Adjustments 

For ease of discussion, the term propagation effects will also include flight effects and the 
number of engines. Propagation effects must be incorporated into the noise predictions 
as a function of frequency, emission angle, and radiation distance. The adjustments 
required to predict the noise at the observer are different for each of the five jet noise 
components. The jet mixing sources for example include geometric near-field and 
acoustic near- field effects, which are not a part of the shock noise model. This section 
describes the propagation effects for the five jet noise components. 

Geometric Near-Field 

The jet mixing noise components include adjustments for acoustic near-field and 
geometric near-field effects. When the observer (microphone) is close to the jet, the 
acoustic field is simulated better by a line source than a point source. The transition from 
the near-field to the far-field is made by the following empirical formula for the mean- 
square acoustic pressure 

(p 2 ) cc (D/r) 2 {l + b/r)~ l (136) 

The parameter D is the jet diameter for each component of the jet, r is the distance from 
the apparent source to the observer, and b is an empirically defined source length. 
Converting the acoustic pressure to SPL, one finds that 20 Log l0 (D/r) is the far-field 
spherical divergence and -10Log 10 (l + b/r) is the near-field correction which decreases 
the extrapolated far-field SPL in the geometric near-field. The source length b is 
empirically determined and is listed in tables 3.9.2 and 3.9.3. The empirical geometric 
near-field effect is calculated for each jet component (primary internal, primary external, 
mixed) in the prediction model. 

Acoustic Near -Field 


The acoustic near-field is also determined empirically. In this case, the wave length is used to 
measure the acoustic near-field distance of the observer from the source. The mean-square 
acoustic pressure is approximated by 


(. P 2 ) a {A/rf l + l.3{A/r) 2 


(137) 


where A =c a /f is the wave length in air at rest. When the mean-square acoustic pressure is 
converted to SPL, the term 20 Log w (A/r) is the far-field spherical divergence. The near-field 


correction is 10 Log l0 


1 + 1-3 (y/rY 


70 



Propagation Effects for the Primary Internal and External Jets 

The following equation defines the propagation effects for the primary internal and 
external jets. 


A SPL. = 20 Log 


10 


f D ^ 


( b^ 



fc 1 

2 

r J 

- 10 Log l0 

1 + — 

r l 

+ 10 Log l0 

1 + 1.3 

{ r pf) 


V p ) 


^ p ) 





( 138 ) 


+ 20 Log 


f P A 


10 


\ P std J 


+ 20 Log 


r P pi + Pa^ 


10 


2 P. 


10Log, 0 (lV ) 


-| A„ t ,(f,rh)\r r +G-EGA 


Table 3.9.2 defines the terms in the equation 138. 

Table 3.9.2 - Definition of A SPL Adjustments for 
Primary Internal and External Jet Noise Components 


Adjustment for: 


Term 


Spherical spreading 

20 Log l0 

( ° P " 

r 

V P J 


Geometric near-field 

-10 Log l0 
where b = 2D 

u±-j 

r 

\ p J 

JL 

C a 

X 1/2 
/ 

' I / . 

Acoustic near-field 

10 Log m 

1 + 1.3 Cfl 

r f 

V P J J 

2 


Ambient pressure 

20 Log m 

f P 

a 

p 

K r std y 

' 

Density 

onr (Ppi + Pa 

20 Log io 

l 2 Pa , 


Number of sources 

10 Log m [N e J 

Atmospheric absorption 

-\ A absU pi ,rh)\r p 

f = f 

p 1-M a cos0 p 

Ground reflections 

G 

Extra Ground Attenuation 

EGA 
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Propagation Effects for the Mixed Jet 


The following equation defines the propagation effects for the mixed jet component 


A SPL m = 20 Log l0 


f D ^ 


r ^ 



~ 10 Log w 

i+— 

+ 10 Log m 

v r m ) 


l r pj 



1 + 1.3 


f \ 2 
c„ 


+20 Log, 


10 
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+20 Log 

V P. 1 
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P,n + Pa 


10 


^ 2 * ) 


K r rfj 
10 Log 10 (N eng ) 


(139) 


- \A a bs (f ’ rh)\ r s +G - EGA 


Table 3.9.3 defines the terms in the equation above. 

Table 3.9.3 - Definition of A SPL Adjustments 
for Mixed Jet Noise Components 


Adjustment for: 

Term 

Spherical spreading 

20 Log 10 

f D \ 

l On J 


Geometric near-field 

- 10 Log m 
where b = 2 D m 

r 

1 

V 

+ 

b 

+ — 

V 

f Dj 

{ f 

\ 

\ 

''a 

) 

1/2 

Acoustic near-field 

10 Log w 

1 + 1.3 

r 

c a 

r f 

v p j 

y“ 

) 


Ambient pressure 

20 Log l0 

f P \ 

P 

\ r std ) 


Density 

f 

20 Log l0 

V 2 P a 

\ 

) 

Number of sources 

10 Log m {N e J 

Atmospheric absorption 

-\ A ab,(fm’ rh )\ r m 

f = f 

Jm l-M a co&0 m 

Ground reflections 

G 

Extra Ground Attenuation 

EGA 
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Propagation Effects for the Broadband Shock Components 

The following equation defines the propagation effects for the mixed jet component 


ASPL, r , oct .l0Log„(N, nl )+20Log, a 


T(-.) 


-| A^f,rh)\r,+G-EGA (140) 


Table 3.9.4 defines the terms in equation 140. 

Table 3.9.4 - Definition of A SPL Adjustments 


Adjustment for: 

Term 

Number of sources 

lOios.olA'.J 

Ambient Pressure 

OQ L02 f ^ ^ ^ 

ZAJL/Og IQ / \ 2 / \ 

Atmospheric absorption 

-\ A abs(f,rh) \r s 

Ground reflections 

G 

Extra Ground Attenuation 

EGA 
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Total Jet Noise 


The total jet noise is the logarithmic sum of the five jet noise components, 

SPL total (/ , 9, (j), r s ) = 10Log 1() 2 (141) 

where 


2 __ 10 ^V 10 _|_ 1 Q ^ jq 5PL w /10 ^ J 0 ^^ shock,pi / 10 10 5PW ’ m/10 


(142) 





1/3 Octave Band SPL 


Component noise spectral shape 



0.5 1 1.5 2 

Log (Strouhal Number) 


Figure 3.9.2 Effects of parametric variation on component spectra 
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Reference 

Center 

Station 


(a) Model in-flow coordinates 



(b) Full-scale flyover coordinates 


Figure 3.9.3 Schematic of jet noise propagation from a distributed source location 
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Source position at 
Source position at time of reception 

time of emission M a 



Observer 


Figure 3.9.4 Schematic of emission and reception coordinates 
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3.10 Jet Noise (Stone Methodology) 


The method described in this section was developed during the HSR program by Modem 
Technologies Corporation (MTC) for General Electric Aircraft Engines [14]. The 
method predicts jet mixing and broadband shock noise for a 2-D mixer ejector nozzle 
system. For purposes of analysis, the jet is conceptually divided into three regions: the 
internal region within the ejector, the premerged region immediately down stream of the 
ejector, and the merged region far enough down stream that primary and secondary 
regions of the flow can no longer be identified. Five noise sources may result from these 
regions: internal mixing noise, internal broadband shock noise, premerged mixing noise, 
external broadband shock noise, and merged mixing noise. The method consists of semi- 
empirical correlations of component noise characteristics in terms of geometric and 
aero/thermodynamic parameters. Each source is modeled by analogy to a round 
convergent nozzle, with different characteristic velocities and lengths along with 
appropriate physical properties. The model was developed from experimental data taken 
in GE's Cell 41 anechoic free -jet facility. The code was verified by comparisons with 
GEAE databases over a range of variables including simulated flight Mach numbers from 
0 to 0.32 and primary jet velocities from 1 100 to 2800 ft/sec. The source code provided 
by MTC has been incorporated into the HSRNOISE code with minor modifications to 
accommodate the HSRNOISE code design and coordinate system. 


Symbols 

A, 

\n 

Aec 

Ah 


^ amb 

CER 

r 

source 

D e q,ej 
D 'h.pri 

DM 

f 

fpk 

F 

S 

h„; 


'bulk, end 


“‘bulk, start 


perf , end 


2 

ejector exit area, ft - 

primary nozzle flow area, ft - 
secondary nozzle flow area, ft - 

primary nozzle throat area, ft - ( input parameter- APT) 

sound speed, ft/sec 

ambient sound speed, ft/sec 

core expansion ratio, ( input parameter - CER) 

coefficient of component noise sources 

ejector nozzle equivalent diameter, ft {^4A ej /jr) 

primary nozzle hydraulic diameter, ft ( input parameter - DHP) 
degree of mixing, ( input parameter - DM) 
one-third-octave band center frequency 

peak frequency of broadband suppression, ( input parameter - FPK) 
directivity function, dB 

2 

gravitational constant, 32.174 ft/sec 
ejector height, ft ( input parameter - HEX) 
ejector length, ft ( input parameter — LE) 

length from primary exit to end of bulk treatment, ft ( input parameter - LBE) 
length from primary exit to start of bulk treatment, ft ( input parameter - LBS) 
length from primary exit to end of perforate treatment, ft ( input parameter - LPE) 
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perf, start 

M a 

M J 

M d 

M 

M pri 

p 

amb 


t,pri 

PEN 

r s 

R 

S 

St 

St pk 
SAR 
T 

amb 


t,pri 


V 

V a 

V, 


w 

« 7 

a, 

A dB, 


bulk,pk 


A dB 


perf.pk 


pri 


P 

6 

co 


Subscripts 

amb 

bulk 


length from primary exit to start of perforate treatment, ft ( input parameter - LPS) 
aircraft Mach number 
jet Mach number 

primary nozzle design Mach number 
molecular weight 

molecular weight of primary jet, ( input parameter - MMC) 

ambient pressure, lb/ff ( input parameter - PAMB) 

2 

static pressure, lb/ff 
total pressure, lb/ff 

total pressure of the primary jet upstream of the nozzle, lb/ff ( input parameter - PC) 

mixer spoke penetration to nozzle full scale height ratio ( input parameter - PEN) 

distance from source to observer, ft 

gas constant, 53.35 ft-lbf / lbm °R 

spectrum directivity function 

Strouhal number 

peak Strouhal number 

suppressor area ratio ( input parameter - SAR) 
ambient temperature, °R 
static temperature, °R 
total temperature, °R 

total temperature of the mixed jets, °R ( input parameter - TEX) 

total temperature of the primary jet upstream of the nozzle, °R ( input parameter - TC) 

velocity, ft/sec 

flight velocity, ft/sec ( input parameter - VI) 

exit velocity of mixed jets, ft/sec ( input parameter - VEX) 

pumping ratio ^ = ™ P n ) ( input parameter - WSWP) 

ejector width, ft ( input parameter - WEX) 

angle of misalignment between primary nozzle and ejector flow, degrees 
jet angle-of-attack, degrees 

peak bulk suppression at FPK per foot of treatment ( input parameter - SUPPK) 

Peak perforate suppression, dB/wavelength ( input parameter - DL) 
ratio of specific heat of ambient flow 

primary jet ratio of specific heat ( input parameter - GAMMAC) 
azimuth angle, degrees 
density, lbm/ft 3 

polar directivity angle, degrees 
density exponent 


ambient 

bulk 
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D 

Doppler 

e 

effective 

end 

end 

eq 

equivalent 

exit 

exit 

fe 

fully expanded 

guess 

guess 

IM 

internal mixing noise 

IS 

internal shock noise 

mix 

mixed 

MM 

merged mixing noise 

perf 

perforate 

pk 

peak 

pri 

primary 

PM 

premerged mixed noise 

PS 

premerged shock noise 

s 

static 

sec 

secondary 

start 

start 

sup 

suppression 

t 

total 

th 

throat 

Superscripts 

* 

nondimensional 
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Component Formulation 

For purposes of this analysis, the jet is conceptually divided into three regions as depicted in 
figure 3.10.1. The regions are; the internal region within the ejector, the premerged region 
immediately downstream of the ejector, and the merged region far enough downstream that the 
primary and secondary regions of flow can no longer be identified. Five noise sources may 
result from these three regions; internal mixing noise, internal broadband shock noise, premerged 
mixing noise, premerged broadband shock noise, and merged mixing noise. The formulations 
described in this section represent free-field lossless sound pressure levels. Since propagation 
and flight effects for each source are different, they are also defined in this section. 


Primary Flow 


Secondary Flow 




Premerged Merged 

Region Region 



Figure 3.10.1 Two-dimensional mixer ejector noise generating regions 
Iteration Procedure 


Prior to computing the component acoustic pressures, thermodynamic properties of the primary 
and secondary flow regions are computed using an iterative procedure. The procedure is based 
on conservation of mass, momentum, and energy and begins with an initial guess of the total-to- 
static pressure ratio of the secondary flow. From this the mass flow rates of the primary and 
secondary flows are computed. The ratio of the computed mass flow rates is compared to the 
input pumping ratio. If the initial pass does not converge, the pressure ratio guess is revised, and 
the procedure is repeated, until the procedure converges. In some cases the procedure may not 
converge. Should this occur, it is an indication that a consistent set of parameters cannot be 
computed from the input parameters provided. 


Step 1. Assume an initial guess for the total-to-static pressure ratio of the secondary 
stream. 


guess, sec 


t, sec 


amb 


= 1.108 


( 1 ) 


Assume an initial guess for the primary flow area is equal to the primary nozzle exit area. 


A = A 

pri,exit pri 


( 2 ) 


The primary nozzle area is obtained from the primary nozzle throat area A th and the core 
expansion ratio CER . 
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( 3 ) 

( 4 ) 

( 5 ) 

(6) 

( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 
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Step 3. Compute the primary flow parameters. 


Primary flow static temperature: 


T =T 

s,pri t,pri 


r y-'/JV; 

1 


p 

V >7 


Total to static pressure ratio: 


P = P .IP 

pri t,pri s, sec 


( 13 ) 


(14) 


Ratio of specific heat: 


y pri = 1.42427 - 1 .07556x10 “ 4 AT + 2.67072x10 ~ 8 AT 2 
- 2.96091 lxl0~ 12 A T 3 - 3.82470xl0“ 17 AT 4 


(15) 


A T = T Lpri -T amh (16) 

Sound speed: 

C pri = pri S R P Js. pri ( i7 ) 


R pri = 


1544 


M pri = 28.9691 + 3.3221xl0' 12 AT + 2.77821x10“'° AT 2 
Primary jet Mach number: 


(18) 

(19) 


M 


J’Pn 



( 20 ) 


Primary flow velocity: 


V . = M ■ c 

pri j, pri pri 


Primary flow density: 


P pri 


R T 

pri ± s,pri 


( 21 ) 


( 22 ) 
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Primary flow exit area: 

1 (y P n+ 1) 

2 {jvhI 

(23) 

Primary weight flow rate: 

w . = p A V . (24) 

rK /?n r pri pri v pri V / 

Step 4. Compute secondary mass flow rate from pumping ratio supplied as input. 

w ^(calculated) = w* w pri (25) 

Step 5. Compare calculated secondary weight flow rate to secondary mass flow rate 
supplied as input. If |r -l|> 0.0001 where 



vi’ sc c (calculated ) 
w scc (input) 


(26) 


Set new guess for pressure ratio to 

P;,,Jnew).0.gp;,,Jold) + 02(p;,,Jold)f (27) 

and repeat steps 2 through 5. If |F — 1| < 0.0001 the procedure converged. 

Mixed Jet Parameters 

Two options are provided to determine the properties of the mixed jet. Option 1 allows 
the user to input the total temperature and jet velocity of the mixed jet from which other 
properties are computed. Option 2 computes the mixed jet properties on a mass average 
basis of the primary and secondary streams. 

Option 1 : When T t imx and V mix are provided as input 




(28) 


where the ratio of specific heat of the mixed flow y m is the mass average of the ratio of 
specific heats from the primary and secondary streams. 


Ypri + ™ Y 

1 + w 


( 29 ) 
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R g 44 

M .... 


rr M pri + ™ M «> c 

M mix = ~ ^ 77 

1 + w 


( 30 ) 


(31) 


Option 2: When T t mix and V mix are not provided as input, the static temperature is 
computed as a mass average of the two streams. The static temperature at the ejector exit 
is given by 


f T . 

+ w T 

\ 

f p 

X 

s,pri 

s 

sec 

amb 


^ 1 + w* 

7 

p 

V sec 

7 


l-l/JV. 


(32) 


The mixed jet velocity and total temperature are computed internally by the following 
equations. 


w . + w 

y pri sec 


P mix / 


(33) 


T t .mix ^ s .mix 


1- 


1 


Y mix ) 


V. 

mix 


2 gR n 


(34) 


The mixed jet density, sound speed, and Mach number are computed by 


P mix 


amb 


R T 

' mix 1 s.mix 


c = I gy A 1 T 

mix ■y o / mix mix s.mix 


(35) 

(36) 


V . 

M = m,x 

j,mix 


(37) 


Mixer Ejector Correlation Parameters 

In addition to the thermodynamic properties in the three noise generating regions, 
parameters such as the degree of mixing and the suppressor area ratio were established to 
correlate the noise characteristics of the mixer ejector nozzle system. 
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Degree of Mixing 

The degree of mixing parameter enables noise characteristics to be empirically correlated 
to measured flow field characteristics. The characteristic velocity of the premerged 
region is bounded by the fully mixed jet velocity at the low end and by the fully mixed 
primary jet at the high end. The degree of mixing parameter is defined as the ratio of the 
difference between the fully expanded primary jet velocity and the characteristic 
secondary velocity of the premerged region, to the difference between the primary 
velocity and the fully mixed jet velocity. An empirical relation for DM was obtained as a 
function of the suppressor area ratio, the penetration, and the ratio of the ejector length to 
the hydraulic diameter of the primary nozzle and the primary nozzle pressure ratio. 
Algebraically, the degree of mixing is given by 


DM 


0.93 Z 
2.3 + Z 


(38) 


where 

, SAR(PENf{L,/D t J 

(. npr) l 5 


(39) 


Suppressor Area Ratio 

The suppressor area ratio SAR is defined as the ratio of the total mixing area (i.e. sum of 
the primary and secondary flow areas) to the throat area of the primary jet: 

SAR = Ap " + Asec (40) 

A,. 


Mixer Area Ratio 


The mixer area ratio and the geometric factor are also used to correlate the noise of the 
mixer ejector. The mixer area ratio is defined as the ratio of the mixer ejector exit area to 
the total mixing area: 


MAR 


\ri + Aec 


(41) 


Geometric Factor 


The geometric factor is used to define the azimuthal directivity effects of the mixer 
ejector nozzle. The geometric factor GF is a function of the aspect ratio and the azimuth 
angle and is defined as 
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GF = % kei ! Wei ^ C ° S ^ ^ “ ^ cor 

1 sin^ (p < </> cor 

where the comer angle is given by 

Kr = tan -'(^./wj 

Figure 3.10.2 defines the mixer ejector coordinate system. 

w 0 




( 42 ) 


(43) 


ej 


► I 


<\>=0 



h e j 


Figure 3.10.2 Two-dimensional mixer ejector nozzle coordinate system 
Internal Mixing Noise 

The free-field lossless sound pressure level for the primary internal mixing noise 
component is given by 

SPL im = 147 + 10 Log w (A pi ,/r ; )+ 10 'a IM Log 10 (p pri / p sec ) 

+ 60l °£io ( y e , pri /c a )- A w - 19 F us 0 9)Log w (GF) 


- F 2 ,s(0') + 60 Lo g w 

- S IM (f,6)-AdB sup 


1 + M 


e 

V 150y 


- 40 Log w {\- M a cos 6 D ) 


(44) 


The amplitude of the internal mixing noise is determined by the primary flow area, the 
density and jet velocity of the primary and secondary streams, and the length and 
equivalent diameter of the ejector. By analogy to circular jets, temperature effects are 
determined by the density ratio and the density ratio exponent co IM . The density ratio 
exponent is given by 


COjm 


0.6 + (V e ,, M /c 0 ) 3 ' 5 


(45) 
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where V e 1M is the internal mixing noise effective velocity defined as 


V,.m -V,Jl-('WO cosa i] W (46) 

The angle a I represents the angle of misalignment between the primary nozzle and the 
ejector flow. 

The hardwall shielding term A flw is defined as the ratio of the ejector length to the 
equivalent exit diameter of the ejector. The hardwall shielding term is given by 

A™ -0.5 LJD„„ (47) 

where the equivalent exit diameter of the ejector is 

J ~4A~ 

(48) 

and the exit area of the mixer ejector is 

A ej = W ejKj (49) 

Internal Mixing Noise Directivity 

The internal mixing noise directivity in the absence of flight effects is defined by two 
terms. The term 19 F, lM (O)Log U) (GF ) (in equation 44) combines azimuthal and polar 

angle effects. The function fj IM ( 6 ) is given by 


Fum (0) = 


o: 


,M e\ M ^ 125 


, _ (ftjM — [ 25 ) 25 

90 lM 


where the effective directivity angle 6' m is defined as 

y V ' 1 
e' 


V C a j 


The function F 2 IM {(f ) is obtained from table 3.10.1. 


(50) 


( 51 ) 



Table 3.10.1 Internal Mixing Noise Directivity Function 


e\ M 

F 

1 2,/M 

e\u 

F 

1 2,1 M 

0 

19.4 

110 

6.1 

10 

18.6 

120 

5.4 

20 

17.8 

130 

5.6 

30 

17.0 

140 

10.2 

40 

16.2 

150 

14.3 

50 

15.4 

160 

19.0 

60 

14.6 

170 

23.2 

70 

12.9 

180 

27.7 

80 

11.1 

190 

31.0 

90 

9.4 

200 

34.4 

100 

7.9 

250 

50.6 


Internal Mixing Noise Flight Effects 


Adjustments to the amplitude of the internal mixing noise level due to flight effects are 
determined by two terms. The term 60 Log w [l + M fl (^/150)'° is very small except for 


angles greater than 120 degrees. The convective amplification term 

-40Log 10 [l - M a cos 6^] is also used to describe the change in noise levels due to the 

motion of the aircraft. 


Internal Mixing Noise Spectrum Directivity Relations 

The internal mixing noise spectrum is correlated as a function of a modified Strouhal 
number, which includes azimuthal and aspect ratio effects as well as temperature effects. 
The length scale for the Strouhal number is the hydraulic diameter of the primary jet. 

The characteristic velocity is the internal mixing noise effective velocity. The internal 
mixing noise Strouhal number is defined as 


St IM = 


fD 


h,pri 


V 


4gf 


f y0.4(i+cos e\ u 

* t,pri 


e,IM 


T 

V s,pri J 


(l “ M a cosd D ) 


(52) 


and includes a Doppler shift. The internal mixing noise peak Strouhal number varies 
only slightly with the effective directivity angle as indicated by the values presented in 
table 3.10.2. The internal mixing noise spectrum function is obtained from table 3.10.3 
as a function of the effective directivity angle and of Log w (St' IM ) given by 

Log 10 {St IM ) = Log 10 ( St IM ) - Log 10 (St pk IM ) (53) 
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Table 3.10.2 Internal Mixing Noise Peak Strouhal Number 


e\u 

L °8m( St P kj M ) 

e\ M 

L °8m( St pkj M ) 

0 

0.6 

110 


10 

0.6 

120 



0.6 

130 



0.6 

140 



0.6 

150 



0.6 

160 



0.6 

170 



0.6 

180 


80 

0.5 

190 



0.5 



100 

0.5 

250 



90 






Table 3.10.3 Internal Mixing Noise Spectrum Directivity Function, S IM (f, 6) 


Log l0 (St, M ) 

0.0 

70.0 

80.0 

90.0 

Effective Directivity Angle, d' 1M , degrees 
100.0 110.0 120.0 130.0 140.0 150.0 160.0 

250.0 

- 3.6 

99.0 

99.0 

91.9 

97.1 

95.7 

77.7 

66.0 

57.8 

57.1 

56.8 

54.9 

54.9 

- 2.3 

43.7 

43.7 

45.1 

47.7 

47.0 

38.7 

34.2 

30.2 

32.4 

32.7 

31.5 

31.5 

- 2.2 

39.6 

39.6 

41.5 

43.9 

43.3 

35.7 

31.8 

28.2 

30.5 

30.8 

29.7 

29.7 

- 2.1 

36.2 

36.2 

38.1 

40.3 

39.7 

32.8 

29.4 

26.3 

28.6 

29.0 

28.0 

28.0 

- 2.0 

32.9 

32.9 

34.8 

36.8 

36.2 

30.0 

27.0 

24.4 

26.8 

27.3 

26.3 

26.3 

- 1.9 

29.7 

29.7 

31.6 

33.5 

32.9 

27.3 

24.8 

22.5 

25.1 

25.6 

24.6 

24.6 

- 1.8 

26.7 

26.7 

28.6 

30.4 

29.8 

24.7 

22.6 

20.7 

23.4 

23.9 

23.0 

23.0 

- 1.7 

23.9 

23.9 

25.8 

27.4 

26.8 

22.3 

20.5 

19.0 

21.7 

22.3 

21.5 

21.5 

- 1.6 

21.2 

21.2 

23.1 

24.6 

23.9 

20.0 

18.5 

17.3 

20.1 

20.8 

20.0 

20.0 

- 1.5 

18.6 

18.6 

20.5 

21.9 

21.2 

17.8 

16.6 

16.6 

18.6 

19.2 

18.5 

18.5 

- 1.4 

16.2 

16.2 

18.1 

19.4 

18.7 

15.7 

14.7 

14.0 

17.0 

17.7 

17.0 

17.0 

- 1.3 

14.0 

14.0 

15.9 

16.9 

16.3 

13.7 

13.0 

12.5 

15.5 

16.1 

15.5 

15.5 

- 1.2 

11.9 

11.9 

13.8 

14.7 

14.0 

11.8 

11.3 

10.9 

13.9 

14.6 

14.0 

14.0 

- 1.1 

9.9 

9.9 

11.8 

12.6 

11.9 

10.0 

9.6 

9.4 

12.4 

13.0 

12.5 

12.5 

- 1.0 

8.1 

8.1 

9.9 

10.6 

10.1 

8.5 

8.2 

8.1 

11.0 

11.6 

11.1 

11.1 

- 0.9 

6.6 

6.6 

8.1 

8.8 

8.4 

7.1 

7.0 

6.9 

9.7 

10.3 

9.9 

9.9 

- 0.8 

5.6 

5.6 

6.5 

7.2 

7.0 

6.0 

6.0 

6.0 

8.5 

9.2 

8.8 

8.8 

- 0.7 

4.8 

4.8 

5.3 

5.9 

5.9 

5.2 

5.2 

5.2 

7.5 

8.1 

7.8 

7.8 

- 0.6 

4.0 

4.0 

4.4 

4.9 

4.9 

4.5 

4.5 

4.6 

6.8 

7.1 

6.8 

6.8 

- 0.5 

3.4 

3.4 

3.8 

4.1 

4.1 

3.9 

4.0 

4.1 

6.2 

6.0 

5.8 

5.8 

- 0.4 

3.0 

3.0 

3.2 

3.4 

3.4 

3.3 

3.4 

3.6 

5.3 

5.0 

4.8 

4.8 

- 0.3 

2.4 

2.4 

2.5 

2.7 

2.7 

2.7 

2.8 

2.8 

4.0 

3.7 

3.5 

3.5 

- 0.2 

1.8 

1.8 

1.8 

1.9 

2.0 

2.0 

2.0 

2.0 

2.7 

2.6 

2.5 

2.5 

- 0.1 

0.9 

0.9 

0.9 

1.0 

1.0 

1.0 

1.0 

1.0 

1.4 

1.3 

1.3 

1.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.4 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

4.0 

0.5 

5.0 

5.0 

5.0 

5.0 

5.0 

5.5 

5.5 

5.5 

5.5 

5.5 

5.5 

5.5 

0.6 

6.0 

6.0 

6.0 

6.0 

6.0 

7.2 

7.2 

7.2 

7.2 

7.2 

7.2 

7.2 

0.7 

7.0 

7.0 

7.0 

7.0 

7.0 

9.2 

9.2 

9.2 

9.2 

9.2 

9.2 

9.2 

0.8 

8.0 

8.0 

8.0 

8.0 

8.0 

11.5 

11.5 

11.5 

11.5 

11.5 

11.5 

11.5 

0.9 

9.0 

9.0 

9.0 

9.0 

9.0 

14.0 

14.0 

14.0 

14.0 

14.0 

14.0 

14.0 

1.0 

10.0 

10.0 

10.0 

10.0 

10.0 

16.5 

16.5 

16.5 

16.5 

16.5 

16.5 

16.5 

3.6 

36.0 

36.0 

36.0 

36.0 

36.0 

81.5 

81.5 

81.5 

81.5 

81.5 

81.5 

81.5 
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Internal Mixing Acoustic Suppression 

Acoustic suppression of the internal jet noise is provided by perforate and bulk absorber 
liners. The effectiveness of the liners depends on the source location and the peak 
attenuation frequency of the liner. Numerically, the effective length of the perforate liner 
is defined as 

^ perf ^ perf ,end max (a s , L per f start ) (54) 

where L perf start is the length from the primary exit to the start of perforate treatment, 
and L perf end is the length from the primary exit to the end of perforate treatment. The 
source location x s is given by 


h . [>n 'k/ 


(55) 


where x d is the source location per unit nozzle diameter and is defined in table 3.10.4 as a 
function of the Strouhal number given by 

St = ^±Ei. (56) 

V.JM 

The length scale is the hydraulic diameter of the primary jet and the characteristic 
velocity is the internal mixing noise effective velocity previously defined. 

Similarly, the effective length of the perforate liner is 

^bulk ^ bulk . end max(x s , L bu i kstart ) (57) 


where L hldk start is the length from the primary exit to the start of the bulk treatment, and 
L hulk end is the length from the primary exit to the start of bulk treatment. 
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Table 3.10.4 Normalized Source Location for Internal Mixing Noise 


St 


St 


0.007147 

40.0 

0.768205 

3.5 

0.011104 

35.0 

0.924181 

3.0 

0.017610 

30.0 

1.002430 

2.8 

0.021333 

28.0 

1.092880 

2.6 

0.025976 

26.0 

1.198580 

2.4 

0.031817 

24.0 

1.323700 

2.2 

0.039245 

22.0 

1.474060 

2.0 

0.048810 

20.0 

1.658080 

1.8 

0.061320 

18.0 

1.888390 

1.6 

0.077999 

16.0 

2.026690 

1.5 

0.088468 

15.0 

2.184830 

1.4 

0.100790 

14.0 

2.367380 

1.3 

0.115416 

13.0 

2.580460 

1.2 

0.132952 

12.0 

2.832390 

1.1 

0.154224 

11.0 

3.134820 

1.0 

0.180389 

10.0 

3.504580 

0.9 

0.213125 

9.0 

3.966940 

0.8 

0.254950 

8.0 

4.561570 

0.7 

0.309823 

7.0 

5.354600 

0.6 

0.344058 

6.5 

6.465100 

0.5 

0.384350 

6.0 

8.131150 

0.4 

0.432363 

5.5 

10.908300 

0.3 

0.490427 

5.0 

16.463200 

0.2 

0.561909 

4.5 

33.129299 

0.1 

0.651856 

4.0 




Perforate and Bulk Liner Acoustic Suppression 

Acoustic suppression is modeled as the sum of perforate wall effects and bulk absorber 
effects. The total acoustic suppression is 

A^ sup = A dB perf + A dB hulk (58) 


Suppression provided by the perforate liner is given by 


MB Perf = [ AdBpelfpk L pelf F(e)/A 


Lperf <0.001 
Lperf a 0.001 


(59) 


where A dB f k is the peak attenuation per unit wave length, and L perf is the effective 
length of the perforate liner for a given frequency. The function f(o) is defined as 
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( ri30-i9A 

f(&) = min 

1.0, 0.9^ 10 > 


( 60 ) 


v J 

The wavelength A = c sec / / is based on the static temperature and sonic velocity of the 
expanded secondary flow. 

Suppression provided by the bulk liner is given by 


^dB bu i k - C bulk A dB bulk pk L bulk /(I + M mjx ) 


(61) 


where C bulk is obtained from table 3.10.5, and AdB hulk pk is the peak attenuation per unit 
length (default 1.6628). 


Table 3.10.5 Bulk Absorber Treatment Suppression for Internal Noise, C bulk 


L° 8 u)(f / f P k ) 

0 

90 

100 

110 

0 . deg 
120 

130 

140 

150 

180 

-3.6 

1.00000 

1.00000 

0.89100 

0.74100 

0.64600 

0.60300 

0.57500 

0.49000 

0.49000 

0.0 

1.00000 

1.00000 

0.89100 

0.74100 

0.64600 

0.60300 

0.57500 

0.49000 

0.49000 

0.1 

0.79400 

0.79400 

0.70800 

0.58900 

0.51300 

0.47900 

0.45700 

0.38900 

0.38900 

0.2 

0.63100 

0.63100 

0.56200 

0.46800 

0.40700 

0.38000 

0.36300 

0.30900 

0.30900 

0.3 

0.50100 

0.50100 

0.44700 

0.37200 

0.32400 

0.30200 

0.28800 

0.24500 

0.24500 

0.4 

0.39800 

0.39800 

0.35500 

0.29500 

0.25700 

0.24000 

0.22900 

0.19500 

0.19500 

0.5 

0.31600 

0.31600 

0.28200 

0.23400 

0.20400 

0.19100 

0.18200 

0.15500 

0.15500 

0.6 

0.25100 

0.25100 

0.22400 

0.18600 

0.16200 

0.15100 

0.14500 

0.12300 

0.12300 

0.7 

0.20000 

0.20000 

0.17800 

0.14800 

0.12900 

0.12000 

0.11500 

0.09770 

0.09770 

0.8 

0.15800 

0.15800 

0.14100 

0.11700 

0.10230 

0.09550 

0.09120 

0.07760 

0.07760 

0.9 

0.12600 

0.12600 

0.11200 

0.09330 

0.08130 

0.07590 

0.07240 

0.06170 

0.06170 

1.0 

0.10000 

0.10000 

0.08910 

0.07410 

0.06460 

0.06030 

0.05750 

0.04900 

0.04900 

1.1 

0.07940 

0.07940 

0.07080 

0.05890 

0.05130 

0.04790 

0.04570 

0.03980 

0.03980 

1.2 

0.06310 

0.06310 

0.05620 

0.04680 

0.04070 

0.03800 

0.03630 

0.03090 

0.03090 

1.3 

0.05010 

0.05010 

0.04470 

0.03720 

0.03240 

0.03020 

0.02880 

0.02450 

0.02450 

1.4 

0.03980 

0.03980 

0.03550 

0.02950 

0.02570 

0.02400 

0.02290 

0.01950 

0.01950 

1.5 

0.03160 

0.03160 

0.02820 

0.02340 

0.02040 

0.01910 

0.01820 

0.01550 

0.01550 

1.6 

0.02510 

0.02510 

0.02240 

0.01860 

0.01620 

0.01510 

0.01450 

0.01230 

0.01230 

1.7 

0.02000 

0.02000 

0.01780 

0.01480 

0.01290 

0.01200 

0.01150 

0.00977 

0.00977 

1.8 

0.01580 

0.01580 

0.01410 

0.01170 

0.01023 

0.00955 

0.00912 

0.00776 

0.00776 

3.6 

0.00025 

0.00025 

0.00022 

0.00019 

0.00016 

0.00015 

0.00014 

0.00012 

0.00012 
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Internal Shock Noise 


The free-field lossless sound pressure level for the internal shock noise component is 
given by 

SPLjs = Cjs + 10 Log m [A pri j rf )+ 10Log 10 [/^/(l + A.0frf s )_ 

+ 10Log 10 (l - D h pri /D eq pri )~ A hw + F l IS (84 )-F 2 js(0) (62) 

- F 3js(8)-40Log m (l-M a cos 8 D )-S IS (f,e,<p)- AclB„ p 

The internal shock noise coefficient C IS was found to be 8.5 dB higher for under- 
expanded jets than for over-expanded jets. C IS is given by 

[154.5 for M j pri <M d (Over - expanded) 

IS [163.0 for Mj pri >M d (Under - expanded) ' ' 


The area correlation factor for the internal shock noise is the primary jet flow area A pn . 


The shock cell strength parameter /3 IS is defined as the square root of the difference 

between the square of primary jet Mach number and the square of the design Mach 
number. If the jet is over-expanded the design Mach number is assumed to be one. The 
shock cell strength parameter is 


Pis ~ \ 


] M ), P n ~ 1 

' M lpri ~ M d 


for M j pn < 1.10 M d 
for M j pri > 1.10 M d 


(64) 


For non-circular nozzles, the ratio of the hydraulic diameter to the equivalent circular 
nozzle diameter is used to account for geometric effects. The hydraulic diameter D h pn is 
equal to four times the primary exit area divided by the perimeter. The equivalent 
diameter D eq pn is equal to the square root of 4 times the primary exit area divided by jt . 

For rectangular nozzles, D h pn / D eq pri can be expressed as a function of aspect ratio as 

yjjiAR /(I + AR ) . Here the aspect ratio is defined as the ratio of the nozzle width to the 
nozzle height. Figure 3.10.3 shows that increasing aspect ratio increases the internal 
shock noise. 
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Figure 3.10.3 Effect of aspect ratio of rectangular nozzles on internal shock noise 
Internal Shock Noise Directivity 

The internal shock noise directivity (in the absence of flight effects) is defined by three 
functions. Function F 1 IS combines polar and azimuthal effects and is given by 


„ [l l(#/l 20)L°g 10 (GF ) 

11 Lo gw (GF) 

F 2 is and F 3 is are explicit functions of the 
S; s -180-sin-'(l-M pr ,) 

F 2 is is given by 


for 6 < 120 
for 6 > 120 

Mach angle 6\ s defined by 


J 0 for 

F 2,is ( 0 ) - jo 25(0 _ g' [S ) for 

and F, IS is defined by the values given 


d<d' [S 

e>e'is 

in table 3.10.6. 


(65) 


( 66 ) 


(67) 
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Table 3.10.6 Internal Shock Noise Directivity Function 


e 

K, >1-1 M d ) 

KM 

sl.l M „) 

0 

23.7 

30.8 

10 

21.6 

27.2 

20 

19.6 

24.8 

30 

17.7 

21.8 

40 

15.9 

18.8 

50 

14.2 

15.8 

60 

12.6 

14.2 

70 

11.1 

12.7 

80 

9.6 

10.0 

90 

8.1 

8.1 

too 

8.5 

6.1 

110 

8.9 

5.6 

120 

9.6 

6.6 

130 

11.0 

7.6 

140 

13.0 

9.1 

150 

15.5 

10.6 

160 

19.5 

12.1 

170 

26.5 

13.6 

180 

34.5 

15.6 


Internal Shock Noise Flight Effects 


The effect of flight on the amplitude of the internal shock noise component is defined by 
the convective amplification term -40Log l0 (l - M a cos 0 D ). 


Internal Shock Noise Spectrum Directivity Relations 

The internal shock noise spectrum is correlated as a function of a modified Strouhal 
number, which includes azimuthal and aspect ratio effects. The length scale for the 
Strouhal number is the hydraulic diameter of the primary jet. The characteristic velocity 
is the internal shock noise is the primary jet velocity. The internal shock noise Strouhal 
number is defined as 


f'JN 

Log ,«(*«)- ftp (M - 1 )(l - M„ cos e D ) (68) 

and includes a Doppler shift. The internal shock noise peak Strouhal number varies 
significantly as a function of the Mach angle and whether the jet is over-expanded or 
under-expanded. Table 3.10.7 provided a listing of the logarithm of the peak Strouhal 
number for over-expanded and under-expanded jets. The internal shock noise spectrum 
function is obtained from table 3.10.8 as a function of Log l0 (St' PS ) given by 

L °g l0 {K ) = Lo § l0 {St is ) ~ Log w (St pk is ) (69) 
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Table 3.10.7 Internal Shock Noise Peak Strouhal Number 


0' 

£»*■„(** ,/s ) 

( MUSLIM ,) 

^ t-H 

A 

o 

too ^ 

O ^ 

0 

0.0 

-0.3 

10 

0.0 

-0.3 

20 

0.0 

-0.3 

30 

0.0 

-0.2 

40 

0.0 

-0.2 

50 

0.0 

-0.2 

60 

0.0 

-0.2 

70 

0.0 

-0.1 

80 

0.0 

-0.1 

90 

0.1 

0.0 

too 

0.1 

0.0 

110 

0.1 

0.1 

120 

0.2 

0.2 

130 

0.3 

0.3 

140 

0.4 

0.3 

150 

0.5 

0.3 

160 

0.6 

0.3 

170 

0.7 

0.3 

180 

0.8 

0.3 


Table 3.10.8 Internal Shock Noise Spectrum Function 


L °S\o( St is) 

s B 

-3.6 

176.0 

-1.0 

1.0 

0.1 

0.5 

3.6 

35.5 


Internal Shock Noise Suppression 

The effect of acoustic suppression on the internal shock noise is computed in the same 
manner as described for the internal mixing acoustic suppression. The source location 
is replaced with 



( 70 ) 


for the internal shock noise suppression analysis. 



Premerged Mixing Noise 


The free-field lossless sound pressure level for the premerged mixing noise component is 
given by 

(r s ,f,0,<ft) = C PM + 10 Log 10 (a pm / r s ) + 1 0 co PM Log 10 ( p PM j p a ) 

+ 75 Log 10 (v ePM /c a )+ F l pu (0,<p)-F 2 PM (0) 

~ 15Log 10 [(l + M cPM cos#) 2 + 0.04 M) pm ] 

-lOLog lo (l-M a cos0 D )-S PM (f,0) 


The premerged mixing noise coefficient C PM has a minimum value of 137 and a 
maximum value of 142 depending on the primary nozzle pressure ratio, the suppressor 
area ratio, and the lobe penetration. C PM is 


C PM = mini 


142, max 


137, 112.0 + 12.5 


t.pri , 


^SAR^PEN 


(72) 


The area correlation factor for the premerged mixing noise is the effective premerged 
areaA FM . A pu is a function of the degree of mixing, the mixer area ratio, the primary 
flow area, and the ejector exit area. A pu is given by 


A pm = (1 - DM ) — + (DM) A 

PM v ’mar v 7 e 


(73) 


The density correction factor for the premerged mixing region is 

3(V.™/0 55 


m pm 


.0.6 + (V,, P „/O“ 


1-1 


where V e PM is the effective premerged velocity given by 

K,i>» = V PM ~ a A 2 'pm ) C0S a j 

and the premerged jet velocity V PM is 

V PM .(l -DM)v f , r , + (DMy m „ 

The fully expanded velocity of the primary jet V fe pn is given by 


(74) 


(75) 


(76) 


V = M c 

V fe.pri m fe,pri C fe 


( 77 ) 
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The fully expanded velocity of the primary jet is obtained from the fully expanded Mach 
number 


M 


fe.pri 



and the sound speed at the exit of the mixer ejector 

= V Y priS^pri ^sje 


' fe 


The fully expanded static temperature is 


T r =T 

s,je A t,pri 


( p 

t,pri 

Lamb J 


- 1 +- 


Y pri 


(78) 


(79) 


(80) 


The density of the premerged mixing region is determined by the degree of mixing, the 
primary jet density, and the mixed jet density. The density of the premerged mixing 
region is 

p rM -(l-DM) PrH + (DM)p^ (81) 


Premerged Mixing Noise Directivity 

Two functions define the shape of the premerged mixing noise directivity in the absence 
of flight effects. F t PM {(), <p) defines the relationship between the polar and azimuthal 
angle effects. F lPM {(:), (/)) is given by 


-4.5(0' pm / 9O) 2 Log 10 {GF) 


l,PM 


(M- 


-4.5 


0.5 


®PM 

90 


+ 1.5 


^Pre 

120 


Log io {gf) 


-13.5(^ m /120 ) 2 Log w {GF) 


for 6 < 90 
for 90 <6 < 100 
for 6 > 100 


(82) 


where 


d' P M=d 


Yy A 0 ' 1 

y PM 

V C a J 


(83) 


F 2 PM (p) is a simple function of polar angle and is presented in table 3.10.9. It should be 
noted that the effects of source location have been removed from the function F l PM and 
f 

1 2,PM- 
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Table 3.10.9 Premerged Mixing Noise Directivity Function 


0 

F 2 ,Pm(P) 

6 

F 2 ,Pm(0) 

0 

8.7 

110 

7.2 

10 

8.7 

120 

6.7 

20 

8.7 

130 

9.6 

30 

8.7 

140 

11.5 

40 

8.7 

150 

15.7 

50 

8.7 

160 

18.4 

60 

8.7 

170 

19.4 

70 

8.7 

180 

20.4 

80 

8.7 

190 

21.4 

90 

8.7 

200 

22.4 

too 

7.7 

250 

27.4 


Premerged Mixing Noise Flight Effects 


The effects of source motion are defined by two functions. The first incorporates the 
effects of the source convection within the jet. The convection Mach number is defined 
as 


M c,pm =0.62(V pm - V a cos a , ) / c amb 


(84) 


The angle a, is the angle-of-attack of the jet. In the HSRNOISE code, the aircraft angle- 
of-attack is used. The second effect of flight on the amplitude of the premerged mixing 
noise component is defined by the term - \QnLog w (\ - M a cos 6 D ) . Note that n = 1 for 
the premerged mixing noise as compared to 4 for the internal mixing noise component. 

Premerged Mixing Noise Spectral Directivity 

The premerged mixing noise spectrum is correlated as a function of a modified Strouhal 
number, which includes temperature effects. The length scale for the Strouhal number is 
the hydraulic diameter of the primary jet. The characteristic velocity is the premerged 
mixing noise effective velocity. The internal mixing noise Strouhal number is defined as 


St - 

PM 




h,pri 


frp ^0.4(l + cos<9) 

t.pri 


V 


e,PM 


T 

V x amb J 


[\-M a cos 0 D ) 


(85) 


and includes a Doppler shift. The peak Strouhal number is a function of the Mach angle. 
For angles less than 110 degrees the peak remains constant. As the Mach angle 
increases, the peak shifts to lower Strouhal numbers. The peak Strouhal number for the 
premerged mixing noise component is presented in table 3.10.10. The premerged mixing 
noise spectrum function is presented in table 3.10.11 as a function of Log w (St' PM ) given 
by 

Log w (St PM )= Log l0 (St 

PM )- Logi o ( S t pk ,PM ) (86) 
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Table 3.10.10 Premerged Mixing Noise Peak Strouhal Number 


e 

L°gio(St pk pM ) 

e 

L°8io(§t p k,PM ) 

0 

0.1 

110 

0.1 

10 

0.1 

120 

0.0 

20 

0.1 

130 

- 0.1 

30 

0.1 

140 

- 0.2 

40 

0.1 

150 

- 0.3 

50 

0.1 

160 

- 0.4 

60 

0.1 

170 

- 0.5 

70 

0.1 

180 

- 0.5 

80 

0.1 

190 

- 0.6 

90 

0.1 

200 

- 0.6 

too 

0.1 

250 

- 0.7 


Table 3. 10. 1 1 Premerged Mixing Noise Spectrum Function 


L°g\o(St PM ) 

°PM 

- 3.6 

163.0 

- 0.9 

22.5 

- 0.8 

18.0 

- 0.7 

14.0 

- 0.6 

10.5 

- 0.5 

7.5 

- 0.4 

5.0 

- 0.3 

3.0 

- 0.2 

1.5 

- 0.1 

0.5 

0.0 

0.0 

0.1 

0.4 

0.2 

1.0 

0.3 

2.5 

3.6 

52.0 


Premerged Shock Noise 


The free-field lossless sound pressure level for the premerged shock noise component is 
given by 


SPL„(r„f,e) - C ps + 10 Log,„(A pu jr; )+ \0Log m \fi‘ K /( 1 + fi‘ K )J 
- ■F'i.ps (S) ~ r rs (e)-40Lo g Jl -M, cos9„ ) - S„ (/) 


The coefficient C PS is assigned a value of 157 for nozzle pressure ratios of 4 or greater. 
Nozzle pressure ratios less than 4 are computed by C PS = 11 + 20 (P tpri /P amb )• 

C ps - min[l57,77 + 2o(^„./.P„J] (88) 
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The area correlation parameter is the effective premerged area A PM . 
The shock cell strength parameter f> PS is defined as 


Pps = - 1 ( 89 ) 

2 

where the Mach number M . PS is defined as the product of the premerged effective Mach 

2 

number and the mixed jet Mach number. Mj PS is given by 


M 


2 

j,PS 


( V ) 

v PM 

(v 

\ c pm y 

V C mix y 


(90) 


where the premerged sound speed is 


'PM 




R T 

mix mix 1 s,PM 


(91) 


The static temperature of the premerged jet is 


T = 7 

1 s,PM 1 t.PM 


1 - 


1 


Y mix ) 


v; 




where the total pressure of the premerged region is 

t,,pm ={i-dm )r tpri + (dm )r tmix 


(92) 


(93) 


Premerged Shock Noise Directivity 

The directivity of the premerged shock noise is described by two functions. F 1PS is a 
function of the polar angle and the Mach angle 0 PS . The premerged directivity function 
is given by 


J 0 for 0 < d' PS 

F hPS (V) ~ j OJ5 (0 _ e > ps ) for q > Q' ps 


(94) 


where the Mach angle is defined in terms of the premerged jet Mach number as 

6' PS =180-sin- 1 (l-M iP J (95) 

The second directivity function is presented in table 3.10.12. 
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Table 3.10.12 Premerged Shock Noise Directivity Function 


9 

f 

1 2, PS 

9 

F 

1 2,PS 

0 

8.7 

too 

9.2 

to 

7.2 

110 

10.7 

20 

6.7 

120 

12.2 

30 

6.2 

130 

13.7 

40 

5.7 

140 

14.2 

50 

5.2 

150 

14.7 

60 

6.0 

160 

15.2 

70 

6.6 

170 

15.7 

80 

7.1 

180 

16.2 

90 

7.7 




No variations due to azimutal angle or aspect ratio are included in the prediction of the 
premerged shock noise. 


Premerged Shock Noise Flight Effects 


Correction to the premerged shock noise amplitude due to flight effects are given by the 
convection amplification term 40Log 10 (l- M cos<9 D ). 


Premerged Shock Noise Spectrum Directivity Relations 

The premerged shock noise spectrum is correlated as a function of a modified Strouhal 
number, which includes the strength of the shock cells and the premerged jet velocity. 
The length scale for the Strouhal number is the equivalent diameter of the primary jet. 
The characteristic velocity is the premerged jet velocity. The premerged shock noise 
Strouhal number is defined as 

St K -(iZ>„,,/0.7V p „K s (l-M„cos^)x 

l/il + MWc J COS tf} 2 +0.0196(V FM / C J 2 

and includes a Doppler shift. The premerged shock noise peak Strouhal number is given 
by 


LogioSt pk PS 


f- 0.3 9 <50 

[-0.2 50 < # < 180 


(97) 


The premerged shock noise spectrum function is presented in table 3.10.13 as a function 
of Log w (St' PS ) given by 

Log io (.St ps ) = Log 10 (St ps ) — Log 10 (st pk PS ) (98) 
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Table 3.10.13 Premerged Shock Noise Spectrum Function 


1 o ( St ps ) 

Sps 

- 3.6 

138.5 

- 0.2 

5.9 

- 0.1 

2.0 

0.0 

0.0 

1.0 

1.0 

3.6 

36.0 


Merged Mixing Noise 

The free-field lossless sound pressure level for the merged mixing noise component is 


SPLmM i r s F MM + lQL°8lo[ A e j / r s )+ 1 ® 0) MM P°8\ 0 {Pmix / Pamb ) 

+ 15Lo 8loi V MM / c a)- F l{ 0 ’<P)- F 2{&) 


- 15Log 10 [(l + M cMM cos #) 2 + 0.04.M 

- lOLogiofl - M a cos 0 D ) - S MM ( f,0 ) 


2 

c.MM 


(99) 


The merged mixing noise coefficient C MM is given by 


C mm = min 


142.5, max 


138, 131.7 + 6.0 


P ■ IP u ^ 

?,pn / amb 

SARy/PEN 


( 100 ) 


C MM has a minimum value of 138 and a maximum value of 142.5. The primary nozzle 
pressure ratio, the suppressor area ratio, and the lobe penetration determine the actual 
value of the coefficient. 

The density correction for the merged mixing region is 




3(K,mm/cJ 


3.5 


0-6 +(K,mm/0 


3.5 


1-1 


( 101 ) 


where V e MM is the effective premerged velocity given by 

v .,«« -r.WMr./r.Jcosa, (102) 

The angle a r is the angle-of-attack of the jet. In the HSRNOISE code, the aircraft angle- 
of-attack is used. 
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Merged Mixing Noise Directivity 


Two functions define the shape of the merged mixing noise directivity in the absence of 
flight effects. F l MM {(), (j>) defines the polar directivity and azimuthal directivity effects. 

^1 ,MM (e,t) is given by 


1 MM 


where 


0 

< Ud' M -120) " 
8 


Log 10 {GF) 


for 6' m <120 
for 0' M >120 


(103) 


0L = o 




v o.i 


V C a J 


(104) 


The function F, MM {0) is obtained from table 3.10.14. The original formulation for the 
independent variable 6 in reference 14 is a function of the source location and may 
obtain values greater than 180 degrees. Entries in table 3.10.14 greater than 180 degrees 
can be ignored. 


Table 3.10.14 Merged Mixing Noise Directivity Function 


0 

1 2, MM 

0 

F 

l 2 ,MM 

0 

10.9 

110 

4.4 

10 

10.5 

120 

3.2 

20 

10.1 

130 

2.0 

30 

9.7 

140 

3.2 

40 

9.3 

150 

3.9 

50 

8.9 

160 

10.9 

60 

8.5 

170 

17.9 

70 

8.1 

180 

24.9 

80 

7.9 

190 

31.9 

90 

7.7 

200 

38.9 

too 

4.9 

250 

73.9 


Merged Mixing Noise Flight Effects 


The effects of source motion on the merged mixing noise component is defined by two 
functions. The first incorporates the effects of the source convection within the jet. The 
convection Mach number is defined as 


M 


c,MM 


= 0.62(V ( 


e.MM 


-V 


cos a,)/ c 


amb 


(105) 


The second fight effect term is the convection amplification to the first power. 
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Merged Mixing Noise Spectrum Directivity Relations 

The merged mixing noise spectrum is correlated as a function of a modified Strouhal 
number, which includes temperature effects. The length scale for the Strouhal number is 
the equivalent diameter of the primary jet. The characteristic velocity is the merged 
mixing noise effective velocity. The Strouhal number for the merged mixing region is 
given by 


St = 

L> ‘ MM 


fD, 


ej,eq 


f J ^ 0.4(l +COS0' 


V 


e,MM 


T 

V ± amb J 


(l-M fl costf D ) 


(106) 


The peak Strouhal number St pk is obtained from table 3.10.15 as a function of the 
directivity angle. For angles less than 70 degrees the peak remains constant. As the 
Mach angle increases, the peak shifts to higher Strouhal numbers. The peak Strouhal 
number for the merged mixing noise component is presented in table 3.10.15. The 
merged mixing noise spectrum function is presented in table 3.10.16 as a function 
of Log l0 (St' MM ) given by 

Log l0 (St MM ) = Log^St MM )~ L°g m {Stp kMM ) ( 107 ) 


Table 3.10.15 Merged Mixing Noise Peak Strouhal Number 


0 

L°g\o(st p k .mm) 

0 

L°g\{)(St p k MM ) 

0 

0.3 

110 

0.5 

to 

0.3 

120 

0.5 

20 

0.3 

130 

0.5 

30 

0.3 

140 

0.5 

40 

0.3 

150 

0.6 

50 

0.3 

160 

0.6 

60 

0.3 

170 

0.7 

70 

0.3 

180 

0.7 

80 

0.4 

190 

0.8 

90 

0.4 

200 

0.8 

too 

0.4 

250 

1.0 
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Table 3.10.16 Merged Mixing Noise Spectrum Function S MM (f,0 ) 


Log(St' MM ) Directivity angle, 0 degrees 



0 

120 

130 

250 

- 3.6 

157.8 

157.8 

157.8 

157.8 

- 1.0 

25.2 

25.2 

25.2 

25.2 

- 0.9 

20.8 

20.8 

20.8 

20.8 

- 0.8 

16.8 

16.8 

16.8 

16.8 

- 0.7 

13.2 

13.2 

13.2 

13.2 

- 0.6 

10.0 

10.0 

10.0 

10.0 

- 0.5 

7.2 

7.2 

7.2 

7.2 

- 0.4 

4.8 

4.8 

4.8 

4.8 

- 0.3 

3.0 

3.0 

3.0 

3.0 

- 0.2 

1.6 

1.6 

1.6 

1.6 

- 0.1 

0.6 

0.6 

0.6 

0.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

0.8 

0.8 

0.9 

0.9 

0.2 

2.4 

2.4 

2.7 

2.7 

0.3 

4.4 

4.4 

5.0 

5.0 

0.4 

6.8 

6.8 

7.8 

7.8 

0.5 

9.6 

9.6 

11.4 

11.4 

0.6 

12.8 

12.8 

14.6 

14.6 

0.7 

16.4 

16.4 

18.6 

18.6 

0.8 

20.4 

20.4 

23.0 

23.0 

0.9 

24.8 

24.8 

27.4 

27.4 

3.6 

143.6 

143.6 

146.2 

146.2 


Mixer Ejector Propagation Effects 


Propagation effects for the method described in this section include ambient pressure, 
atmospheric attenuation, ground reflection attenuation, and extra ground attenuation 
(EGA). For convenience, source modification for the number of sources is also included 
with the propagation effects. The free-field lossless spectra is modified for propagation 
effectes by 

SPL = SPLj ree _j ield i oss j ess + A SPL (108) 


where AS PL is given by 


A SPL = 20 Log 


10 


Pa(z,Ya(z s ) 

p a {zoVdz 0 ), 


~ \bs r s + G - EGA + 10 Log 10 (n eng ) 


(109) 


Additional details concerning propagation effects can be found in section 3.13. 
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3.11 Airframe Noise 


The airframe noise prediction method adopted for the HSR program is based on the 
method developed by Fink [15]. Airframe noise is produced by elements of the wing, 
fuselage, and landing gear. The noise generating mechanisms are divided into the 
following eight components: 

1. Wing Noise 

2. Leading Edge Noise 

3. Horizontal Tail Noise 

4. Vertical Tail Noise 

5. Inboard Flap Noise 

6. Outboard Flap Noise 

7. Main Gear Noise 

8. Nose Gear Noise 

Fink developed empirical functions for each noise component based on flyover data from 
a wide variety of aircraft including a sail plane, general aviation aircraft, a business jet, a 
Boeing 747 commercial transport, and a delta wing Convair F106B. Since Fink's method 
was developed from primarily conventional subsonic airplane designs, a model scale 
airframe noise test was conducted in the Boeing Fow-speed Aeroacoustic Facility. The 
objective was to provide estimates of the high-speed civil transport (HSCT) airframe 
noise levels. Noise data were collected on a full span 3% scale model of a double delta 
HSCT configuration. From this data, spectral corrections for the Fink method were 
developed. 

Symbols 

A, in wing inboard trailing edge flap area, ft' ( input parameter - AFIN) 

A f out wing outboard trailing edge flap area, ft 2 ( input parameter - AFOUT) 

A h horizontal tail area, ft' ( input parameter - AH) 

A, vertical tail area, ft' ( input parameter - AV) 

A w wing area, ft 2 ( input parameter - AW) 

b f in inboard flap span, ft ( input parameter - BFIN) 

b f out outboard flap span, ft ( input parameter - BFOUT) 

b h horizontal tail span, ft ( input parameter - BH) 

b v vertical tail span, ft ( input parameter - BV) 

b w wing span, ft ( input parameter - BW) 

c a ambient sound speed, ft/sec ( input parameter - CA) 

d mg wheel diameter main gear, ft ( input parameter - DMG) 

d ng wheel diameter nose gear, ft ( input parameter - DNG) 

D directivity function 

/ one-third-octave band center frequencies, Hz 

F spectrum function 
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Sc 

G 

k 

L 

mg 

L ng 


N. 


mg 


N, 


ng 


N 


w,mg 


K 

N 


f.in 

Nf,out 
P ref 


P 

r s 

* 

r 

K 

St 

T a 

P f An 

P /.out 

6 * 

<P 

Pa 

n* 

e D 

Pa 


2 

gravitational constant, 32.174 ft/sec" 
general function 
constant coefficient 

length of main gear strut, ft ( input parameter - RLMG) 

length of nose gear strut, ft ( input parameter - RLMG) 

aircraft Mach number ( input parameter - MA) 

number of main gear ( input parameter - NMG) 

number of nose gear ( input parameter - NMG) 

number of wheels per main gear ( input parameter - NWMG) 

number of wheels per nose gear ( input parameter - NWMG) 

number of inboard training edge slots ( input parameter - NFS IN) 

number of outboard training edge slots ( input parameter - NFSOUT) 

reference pressure, 4.177 x 10'^ (lb/ft 2 ) 
nondimensional mean-square acoustic pressure 
source radius, r s = 150 ft 

& 

nondimensional source radius, r = r s / b w 
Reynolds number 
Strouhal number 
ambient temperature, °F 

inboard flap deflection angle, degrees ( input parameter - FLAIN) 
outboard flap deflection angle, degrees ( input parameter - FLAOUT) 

nondimensional turbulent boundary thickness 
azimuthal angle, degrees 
ambient dynamic viscosity, lbm/ft-sec 
nondimensional acoustic power 
Doppler angle, degrees 

ambient density (lbm/ft ), ( input parameter - RHOA) 


Method 

The airframe prediction equations for each component are cast in the following general 
form. The nondimensional mean-square pressure (p 1 ) is 

{ p 2 )' -£^ d (»o,<pMsi) (d 

where II*is the nondimensional acoustic power, D is the directivity function, and F is 
the spectrum function. The general form of the nondimensional acoustic power is 


no 



rf = kM n G 


( 2 ) 


where k is constant, M" a is the aircraft Mach number raised to a power, and G is 
uniquely defined for each airframe source. G is a function of the boundary layer 
thickness for airfoil related noise sources such as wing and flap noise. For the landing 
gear sources, G is a function of the geometric properties of the landing gear. 

The free-field lossless sound pressure level for a given component is obtained from 

2 

SPL = 10 LogJp 2 )* + 20 Log l0 ^ (3) 

Pref 


Each airframe noise component is predicted at a source radius of 150 feet. The total 
airframe noise is the sum o f the mean- square pressures of the individual components. 

Wing Noise 

The primary noise generating mechanism of the wing is the convection of the turbulent 
boundary layer past the trailing edge. The acoustic power is given by 

n* = kMX (4 ) 


[4.464x10 5 CleanWing 
[7.075 xl0“ 6 SuperCleanWing 


For most commercial transports including the HSCT, k 


The nondimensional turbulent boundary layer thickness 
equation developed by Prandtl. 


Sl -0.31 

b w v 


PqM„C a \. 

PaK J 


, - 0.2 


The ambient dynamic viscosity is given by 
T 15 xl0“ 8 


fx a = 2.279 


(r, +198.6) 


= 7.075 x 10 -6 . 

is computed from the flat plate 

(6) 

(7) 


where the ambient temperature T a is 


ill 




The directivity function is assumed to be aligned with the lift dipole and is given by 

D(d D ,<p) = 4 cos 2 (^) cos 2 / 2) (9) 

The spectrum function is an empirical function of the Strouhal number. The spectrum 
function for a delta wing is given by 

F{St w ) = 0A&5(l0Stjll0StJ s + 0.5} 4 (10) 

where the length scale for the Strouhal number is based on the wing span. 

( 11 ) 

M aC a 

Leading Edge Slat 

Leading edge slat noise is caused by two mechanisms. First, the slat produces an 
increment of wing trailing edge noise due to the impact on the boundary layer of the 
wing. Second, the leading edge slat itself produces trailing edge noise. The equations for 
the acoustic power, boundary layer thickness, and directivity function are similar to the 
equations for wing noise. The constant k for the leading edge slat was derived for a 


rectangular wing. 


ir = kM s x 

(12) 

k = 4.464 xlCT 5 

(13) 

r \-°- 2 

0 3 n A "’ P “ M « C ° A " 

w bl{ f*Jb w J 

(14) 

D(d D ,<p) = 4 cos 2 (</)) cos 2 (0 D / 2) 

(15) 


As a result of the two noise generating mechanisms, two spectrum functions F l and F 2 
are required. 


F^St J = 0.485(l0St J 4 [(lOSr J 35 + O. 5} 4 
F 2 (StJ = 0.613(2.19StJ 4 [(2.19StJ 1 ' 5 +0.5^ 
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These spectrum functions were empirically derived based on the wing Strouhal number 


St 


w 


4^[l-M„cos(0j] 


(18) 


Total leading edge slat noise is the sum of the mean-square pressures from the two noise 
generating mechanisms 


2 



n 

4 jit* 2 


d(o d ,Af,(s w ) + F 2 (sJ 


(19) 


Horizontal and Vertical Tail 


The noise generated by the horizontal and vertical tails is also similar to the noise 
generated by the wing. The noise generating mechanism is the same, e.g. the convection 
of the turbulent boundary layer past the trailing edge. Since the equations for the 
horizontal and vertical tails are identical except for the directivity functions, only the 
formulation for the horizontal tail is presented. To obtain the formulation for the vertical 
tail, replace the subscript h with v in the following equations. 

n* = kM 5 X (20) 


k = 4.464 xl0“ 5 


( 21 ) 


K - 0 . 37 ^ 

b,. 


s \- 0.2 

PMA A 

f A 


( 22 ) 


D(0 D ,t)- 


|4cos 2 (^)cos 2 (^ £) ll) 
|4sin 2 (^)cos 2 (6' £) / 2) 


Horizontal Tail 
Vertical Tail 


F(St h ) = 0.613(105^ ) 4 [{lOSt h ) L5 + 0.5 


(23) 

(24) 


s '*-77 

M r 


(25) 


Inboard and Outboard Flap Noise 

The noise generated by the inboard and outboard flaps is also the result of the convection 
of the turbulent boundary layer past the trailing edge. Since the equations for the inboard 
and outboard flaps are identical, only the formulation for the inboard flap is presented. 
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To obtain the formulation for the outboard flap, replace the subscript fin with fout in the 
following equations. 


n* = kM 6 a G 

(26) 

f 2.787 xl0“ 4 Single or Double Slotted Flaps 

k = < , 

(27) 

[3.509x10 4 Triple Slotted Flaps 




G= , ' sin (/7 /m ) 

U w 

(28) 

D{d D ,(p)= 3(sin 6 f in cos d D + cos d f in sin d D cos (p) 2 

(29) 


The spectrum function for single or double slotted flaps is given by 



^0.0480 St fM 
0.1406 St?£ 5 
216.49 Stf M 


St fjn < 2 

2 - St f in < 20 
Stf M >20 


where the Strouhal number is 




f A fM 
M c b f 

ci ci j ,in 


[l-M a cos(0 fl )_ 


The spectrum function for triple slotted flaps is given by 



' 0.0257 St fjn 
0.0536 Stffr 05 


St f, in K 2 

2 - St f in <75 
Stfjn >75 


(30) 


(31) 


(32) 


The slot configuration for the HSCT was single slotted flaps. 

Main Gear and Nose Gear Noise 

The noise generating mechanism from the gear is complex and highly dependent on the 
landing gear design. Gear noise is dominated by two mechanisms; noise generated by the 
wheels and noise generated by the strut. Separate formulations for wheel and strut noise 
are provided. 
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Since the equations for the main gear and the nose gear are identical, only the formulation 
for the main gear is presented. To obtain the formulation for the nose gear, replace the 
subscript mg with ng in the following equations. 

Wheel Noise 

The acoustic power, directivity, and spectrum functions for wheel noise are given by 


n* = kM 6 a G 
k = 4.349 xl0“ 4 


D (^^) = | sin2 (^) 

13.59*, ; s ( 12.5 + 

- 7T^[ ] - M :Cosl&„)] 

" M aC a 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


Strut Noise 

The acoustic power, directivity, and spectrum functions for strut noise are given by 


n* = kMlG 

(39) 

k = 2.753 xl0“ 4 

(40) 

(d V L 

q mg mg 

V b w J d mg 

(41) 

d{ 6 d ,(/) ) = 3sin 2 (^ 0 )sin 2 (^) 

(42) 

f(*„J-1.280S4(i.06 + s4) J 

(43) 


Total gear noise is the sum of the mean-square pressures from the wheels and the strut 

(44) 


e.g. (p 2 ) = (p 2 ) + (p 2 

° ' / \ / wheel \ I strut 
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or 



n 

Am, 


*2 


{@D ’ P wheels (pt mg ) + D strut i^D ’ *P strut (*^( mg ). 


(45) 


Spectral Adjustment for HSCT Configuration 

Acoustic data were collected in the Boeing Low-Speed Aeroacoustic Facility (LSAF) on 
a 3% full span model of a double delta HSCT wing configuration. The model measured 
approximately 106 inches long with a wing span of approximately 45 inches. The 
inboard wing leading edges were swept approximately 68 degrees and the outboard 
leading edges were swept approximately 48 degrees. The adjustable leading edge flaps 
were deflected 30 degrees, and the adjustable trailing edges were deflected 20 degrees. 
Removable landing gear were used on the model. The main gear was modeled as a 
triple-post strut with six wheels per truck. The nose gear was modeled as a single-post 
with two-wheels. There were no gear cavities or doors, and the model had flow-through 
nacelles. 


Boeing made noise predictions using the Fink method. At the 90 degree angle a set of 
spectral deltas were derived to improve the agreement between prediction and 
measurement. A set of directivity deltas were derived based on the PNLT directivity 
shape. The spectral deltas are computed by 

A = A(0)B(f) (46) 

The deltas are then added to the sound pressure level 

SPLcorrected = SPL Fink A (47) 


Table 3.11.1 provides the decibel corrections that are applied to wing, leading edge flaps, 
horizontal tail, vertical tail, inboard flap, and outboard flap. Corrections were not derived 
for the main gear or the nose gear. 

Total Noise 

The free-field lossless total airframe noise is the sum of the mean-square pressures of the 
eight airframe subcomponents. Propagation effects are defined in section 3.13. 
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Table 3.1 1.1 Airframe Noise Adjustment Functions for the HSCT Configuration, (dB) 


Emission Angle, degrees 




10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

Freq . 

A ( 0 ) 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.0 

0.7 

0.5 

0.3 

0.0 

- 0.2 

- 0.5 

- 0.8 

- 1.0 

Hz . 

B ( f ) 


















50 

- 12.0 

- 16.8 

- 16.8 

- 16.8 

- 16.8 

- 16.8 

- 16.8 

- 16.8 

- 16.8 

- 12.0 

- 8.4 

- 6.0 

- 3.6 

0.0 

2.4 

6.0 

9.6 

12.0 

63 

- 11.0 

- 15.4 

- 15.4 

- 15.4 

- 15.4 

- 15.4 

- 15.4 

- 15.4 

- 15.4 

- 11.0 

- 7.7 

- 5.5 

- 3.3 

0.0 

2.2 

5.5 

8.8 

11.0 

80 

- 10.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 10.0 

- 7.0 

- 5.0 

- 3.0 

0.0 

2.0 

5.0 

8.0 

10.0 

too 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

ill’ 

4.5 

7.2 

9.0 

125 

- 8.0 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 8.0 

- 5.6 

- 4.0 

- 2.4 

0.0 

1.6 

4.0 

6.4 

8.0 

160 

- 8.0 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 8.0 

- 5.6 

- 4.0 

- 2.4 

0.0 

1.6 

4.0 

6.4 

8.0 

200 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 


4.5 

7.2 

9.0 

250 

- 9.5 

- 13.3 

- 13.3 

- 13.3 

- 13.3 

- 13.3 

- 13.3 

- 13.3 

- 13.3 

- 9.5 

- 6.7 

- 4.8 

- 2.9 

0.0 

1.9 

4.8 

7.6 

9.5 

315 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

400 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

500 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

630 

- 10.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 10.0 

- 7.0 

- 5.0 

- 3.0 

0.0 

2.0 

5.0 

8.0 

10.0 

800 

- 10.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 10.0 

- 7.0 

- 5.0 

- 3.0 

0.0 

2.0 

5.0 

8.0 

10.0 

1000 

- 10.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 14.0 

- 10.0 

- 7.0 

- 5.0 

- 3.0 

0.0 

2.0 

5.0 

8.0 

10.0 

1250 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

1600 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

2000 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

2500 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

3150 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

4000 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

5000 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

6300 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

8000 

- 9.0 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 12.6 

- 9.0 

- 6.3 

- 4.5 

- 2.7 

0.0 

KB 

4.5 

7.2 

9.0 

10K 

- 8.0 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 11.2 

- 8.0 

- 5.6 

- 4.0 

- 2.4 

0.0 

1.6 

4.0 

6.4 

8.0 
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3.12 Jet Shielding 


A change in the received noise spectra caused by the presence of adjacent jets is known as jet 
shielding. The method described in the section is a simplified method that was adopted for the 
HSR system studies. The method requires only geometric input and does not require input 
concerning the exhaust parameters of the adjacent jets or the spacing between the jets. Jet 
shielding can be applied to all aft radiated noise source. 


Symbols 


h 

N 

eng 

SI 


y a 

z a 

P 

<P 

A SPL 


aircraft altitude at visual overhead, ft ( input parameter - ALTE VO ) 

number of engines ( input parameter - NENG) 

sideline distance, ft ( input parameter - SL) 

aircraft position along the flight path, ft 

aircraft position perpendicular to flight path, ft (y„ = 0 ) 

aircraft position above ground level, ft ( input parameter - Z A) 

angle between the flight path vector and the vector for the aircraft to the origin 

flight path angle, degrees (j input parameter - FPA) 

jet shielding suppression level, dB 


The empirical equation which models jet shielding is given by 


ASPL = 10 Log l0 1 



where 


( 1 ) 


cos 2 /? 


{x a cos (p + z a Sin 4>) 2 

2 , 2 
x„ + z„ 


and the aircraft position along the flight path is 


( 2 ) 


X 


a 


Zq ~ k 
tan^ 


(3) 


Equation 1 is only valid for values of f> between 0° and 90°. No shielding effect is 
applied for angles greater than 90°. Note that ASPL = 0 for the case of a single engine 
aircraft ( N eng = 1) and for an observer directly under the flight path (SI = 0). 
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3.13 Propagation Effects 

The fan, core, turbine, and airframe noise sources are modeled as point sources. The 
free-field spectra for each is computed at a reference distance of 150 feet and then 
propagated to the observer assuming straight line propagation. Propagation effects 
include spherical spreading, convective amplification, ambient pressure, atmospheric 
attenuation, ground reflection attenuation, and extra ground attenuation (EGA). For 
convenience, source modification for the number of sources is also included with the 
propagation effects. The free-field lossless spectra is modified for propagation effectes 
by 


SPL SPLj ree _j- ieId i oss i ess + A SPL 


( 1 ) 


where A SPL is given by 


A SPL = 20 Log 


^150 A 


10 


V 7 


-40 Log l0 (l - M a cos d D ) + 10 Log l0 (n ) 


+ 20 Log 


10 


Pa(z s )c 2 a(z s ) 


“ \ bs r s +G - EGA 


vA,b 0 k 2 U J, 

and where each adjustment term is defined in table 3.13.1. 

Table 3.13.1 - Definition of A SPL Adjustments 


Adjustment for: 


Term 


( 2 ) 


Spherical spreading 

onr f 150 1 

20 Log I0 

{ P J 

Convective amplification 

40Log w (l - M a cos0 D ) 

Number of sources 

lOiogdAU) 

Ambient pressure 

on? f PafckkJl 

20 L°g l0 a ) * 

\PaV-oFaV-o)) 

Atmospheric absorption 

\ A absU pi ,rh)\r s 

f = f 

p 1 - M a cos 6 p 

Ground reflections 

G 

Extra Ground Attenuation 

EGA 


Atmospheric absorption is computed by either the SAE ARP866A method or the Shields 
and Bass method [16]. Ground effects include attenuation of the spectra due to 
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reflections from the surface and extra ground attenuation due to temperature and wind 
gradients. Ground reflection is predicted by the method developed by Chien and Soroka 
[17], EGA is predicted by SAE AIR 923 [18], 

The jet noise source is modeled as a distributed source. Special propagation effects are 
applied to the jet noise source to account for acoustic and geometric near field effects. 
These effects are described in detail in the sections 3.9 and 3.10. 

3.14 Noise Metrics 

One-third-octave band noise levels are predicted at a user defined location, typically one 
of the FAA certification points. Output from each noise source is provided for the range 
of frequencies and emission angles requested by the user. HSRNOISE computes the 
following noise metrics; OASPL, PNL, PNLT, and EPNL. 
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4.0 Code Verification 


The HSRNOISE code was compared with noise levels predicted by the MCP code for the 
Technology Concept (TC) aircraft. The TC aircraft is a Mach 2.4 supersonic transport 
with a takeoff gross weight of 650,000 pounds and a range of 5000 nautical miles. The 
propulsion system is a mixed flow turbofan with a two-dimensional bifurcated inlet and a 
mixer ejector suppressor nozzle. The input parameters for this analysis are from the 
3870.47 engine cycle. 

One-third-octave band spectra, OASPL, PNL, PNLT, and EPNL noise predictions were 
made at the FAA approach, sideline, and centerline observer locations. Section 4. 1 
describes the analysis of the EPNL comparisons as well as comparisons of the OASPL, 
PNL, and PNLT noise metrics. Differences between the MCP and HSRNOISE 
predictions are explained. Section 4.2 describes the analysis of the one-third-octave band 
spectra. Free-field lossless and full system spectra are presented. The free-field lossless 
spectra identify differences in the implementation of the noise source methodologies. 

The full system predictions, which include atmospheric attenuation, ground reflections, 
extra ground attenuation, and jet shielding, identify differences in propagation effects. 
The input files used for the free-field lossless and full system predictions are provided in 
section 5.0. 

4.1 EPNL Comparisons 

This section presents comparisons of the EPNL levels predicted by the HSRNOISE and 
MCP codes. The EPNL noise metric was instrumental in determining the weight and 
range of the TC aircraft. Columns A, B, and C of table 4.1.1 show the results of the 
initial EPNL comparisons between the HSRNOISE and MCP codes. Differences in total 
aircraft EPNL of 2.1 at approach, 2.6 on the sideline and 1.0 on the centerline were not 
anticipated. Causes for the differences included different numerical techniques for 
integrating the PNLT time history, different procedures for computing the PNLT time 
history, different ground reflection models, and differences in source levels. 

EPNL Numerical Integration 

The EPNL noise metric is calculated by integrating the PNLT time history between the 
limits of the 10 dB down points. FAA FAR 36 Appendix B dictates that the EPNL 
should be calculated using half-second time intervals. Both the HSRNOISE and MCP 
codes predict the PNLT time history at fixed emission angles, which do not correspond to 
half-second time intervals. Since half-second PNLT data is not available, the question 
becomes what numerical integration technique is appropriate. Numerical integration was 
initially performed in the HSRNOISE code using the FAA procedure or equivalently the 
rectangular rule. Half-second PNLT data was obtained by interpolating the PNLT time 
history. The interpolation procedure was designed to ensure that the maximum PNLT 
was included in the interpolated time history. Numerical integration is performed in the 
MCP code using the trapezoidal rule. 
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To determine if the integration technique caused differences in EPNL levels, both the 
trapezoidal rule and the rectangular rule were implemented in the HSRNOISE EPNL 
algorithm. Column B of table 4.1.1 lists the EPNL levels predicted by the rectangular 
rule. These predictions are identified as method 1. Column D lists the EPNL levels 
predicted using the trapezoidal rule. These predictions are identified as method 2. 
Column F lists the difference between methods 1 and 2. Except for the sources with 
tones, both integration methods predict essentially the same levels. The trapezoidal rule 
was adopted as the EPNL integration technique for the HSRNOISE code because it does 
not require interpolation of the spectra. Interpolating spectra with tones may introduce 
errors in the half second interval data required for the rectangular rule integration. 

Table 4.1.1 EPNL Comparisons between MCP and HSRNOISE Codes 



(A) 

MCP 

(Boeing) 

(B) 

HSRNOISE 
(Method 1) 

(C) 

A1 

(D) 

HSRNOISE 
(Method 2) 

(E) 

A2 

(F) 

A3 

Approach 

Fan 

91.8 

92.2 

0.4 

92.7 

0.9 

-0.5 

Core 

84.4 

84.3 

-0.1 

84.1 

-0.3 

0.2 

Turbine 

94.5 

97.6 

3.1 

97.1 

2.6 

0.5 

Jet 

92.6 

95.0 

2.4 

95.0 

2.4 

0.0 

Airframe 

89.4 

90.8 

1.4 

90.8 

1.4 

0.0 

Total 

99.7 

101.8 

2.1 

101.8 

2.1 

0.0 

Sideline 

Fan 

89.1 

89.0 

-0.1 

89.2 

-0.1 

-0.2 

Core 

74.3 

76.7 

2.4 

76.8 

2.5 

-0.1 

Turbine 

81.3 

85.4 

4.1 

85.3 

4.0 

0.1 

Jet 

100.6 

103.6 

3.0 

103.5 

2.9 

0.1 

Total 

101.3 

103.9 

2.6 

103.8 

2.5 

0.1 

Centerline 

Fan 

87.4 

87.2 

-0.2 

87.6 

0.2 

-0.4 

Core 

79.0 

78.3 

-0.7 

78.3 

-0.7 

0.0 

Turbine 

85.6 

88.9 

3.3 

88.7 

3.1 

0.2 

Jet 

99.4 

100.4 

1.0 

100.2 

0.8 

0.2 

Total 

100.2 

101.2 

1.0 

101.0 

0.8 

0.2 


A1 = B-A A2 = D-A A3 = B-D 


Method 1: EPNL calculated using rectangular rule Method 2: EPNL calculated using trapezoidal rule 


PNLT Calculations 

OASPL, PNL, and PNLT time history comparisons are provided in figures 4.1.1 through 
4.1.9. The PNLT comparisons identified a second difference between HSRNOISE and 
MCP. After discussing the PNLT comparisons with Boeing, it was learned that tone 
corrections below 1250 Hertz and SPL levels below 30 dB were not included the MCP 
predictions. This modification to the PNLT calculations was not incorporated into the 
HSRNOISE code. 
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Ground Reflection Model 

The HSRNOISE and MCP codes use different methods to model the effects of ground 
reflections. The HSRNOISE code uses the Chien and Soroka model and the MCP code 
uses the Shivashankara model. All input parameter to the ground reflection models were 
explicitly defined with the exception of the ground flow resistance. Modeling the ground 
as a soft surface as opposed to a hard reflecting surface has a major impact on the shape 
of the received spectra and the EPNL level. A range of ground flow resistances were 
tested in the HSRNOISE code to determine the impact on the EPNL level. The EPNL 
levels listed in table 4.1.1 were calculated using a ground flow resistance of l.lxlO 6 
lbm/sec-ft , which is representative of an asphalt surface. The EPNL levels listed in table 
4.1.2 were calculated using a ground flow resistance of 9.4xl0 3 lbm/sec-ft 3 , which is 
representative of a grassy field. Interestingly, the value of the ground flow resistance, 
which gave the best EPNL comparison, did not yield the best spectral comparison. 
Figures 4.1.10 and 4.1.11 show the total aircraft spectra at the sideline location with a 
ground flow resistance of l.lxlO 6 lbm/sec-ft 3 and 9.4xl0 3 lbm/sec-ft 3 respectively. The 
spectral shapes predicted by HSRNOISE and MCP are significantly different in figure 
4.1.10, especially at the approach angle of 30 degrees. For this reason the ground flow 
resistance was set to l.lxlO 6 lbm/sec-ft 3 for the spectral analysis. 

Table 4.1.2 EPNL Comparisons between MCP and HSRNOISE Codes 
Using Ground Flow Resistance Typical of a Grassy Field 



(A) 

MCP 

(Boeing) 

(B) 

HSRNOISE 
(Method 2) 

(C) 

A 

Approach 

Fan 

91.8 

90.1 

1.7 

Core 

84.4 

82.9 

1.5 

Turbine 

94.5 

95.8 

-1.3 

Jet 

92.6 

93.1 

-0.5 

Airframe 

89.4 

89.9 

-0.5 

Total 

99.7 

100.0 

-0.3 

Sideline 

Fan 

89.1 

87.0 

2.1 

Core 

74.3 

74.9 

-0.6 

Turbine 

81.3 

82.7 

-1.4 

Jet 

100.6 

101.6 

-1.0 

Total 

101.3 

102.0 

-0.7 

Centerline 

Fan 

87.4 

85.7 

1.7 

Core 

79.0 

77.2 

1.8 

Turbine 

85.6 

87.0 

-1.4 

Jet 

99.4 

99.0 

0.4 

Total 

100.2 

99.6 

0.6 


A = A-B 


Method 2: EPNL calculated using trapezoidal rule 
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Figure 4.1.4 Sideline OASPL comparisons for the TC aircraft 
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Figure 4. 1 .7 Centerline OASPL comparisons for the TC aircraft 








Figure 4.1.10 Total sideline noise prediction for the TC aircraft 
Ground flow resistance 1.1 x 10 6 lbm/sec-ft 3 
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Figure 4.1.11 Total sideline noise prediction for the TC aircraft 

3 3 

Ground flow resistance 9.4 x 10 lbm/sec-ft 
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4.2 Spectral Analysis 

The previous section provided a comparison of the EPNL noise levels for each noise 
source. However, as was demonstrated in the previous section, good EPNL comparisons 
do not necessarily imply good spectral comparisons. This section provides comparisons 
of the component spectra. Free-field predictions were made to remove the effects of 
ground reflection, EGA, and jet shielding. Corrections to the source spectra or "Deltas" 
were also removed from the free-field predictions. After the HSRNOISE code was 
verified with the free-field comparison, full system predictions were made. Free-field 
and full system predictions were made for the approach, sideline, and centerline 
conditions. Input files to the HSRNOISE code used in the spectral analysis are listed in 
Section 5. 

One-third-octave band spectra were plotted at emission angles from 30 to 170 degrees in 
20 degree increments. Each spectrum was divided into three regions as defined below. 

Low Frequencies 50 - 250 Hertz 
Mid Frequencies 315 - 1600 Hertz 
High Frequencies 2000 - 10000 Hertz 

Each region of the spectra was evaluated according to the following criteria 

Excellent 0.0 < |A| < 0.3 

Good 0.3 <|A| < 0.7 

Fair 0.7 < |A| < 1.0 

Poor |A| >1.0 

where A is the RMS error for a given frequency range given by 

( 1 ) 

In equation 1, A is the frequency index. N equals 8 for the low, mid, and high frequency 
evaluations. N equals 24 when all frequencies are evaluated. 



136 



4.2.1 Total Aircraft Noise 


Figures 4.2. 1.1, 4.2. 1.3, and 4.2. 1.5 show the total aircraft free-field spectra for the 
approach, sideline, and centerline observer locations. Figures 4.2. 1.2, 4.2. 1.4, and 4.2. 1.6 
show comparisons of the full system predictions. RMS error values are provided in 
tables 4.2. 1.1 through 4.2. 1.6. Discrepancies between MCP and HSRNOISE predictions 
are apparent in both the free-field and full system spectra. 

Table 4.1.1 (page 123) listed a difference of 2.1 EPNdB in the total aircraft noise for the 
approach case. The free-field spectra for this case is shown in figure 4.2. 1.1. From the 
component analysis it was determined that fan tones dominate the high frequency spectra 
at emission angles of 30 and 50 degrees and turbine noise dominates the high frequency 
spectra for emission angles greater than 90 degrees. While the MCP and HSRNOISE 
codes use the same theoretical models to predict fan and turbine noise, the 
implementation of the methods are different. Refer to Sections 4.2.3 and 4.2.4 for details. 
Differences on the order of 1 to 1.5 dB were detected in the fan tone levels. Differences 
on the order of 1 to 5 dB were detected in the high frequency region of the free-field 
turbine spectra. Figure 4.2. 1.2 shows the full system spectra for the approach case. In 
addition to the discrepancies observed in the free-field spectra, discrepancies are also 
noted in the full system spectra at the frequencies associated with ground dips. 

Better EPNL comparisons were achieved for the sideline and centerline cases where fan 
and turbine noise are less important. Comparisons of the free-field spectra were excellent 
in the low frequency range and mostly excellent to good in the mid frequency range. As 
in the approach case, comparison of the high frequency spectra was poor. Discrepancies 
in the high frequency range are associated with the ground reflection model which will be 
explained in greater detail in the jet noise comparisons. 

Table 4.2.1. 1 Free-field Total Aircraft Approach Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

» 

Low 

0.3 

0.3 

0.3 

0.1 

0.1 

0.1 

0.3 

0.3 

Mid 

0.5 

0.6 

0.6 

0.3 

0.8 

0.9 

0.6 

0.6 

High 

0.7 

0.8 

1.4 

2.0 

3.4 

4.1 

0.8 

1.4 

All 

0.5 

0.6 

0.9 

1.1 

2.0 

2.4 

0.6 

0.9 

Table 4.2. 1.2 Full System Total Aircraft Approach Noise RMS Error Table, dB 

Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

170° 

Low 

0.9 

1.4 

1.8 

0.9 

0.9 

1.6 

1.4 

1.8 

Mid 

0.5 

1.2 

0.9 

1.2 

1.6 

1.4 

1.2 

0.9 

High 

0.9 

0.9 

2.2 

3.2 

4.2 

3.7 

0.9 

2.2 

All 

0.8 

1.2 

1.7 

2.1 

2.6 

2.4 

1.2 

1.7 
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Table 4.2. 1.3 Free-field Total Aircraft Sideline Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

o 

O 

I'' 

Low 

0.0 

0.1 

0.1 

0.0 

0.1 

0.0 

0.1 

0.1 

Mid 

0.1 

0.1 

0.1 

0.0 

0.0 

0.1 

0.1 

0.1 

High 

1.9 

1.0 

0.8 

0.9 

1.5 

1.4 

1.0 

0.8 

All 

1.1 

0.6 

0.5 

0.5 

0.9 

0.8 

0.6 

0.5 

Table 4.2. 1.4 Full System Total Aircraft Sideline Noise RMS Error Table, dB 

Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

170° 

Low 

0.1 

0.3 

0.9 

0.9 

1.5 

1.3 

0.3 

0.9 

Mid 

1.7 

1.1 

0.9 

0.9 

1.2 

1.6 

1.1 

0.9 

High 

5.7 

0.3 

0.9 

1.3 

2.4 

2.8 

0.3 

0.9 

All 

3.4 

0.7 

0.9 

1.1 

1.8 

2.0 

0.7 

0.9 


Table 4.2. 1.5 Free-field Total Aircraft Centerline Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

» 

Low 

0.1 

0.1 

0.1 

0.1 

0.3 

0.0 

0.1 

0.1 

Mid 

0.5 

0.8 

0.4 

0.1 

0.3 

0.3 

0.8 

0.4 

High 

1.7 

1.8 

1.6 

2.5 

3.0 

3.1 

1.8 

1.6 

All 

1.0 

1.1 

1.0 

1.5 

1.8 

1.8 

1.1 

1.0 


Table 4.2. 1.6 Full System Total Aircraft Centerline Noise RMS Error Table. dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

o 

o 

110° 

130° 

150° 

170° 

Low 

1.2 

1.3 

2.1 

1.3 

1.3 

2.2 

1.3 

2.1 

Mid 

1.0 

0.8 

0.6 

0.6 

0.6 

0.7 

0.8 

0.6 

High 

0.9 

0.8 

1.3 

1.8 

2.0 

1.8 

0.8 

1.3 

All 

1.0 

1.0 

1.5 

1.3 

1.4 

1.7 

1.0 

1.5 
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4.2.2 Jet Noise 


Jet noise dominated the spectra for the sideline and centerline cases. Figures 4.2.2. 1, 

4.2. 2. 3, and 4.2. 2. 5 show the free-field one-third-octave band spectra for the approach, 
sideline, and centerline cases. Figures 4.2. 2. 2, 4.2. 2.4, and 4.2. 2. 6 show the full system 
predictions. RMS error values are provided in tables 4.2.2. 1 through 4.2. 2. 6. 
Comparisons of the free-field spectra were excellent for all frequency ranges. 
Comparisons of the full system spectra were good to fair in the low frequency range, fair 
to poor in the mid frequency range, and generally poor in the high frequency range. 

The approach and centerline predictions provide insight into the cause of the 
discrepancies between the HSRNOISE and MCP predictions. The free-field and full 
system predictions differ only by the execution of the ground reflection model. 
HSRNOISE uses the Chien and Soroka ground reflection model while MCP uses the 
Shivashankara model. Although the locations of the ground dips between the two 
methods match very well, SPL differences as large as 4 dB were noted at frequencies 
corresponding to the ground dips. The ground flow resistance parameter was adjusted in 
an attempt to reduce the discrepancies. The default of SIGMA= l.lxlO 6 lbm/sec-ft 3 
provided the best comparison. 

In addition to differences in SPL levels at the ground dips, the Shivashankara model 
attenuates high frequencies more than the Chien and Soroka model. In general, 
differences on the order of 2 dB were noted at all directivity angles for frequencies 
greater than 1000 Hertz. Since the PNLT noise metric weights the high frequency 
spectra, the HSRNOISE code would be expected to predict higher EPNL levels than the 
MCP code. 

The full system sideline predictions include extra ground attenuation (EGA) and jet 
shielding in addition to ground reflections. The methods used to predict EGA and jet 
shielding are identical to those used in MCP. Differences in the SPL are consistent with 
those for the approach and sideline cases. 

It should be noted that the jet noise was predicted using the Boeing JN8C4 methodology. 
Comparisons were not made using Stone’s jet noise methodology because it was a last 
minute addition to the HSRNOISE code. 
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Table 4.2.2. 1 Free-field Approach Jet Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

o 

O 

» 

Low 

0.1 

0.1 

0.2 

0.2 

0.2 

0.2 

0.1 

0.2 

Mid 

0.2 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

High 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

All 

0.2 

0.2 

0.2 

0.2 

0.3 

0.3 

0.2 

0.2 

'able 4.2. 2. 2 Full System Approach Jet Noise RMS Error Table, dl 

Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

170° 

Low 

0.8 

1.3 

1.6 

0.9 

0.8 

1.5 

1.3 

1.6 

Mid 

0.6 

1.3 

1.4 

1.3 

1.4 

1.2 

1.3 

1.4 

High 

1.7 

2.1 

2.2 

2.3 

2.3 

2.2 

2.1 

2.2 

All 

1.1 

1.6 

1.8 

1.6 

1.6 

1.7 

1.6 

1.8 


Table 4.2. 2. 3 Free-field Sideline Jet Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

o 

» 

Low 

0.0 

0.0 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

Mid 

0.1 

0.0 

0.0 

0.1 

0.0 

0.1 

0.0 

0.0 

High 

0.2 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

All 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 


Table 4.2. 2. 4 Full System Sideline Jet Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

o 

o 

110° 

130° 

150° 

170° 

Low 

0.1 

0.3 

0.9 

0.9 

1.5 

1.3 

0.3 

0.9 

Mid 

1.7 

1.1 

1.0 

0.9 

1.2 

1.5 

1.1 

1.0 

High 

1.9 

0.5 

1.0 

1.3 

2.1 

2.3 

0.5 

1.0 

All 

1.5 

0.7 

1.0 

1.1 

1.7 

1.7 

0.7 

1.0 
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Table 4.2. 2. 5 Free-field Centerline Jet Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

t- 

Low 

0.0 

0.0 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

Mid 

0.1 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

High 

0.3 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

All 

0.2 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

able 4.2. 2. 6 Full System Centerline Jet Noise RMS Error Table, d 

Emission Angle 
Frequencies 

30° 

50° 

o 

o 

» 

o 

o 

110° 

130° 

150° 

170° 

Low 

0.1 

0.3 

0.9 

0.9 

1.5 

1.3 

0.3 

0.9 

Mid 

1.7 

1.1 

1.0 

0.9 

1.2 

1.5 

1.1 

1.0 

High 

1.9 

0.5 

1.0 

1.3 

2.1 

2.3 

0.5 

1.0 

All 

1.5 

0.7 

1.0 

1.1 

1.7 

1.7 

0.7 

1.0 
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4.2.3 Turbine Noise 


Turbine noise dominated the high frequency spectra during approach for emission angles 
greater than 90 degrees. Figures 4.2.3. 1, 4.2. 3. 3, and 4.2. 3. 5 show the free-field one- 
third-octave band spectra for the approach, sideline, and centerline cases. Figures 
4.2. 3. 2, 4.2. 3.4, and 4.2. 3. 6 show the full system predictions. RMS error values are 
provided in tables 4.2.3. 1 through 4.2. 3. 6. These spectra (both free-field and full system) 
were rated as poor for all frequency ranges. 

Following discussions with engineers at Boeing, three causes for the discrepancies in the 
turbine predictions were identified. First, convective amplification was accounted for 
twice in the MCP turbine module. This causes the MCP spectra to be high for angles less 
than 90 degrees and low for angles greater than 90 degrees. Second, the emission angle 
rather than the Doppler angle was used to Doppler shift the MCP turbine spectra. 

Finally, the algorithm used to predict the turbine noise required a correction to remove 
atmospheric attenuation built into the method. This correction was discovered too late in 
the HSR program to be incorporated into the MCP predictions. 

Spectral differences between 3 and 9 dB were noted in the high frequency range of the 
full system turbine spectra. Implementation differences identified in the free-field 
predictions combined with the differences noted in the ground reflection model discussed 
in the section 4.2.2 caused the EPNL levels predicted by HSRNOISE to be 3 EPNdB 
higher than the levels predicted by MCP. 
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Table 4.2.3. 1 Free-field Approach Turbine Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

o 

O 

Low 

1.6 

1.3 

0.8 

1.0 

0.3 

0.4 

1.3 

0.8 

Mid 

1.8 

1.5 

0.9 

1.2 

0.5 

0.3 

1.5 

0.9 

High 

3.7 

3.1 

3.1 

2.6 

3.0 

3.7 

3.1 

3.1 

All 

2.6 

2.1 

1.9 

1.8 

1.8 

2.2 

2.1 

1.9 

Table 4.2. 3. 2 Full System Approach Turbine Noise RMS Error Table, dB 

Emission Angle 

30° 

50° 

o 

O 

» 

o 

o 

110° 

130° 

150° 

170° 

Frequencies 









Low 

3.0 

2.8 

2.6 

0.9 

1.1 

1.7 

2.8 

2.6 

Mid 

2.5 

2.1 

0.9 

1.5 

2.2 

2.5 

2.1 

0.9 

High 

1.0 

0.8 

2.3 

3.4 

4.9 

6.5 

0.8 

2.3 

All 

2.4 

2.1 

2.1 

2.2 

3.2 

4.1 

2.1 

2.1 

Table 4.2. 3. 3 Free-field Sideline Turbine Noise RMS Error Table, dB 

Emission Angle 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

» 

Frequencies 









Low 

3.6 

2.8 

1.8 

0.1 

1.0 

2.8 

2.8 

1.8 

Mid 

3.5 

2.7 

1.7 

0.3 

1.2 

2.9 

2.7 

1.7 

High 

5.8 

4.1 

3.8 

4.5 

6.4 

8.1 

4.1 

3.8 

All 

4.4 

3.3 

2.6 

2.6 

3.8 

5.2 

3.3 

2.6 

Table 4.2. 3. 4 Full System Sideline Turbine Noise RMS Error Table, dB 

Emission Angle 

30° 

50° 

o 

O 

o 

o 

110° 

130° 

150° 

170° 

Frequencies 









Low 

3.6 

2.6 

2.2 

0.8 

1.9 

3.5 

2.6 

2.2 

Mid 

4.7 

3.4 

1.9 

1.1 

2.2 

4.1 

3.4 

1.9 

High 

12.6 

1.5 

1.5 

3.1 

5.6 

7.2 

1.5 

1.5 

All 

8.0 

2.6 

1.9 

1.9 

3.6 

5.2 

2.6 

1.9 
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Table 4.2. 3. 5 Free-field Centerline Turbine Noise RMS Error Table, dB 


Emission Anele 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

t- 

Low 

4.4 

3.5 

1.5 

0.4 

1.5 

2.2 

3.5 

1.5 

Mid 

4.3 

3.4 

1.3 

0.3 

1.6 

2.4 

3.4 

1.3 

High 

5.1 

3.5 

2.7 

4.0 

5.4 

7.6 

3.5 

2.7 

All 

4.6 

3.4 

1.9 

2.3 

3.4 

4.8 

3.4 

1.9 

Table 4.2. 3. 6 Full System Centerline Turbine Noise RMS Error Table, dB 

Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

o 

o 

110° 

130° 

150° 

170° 

Low 

4.6 

3.8 

2.6 

1.1 

1.4 

2.2 

3.8 

2.6 

Mid 

4.2 

3.0 

0.8 

1.0 

2.6 

3.3 

3.0 

0.8 

High 

5.0 

1.2 

1.9 

3.7 

5.4 

7.2 

1.2 

1.9 

All 

4.6 

2.9 

1.9 

2.3 

3.6 

4.8 

2.9 

1.9 
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4.2.4 Fan Noise 


The high frequency spectrum of the fan noise component was a dominant source during 
approach for emission angles less than 90 degrees. Figures 4.2.4. 1, 4.2.4.3, and 4.2. 4.5 
show the free-field one-third-octave band spectra for the approach, sideline, and 
centerline cases. Figures 4.2.4.2, 4.2. 4.4, and 4.2. 4.6 show comparisons of the full 
system predictions. RMS error values are provided in tables 4.2.4. 1 through 4.2.4.6. 

Discrepancies on the order of 0.5 to 3 dB were noted in the free-field spectra. The largest 
discrepancies occurred in the low and high frequency ranges. Discrepancies on the order 
of 0.3 to 0.5 dB were noted in the mid frequency range. One implementation difference 
between HSRNOISE and MCP was identified. The broadband spectrum index is a 
function of the peak frequency. HSRNOISE uses the peak frequency as defined by 
equation 6 on section 3.6. MCP uses the center frequency of the one-third-octave band in 
which the peak frequency occurs. This implementation difference only affects the 
broadband component. 

Discrepancies were also noted in the free-field tone levels. This discrepancy is most 
evident in the free-field centerline spectra (figure 4.2. 4.5). At 30 degrees, the 
fundamental tone occurs at 4000 Hertz, and is approximately 20 dB above the broadband 
component. A difference of 1.5 dB was noted between the HSRNOISE and MCP 
predicted levels. Unfortunately, the cause of this discrepancy was not resolved. 

Table 4.2.4. 1 Free-field Approach Fan Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

170° 

Low 

0.9 

0.5 

0.9 

1.0 

1.0 

1.1 

0.5 

0.9 

Mid 

0.4 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

0.3 

High 

0.7 

1.1 

1.3 

1.3 

1.3 

1.3 

1.1 

1.3 

All 

0.7 

0.7 

0.9 

0.9 

1.0 

1.0 

0.7 

0.9 

able 4.2.4.2 Full System Approach Fan Noise RMS Error Table, d 

Emission Angle 
Frequencies 

30° 

50° 

70° 

90° 

110° 

130° 

150° 

170° 

Low 

0.9 

1.5 

2.1 

3.5 

3.9 

5.5 

1.5 

2.1 

Mid 

0.6 

1.1 

0.8 

0.8 

1.0 

0.9 

1.1 

0.8 

High 

2.3 

1.8 

1.3 

0.9 

1.0 

1.4 

1.8 

1.3 

All 

1.5 

1.5 

1.5 

2.1 

2.4 

3.3 

1.5 

1.5 
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Table 4.2. 4.3 Free-field Sideline Fan Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

170° 

Low 

1.0 

0.9 

1.0 

0.9 

0.9 

1.0 

0.9 

1.0 

Mid 

0.3 

0.4 

0.4 

0.4 

0.5 

0.5 

0.4 

0.4 

High 

2.9 

2.0 

1.3 

1.2 

1.1 

1.5 

2.0 

1.3 

All 

1.8 

1.3 

0.9 

0.9 

0.9 

1.1 

1.3 

0.9 

"able 4.2. 4.4 Full System Sideline Fan Noise RMS Error Table, dl 

Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

170° 

Low 

1.3 

2.0 

2.0 

1.8 

1.9 

2.7 

2.0 

2.0 

Mid 

2.1 

1.5 

1.1 

1.2 

1.3 

3.5 

1.5 

1.1 

High 

6.1 

2.1 

1.4 

1.0 

1.2 

3.0 

2.1 

1.4 

All 

3.8 

1.9 

1.5 

1.4 

1.5 

3.1 

1.9 

1.5 


Table 4.2.4.5 Free-field Centerline Fan Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

» 

Low 

1.0 

0.9 

1.0 

0.9 

0.9 

1.0 

0.9 

1.0 

Mid 

0.3 

0.4 

0.4 

0.4 

0.5 

0.5 

0.4 

0.4 

High 

2.9 

2.0 

1.3 

1.2 

1.1 

1.5 

2.0 

1.3 

All 

1.8 

1.3 

0.9 

0.9 

0.9 

1.1 

1.3 

0.9 


Table 4.2.4.6 Full System Centerline Fan Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

o 

o 

110° 

130° 

150° 

170° 

Low 

3.1 

3.6 

5.2 

4.1 

4.0 

5.2 

3.6 

5.2 

Mid 

1.4 

1.1 

0.8 

0.7 

0.6 

0.7 

1.1 

0.8 

High 

1.6 

1.1 

1.0 

0.9 

0.7 

0.8 

1.1 

1.0 

All 

2.2 

2.3 

3.1 

2.5 

2.4 

3.1 

2.3 

3.1 
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4.2.5 Airframe Noise 


Airframe noise was only predicted for the approach case. Figure 4.2.5 shows the free- 
field one-third-octave band spectra for the approach, sideline, and centerline cases. 
Figure 4.2. 5. 2 shows the full system prediction. RMS error values are provided in tables 
4.2.5. 1 and 4.2.5.2. 

Comparisons of the free-field spectrum were excellent. Comparisons of the full system 
spectrum were excellent in the low frequency region except for frequencies 
corresponding to the ground dips. Comparisons of the mid-frequency region were good 
to fair. The spectral match was poor in the high frequency region. The discrepancy in 
the high frequency spectra is a result of the ground reflection model. As discussed in the 
section on jet noise, the Shivashankara model attenuates the high frequencies more than 
the Chien and Soroka model. 


Table 4.2.5. 1 Free-field Approach Airframe Noise RMS Error Table, dB 


Emission Angle 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

C- 

Low 

0.3 

0.3 

0.3 

0.3 

0.2 

0.2 

0.3 

0.3 

Mid 

0.3 

0.2 

0.3 

0.3 

0.3 

0.2 

0.2 

0.3 

High 

1.5 

0.8 

0.6 

0.5 

0.5 

0.6 

0.8 

0.6 

All 

0.9 

0.5 

0.4 

0.4 

0.4 

0.4 

0.5 

0.4 

Table 4.2. 5. 2 Full System Approach Airframe Noise RMS Error Table, dB 

Emission Angle 
Frequencies 

30° 

50° 

o 

O 

» 

90° 

110° 

130° 

150° 

170° 

Low 

0.9 

1.5 

1.8 

1.0 

1.0 

1.7 

1.5 

1.8 

Mid 

0.5 

1.1 

1.0 

0.8 

0.9 

0.8 

1.1 

1.0 

High 

1.1 

1.6 

1.7 

1.8 

1.8 

1.7 

1.6 

1.7 

All 

0.9 

1.4 

1.5 

1.3 

1.3 

1.5 

1.4 

1.5 
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4.2.6 Core Noise 


Core noise did not contribute to the total noise of the TC aircraft for the approach, 
sideline, or centerline cases. Figures 4.2.6. 1, 4.2. 6. 3, and 4.2. 6. 5 show the free-field one- 
third octave band spectra for the approach, sideline, and centerline cases. Figures 4.2. 6. 2, 
4.2. 6. 4, and 4.2. 6. 6 show the full system predictions. RMS error values are provided in 
tables 4.2.6. 1 through 4.2. 6. 6. 


An examination of figures 4.2.6. 1 through 4.2. 6. 6 shows a shift in the spectrum predicted 
by the MCP and HSRNOISE programs. This shift is not due to the Doppler effect. In 
fact, the core noise spectrum is not Doppler shifted in either the MCP or the HSRNOISE 
codes. This shift is caused by different implementations of the spectrum function. The 
core noise spectrum function is tabulated as a function of the nondimensional parameter 
?/ defined as 


>1 = 10Tog 10 



l f J 

^ J peak ' 


(Eq. 6 Section 3.7) 


In the MCP code rj is converted to an integer before determining the spectrum level. 
Since the spectrum function is tabulated for integer values, interpolation is not required. 
In the HSRNOISE code, r/ is treated as a real variable and consequently, interpolation of 
the spectrum function is required. This relatively minor implementation difference 
caused the low and high frequency ranges to differ by more than one dB. 


In some cases, the MCP core spectra did not contain SPL levels for the 24 one-third 
octave bands as can be seen in the figure 4.2.6. 1. This caused very large RMS errors to 
be tabulated in the low and high frequency ranges. Because core noise was not a 
significant noise source, this discrepancy was not resolved. 


As with the other noise sources, differences in level were noted at frequencies 
corresponding to the ground dips. The poor match in the high frequency range observed 
in the other sources was not as prevalent in the core comparisons because the shift noted 
in the free-field comparisons tended to mask the increased attenuation of the ground 
reflection model. 
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Table 4.2.6. 1 Free-field Approach Core Noise RMS Error Table, dB 


Emission Angle 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

» 

Frequencies 









Low 

43.1 

42.4 

27.1 

9.6 

0.6 

1.4 

42.4 

27.1 

Mid 

0.5 

0.4 

0.9 

0.6 

0.3 

0.8 

0.4 

0.9 

High 

0.9 

0.8 

0.6 

1.8 

8.6 

25.5 

0.8 

0.6 

All 

24.9 

24.5 

15.7 

5.6 

5.0 

14.7 

24.5 

15.7 

Table 4.2. 6. 2 Full System Approach Core Noise RMS Error Table, dB 

Emission Angle 

30° 

50° 

70° 

90° 

110° 

130° 

150° 

170° 

Frequencies 









Low 

47.1 

42.5 

27.3 

9.8 

1.2 

2.5 

42.5 

27.3 

Mid 

0.5 

1.3 

1.4 

0.7 

1.0 

1.4 

1.3 

1.4 

High 

0.9 

0.8 

1.4 

0.6 

9.2 

25.9 

0.8 

1.4 

All 

27.2 

24.5 

15.8 

5.7 

5.4 

15.0 

24.5 

15.8 

Table 4.2. 6. 3 Free-field Sideline Core Noise RMS Error Table, dB 

Emission Angle 

30° 

50° 

70° 

90° 

110° 

130° 

150° 

170° 

Frequencies 









Low 

48.8 

44.7 

32.4 

5.2 

1.3 

0.2 

44.7 

32.4 

Mid 

1.0 

0.5 

0.8 

0.5 

0.6 

0.8 

0.5 

0.8 

High 

5.1 

1.9 

1.5 

1.1 

22.6 

31.7 

1.9 

1.5 

All 

28.3 

25.9 

18.7 

3.1 

13.1 

18.3 

25.9 

18.7 

Table 4.2. 6.4 Full System Sideline Core Noise RMS Error Table, dB 

Emission Angle 

30° 

50° 

70° 

90° 

110° 

130° 

150° 

170° 

Frequencies 









Low 

51.1 

48.6 

32.5 

5.3 

1.9 

1.4 

48.6 

32.5 

Mid 

2.0 

1.3 

1.1 

0.9 

1.2 

1.2 

1.3 

1.1 

High 

2.8 

1.7 

0.8 

0.9 

21.5 

32.5 

1.7 

0.8 

All 

29.6 

28.1 

18.8 

3.2 

12.5 

18.8 

28.1 

18.8 
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Table 4.2.6. 5 Free-field Centerline Core Noise RMS Error Table, dB 


Emission Anele 
Frequencies 

30° 

50° 

o 

O 

90° 

110° 

130° 

150° 

o 

O 

» 

Low 

50.4 

46.1 

39.7 

12.9 

1.2 

0.3 

46.1 

39.7 

Mid 

1.0 

0.9 

0.7 

0.9 

0.6 

1.1 

0.9 

0.7 

High 

4.4 

1.5 

0.9 

1.3 

14.5 

36.0 

1.5 

0.9 

All 

29.2 

26.6 

22.9 

7.5 

8.4 

20.8 

26.6 

22.9 

Table 4.2. 6. 6 Full System Centerline Core Noise RMS Error Table, dB 

Emission Angle 
Frequencies 

30° 

50° 

70° 

90° 

110° 

130° 

150° 

170° 

Low 

53.8 

49.4 

39.8 

13.1 

1.7 

1.7 

49.4 

39.8 

Mid 

1.2 

0.7 

0.8 

0.5 

0.9 

0.5 

0.7 

0.8 

High 

0.8 

1.0 

0.8 

0.7 

14.0 

36.6 

1.0 

0.8 

All 

31.1 

28.5 

23.0 

7.6 

8.1 

21.2 

28.5 

23.0 
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4.3 Conclusions 


Comparisons have been presented for the SPL, OASPL, PNL, PNLT, and EPNL noise 
levels predicted by the HSRNOISE and MCP codes for the Technology Concept (TC) 
aircraft. It was anticipated that the HSRNOISE would match the EPNLs predicted by the 
MCP within a range of +/- 0.25. Unfortunately, this criterion was not met. Several 
philosophical design differences between the two codes were identified which explain the 
differences in EPNL levels. 

One very significant difference between the two codes is the way in which the PNLT 
noise metric is calculated. In the MCP predictions, tone corrections are not made for 
frequencies below 1250 Hertz and for SPL levels below 30 dB. This modification to the 
PNLT time history resulted in MCP predicting lower EPNL levels than HSRNOISE. 

HSRNOISE and MCP use different ground reflection models. HSRNOISE uses the 
Chien and Soroka ground reflection model; MCP uses the Shivashankara model. The 
Chien and Soroka model predicts deeper ground dips and attenuates the high frequency 
region less than the Shivashankara model. The deeper ground dips caused the PNLT 
correction to be higher in the HSRNOISE code, which also contributed to the higher 
EPNL levels. Of course, modifying the PNLT algorithm to the one used in the MCP 
code would reduce some of this discrepancy. 

Differences were also identified in the implementation of the fan, core, and turbine 
source models. In the case of the fan and core modules, the differences were related to 
spectrum index function. MCP relates the spectrum index to the center frequency before 
calculation of the SPL. HSRNOISE relates the spectrum index to the center frequency 
after calculating the SPL. This causes differences in the magnitude of the predicted SPL. 
The MCP core module did not always predict SPL levels for the 24 one-third-octave 
bands. This causes the PNLT correction to be higher than expected when the spectrum is 
truncated in the high frequency range. An error was identified in the MCP turbine 
module, which had the effect of lowering the EPNL level during approach. It should be 
noted that very good agreement was obtained for the free-field predictions of the jet and 
airframe components. 

The HSRNOISE code was developed to provide members of the HSR program with a 
code to continue research on the HSCT. This report provides a User's Manual, fully 
documents the noise prediction algorithms, and provides comparisons of predicted noise 
levels for the TC aircraft. While the HSRNOISE and MCP codes did not predict 
identical EPNL levels as was anticipated, the main causes for the discrepancies have been 
identified. Because of the termination of the HSR program discrepancies between the 
two codes were not fully resolved. 
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5.0 Input Files 


This section presents the input files used in the verification of the HSRNOISE code. 
Free-field and full system files are listed for the approach, sideline, and centerline cases. 
Free-field input files are included in section 5.1. Full system input files are included in 
section 5.2. 

5.1 Free-field 

The free-field input files include atmospheric absorption but do not include ground 
reflections, EGA or jet shielding. Note that the airframe and jet noise source include 
spectral correction. Spectral corrections were removed from the fan, core, and turbine 
noise sources. 

5.1.1 Free-field Approach 

# 

# Approach test case for the HSRNOISE code 

# 

# Free-field predictions 
$GEOM 


HMIC 

= 

o 

o 


f 

SL 

= 

o 

o 



$END 





$WEATHER 





PSIAMB 

= 

14.7 


r 

RH 

= 

70.0 


r 

TFAMB 

= 

77.0 



$END 





$FLIPATH 





ALTEVO 

= 

394.0 


r 

FPA 

= 

-3.0 


f 

PAE 

= 

7 . 1 


/ 

VAIR 

= 

267.7168 


$END 





$ SOURCE 





IAFRAM 

= 

T, 



IATM 

= 

1, 



IEGA 

= 

F, 



IGRND 

= 

F, 



I FAN 

= 

T, 



ICORE 

= 

T, 



I JET 

= 

T, 



ITUR 

= 

T, 



NENG 

= 

1 



$END 





$DELTAS 





F AND IBB 

= 

17*0.0 

r 


FANSPBB 

= 

24*0.0 

r 


FANBPF 

= 

o 

o 

r 


FAN2BPF 

= 

o 

o 

r 


FAN3BPF 

= 

o 

o 

r 


FAN4BPF 

= 

o 

o 

r 


CORDIBB 

= 

17*0.0 , 
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CORSPBB = 

24*0 

.0 , 





TURDIBB = 

17*0 

.0 , 





TURSPBB = 

24*0 

.0 , 





TURDITN = 

17*0 

.0 , 





TURSPTN = 

24*0 

.0 , 





AIRFDICW 

= 

8*1.4 

, 1.0, 

0.7, 0.5, 0.3, 0.0, -0.2, 

-0.5, 

-0.8, 

1.0, 







AIRFSPCW 

= 

-12 . , 

-11., 

-10., -9., 2*-8 . , -9., 





-9.5, 

3*-9 . 

, 3*-10 . , 9 *-9 . 0 , -8.0, 



AIRFDILE 

= 

8*1.4 

, 1.0, 

0.7, 0.5, 0.3, 0.0, -0.2, 

-0.5, 

-0.8, 

1.0, 







AIRFSPLE 

= 

-12 . , 

-11., 

-10., -9., 2*-8 . , -9., 





-9.5, 

3 *-9 . 

, 3*-10 . , 9*-9 . 0 , -8.0, 



AIRFDITEI 

= 

8*1.4 

, 1.0, 

0.7, 0.5, 0.3, 0.0, -0.2, 

-0.5, 

-0.8, 

1.0, 







AIRFSPTEI 

= 

-12 . , 

-11., 

-10., -9., 2 *-8 . , -9., 





-9.5, 

3 *-9 . 

, 3*-10 . , 9 *-9 . 0 , -8.0, 



AIRFDITEO 

= 

8*1.4 

, 1.0, 

0.7, 0.5, 0.3, 0.0, -0.2, 

-0.5, 

-0.8, 

1.0, 







AIRFSPTEO 

= 

-12 . , 

- 11 ., 

-10., -9., 2 *-8 . , -9., 





-9.5, 

3*-9 . 

, 3*-10 . , 9 *-9 . 0 , -8.0, 



AIRFDIHT 

= 

8*1.4 

, 1.0, 

0.7, 0.5, 0.3, 0.0, -0.2, 

-0.5, 

-0.8, 

1.0, 







AIRFSPHT 

= 

-12., 

- 11 ., 

-10., -9., 2 *-8 . , -9., 





-9.5, 

3*-9 . 

, 3*-10 . , 9 *-9 . 0 , -8.0, 



AIRFDIVT 

= 

8*1.4 

, 1.0, 

0.7, 0.5, 0.3, 0.0, -0.2, 

-0.5, 

-0.8, 

1.0, 







AIRFSPVT 

= 

-12 . , 

- 11 ., 

-10., -9., 2*-8 . , -9., 





-9.5, 

3*-9 . 

, 3*-10 . , 9 *-9 . 0 , -8.0, 



AIRFDIMG 

= 

17*0. 

0, 




AIRFSPMG 

= 

24*0. 

0, 




AIRFDING 

= 

17*0. 

0, 




AIRFSPNG 

= 

24*0. 

0 





$END 


$FANIN 


AFAN 

= 

16.98436 

r 

FANRPM 

= 

3609.0 

r 

NB1 

= 

27 

r 

UT 

= 

923 . 13 

r 

I SHAPE 

= 

0 

r 

IDOOR 

$END 

$COREIN 


0 


AB 

= 

1111.27 

/ 

BL 

= 

0.75 

r 

ICR 

= 

2 

/ 

NFN 

= 

90 

r 

PLD 

= 

0.15 

r 

P3 

= 

92.0 

r 

P4 

= 

85.2 

r 

P7 

= 

20.28 

r 

T3 

= 

940.7 

r 

T4 

= 

1880.5 

r 

T7 

= 

1201.5 

/ 

W3 

= 

141.2 

r 

WFP 4 
$END 

= 

3.9999999E 

-02 
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$TURIN 


CSS 

= 

151.42 

MDOT 

= 

200.20 

NBLADE 

= 

88 

TS 

= 

1172.34 

TURCAL 

= 

o 

o 

TURRPM 

= 

3609.0 

VTR 

= 

937.8076 

$END 



$ JN8C4 IN 



APRI 

= 

8.476 

TPRI 

= 

989.84 

VPRI 

= 

1048.88 

WPRI 

= 

327.87 

ATHP 

= 

7.474 

PLUGD 

= 

o 

o 

VSEC 

= 

369.42 

WSEC 

= 

316.23 

ASEC 

= 

12.352 

TSEC 

= 

514.51 

SPOKE 

= 

20. , 

PEN 

= 

0.925 , 

D I VANG 

= 

o 

o 

EJL 

= 

11.0 

EJD 

= 

5.15 

XMAR 

= 

0.98 

EJLIN 

= 

7.09 , 

FLIN 

= 

2500.0 , 

EJASP 

= 

1.42 

CFG 

= 

0.95 

PSI0 

= 

o 

o 

DELPE 

= 

o 

o 

DELPI 

= 

o 

o 

DELMIX 

= 

o 

o 

$END 



$AIRFIN 




ICWING = 2, 
IVTAIL = 1, 
IHTAIL = 1, 
ILEDGE = 1, 
ITEIN = 1, 
ITEOUT = 1, 
INGEAR = 1, 
IMGEAR = 1, 

AW = 8878.0, 

BW = 154.3, 

AV = 435.0, 

BV = 16.8, 

AH = 726.0, 

BH = 71.0, 

AFIN = 100.0, 
BFIN = 10.0, 
FLAIN = 10.0, 
NFSIN = 1, 
AFOUT = 200.0, 
BFOUT = 20.0, 
FLAOUT = 10.0, 



NFSOUT 
NNG = 
NWNG = 
DNG = 
RLNG = 
NMG = 
NWMG = 
DMG = 
RLMG = 
$END 


= 1 , 
r 

2 , 

• 3 , 
13.7 

I 

8 , 

• 7 , 

16.0 



=tt= =ft= =tt= =H= % in 


1.2 Sideline 

Sideline test case for the HSRNOISE code 
Free-field predictions without Deltas 
$GEOM 

HMIC = 0.0, 

SL = 1476.0 

$END 

$WEATHER 

PSIAMB = 14.7 , 

RH = 70.0 , 

TFAMB = 77.0 

$END 

$FLIPATH 

ALTEVO = 700.0000 , 

FPA = 8.100000 , 

PAE = 17.00000 , 

VAIR = 374.5672 

$END 
$ SOURCE 

IAFRAM = F , 

IATM = 1 , 

IEGA = F, 

IGRND = F, 

I FAN = T, 

ICORE = T, 

I JET = T, 

ITUR = T, 

NENG = 1 

$END 
$DELTAS 

FANDIBB = 17*0.0 , 

FANSPBB = 24*0.0 , 

FANBPF = 0.0 , 

FAN2BPF = 0.0 , 

FAN3BPF = 0.0 , 

FAN4BPF = 0.0 , 

CORDIBB = 17*0.0, 

CORSPBB = 24*0.0 , 

TURDIBB = 17*0.0, 

TURSPBB = 24*0.0 , 

TURDITN = 17*0.0 , 

TURSPTN = 24*0.0 
$END 
$FANIN 

AFAN = 16.98436 , 

FANRPM = 5618.0 , 

NB1 = 27 , 

UT = 1431.22 

I SHAPE = 0 , 

IDOOR = 0 

$END 
$COREIN 

AB = 1111.27 , 
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BL 

ICR 

NFN 

PLD 

P3 

P4 

P7 

T3 

T4 

T7 

W3 

WFP 4 
$END 
$TURIN 
CSS 
MDOT 

NBLADE = 
TS 

TURCAL = 
TURRPM = 
VTR 
$END 

$ JN8C4 IN 
APRI 
TPRI 
VPRI 
WPRI 
ATHP 
PLUGD 
VSEC 
WSEC 
ASEC 
TSEC 
SPOKE 
PEN 

DIVANG = 
EJL 
EJD 
XMAR 
EJLIN 
FLIN 
EJASP 
CFG 
PSIO 
DELPE 
DELPI 
DELMIX = 
$END 


0.75 

2 

330 , 

0.15 

294.5 

274.6 
50.15 

1362.1 , 

3239.4 

1992.8 

340.8 

3 . 9999999E-02 


151.42 

487.50 

88 

1919.03 

0.0 

5618.0 

1508.941 


8.48 
1668.84 
2467 . 14 

698.95 
7.47 
0.0 

564.42 
440.94 

11.95 
546.03 
20 . 

0.925 , 

0.0 
11.0 
5.15 

0.919 , 

7.09 

2500.0 , 

1.42 
0.95 
0.0 
0.0 , 

0.0 , 

0.0 
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5.1.3 Free-field Centerline 


Centerline test case for the HSRNOISE code 
Free-field predictions without Deltas 
$GEOM 

HMIC =0.0 , 

SL =0.0 

$END 

$WEATHER 

PSIAMB = 14.7 , 

RH = 70.0 , 

TFAMB = 77.0 

$END 

$FLIPATH 

ALTEVO = 1273.0 , 

FPA = 2.30 , 

PAE = 11.7 , 

VAIR = 377.7744 

$END 
$ SOURCE 
IAFRAM = F , 

IATM = 1 , 

IEGA = F, 

IGRND = F, 

I FAN = T, 

ICORE = T, 

I JET = T, 

ITUR = T, 

NENG = 1 

$END 
$ DELTAS 

FANDIBB = 17*0.0 , 

FANSPBB = 24*0.0 , 

FANBPF = 0.0 , 

FAN2BPF = 0.0 , 

FAN3BPF = 0.0 , 

FAN4BPF = 0.0 

CORDIBB = 17*0.0 , 

CORSPBB = 24*0.0 , 

TURDIBB = 17*0.0 , 

TURSPBB = 24*0.0 , 

TURDITN = 17*0.0 , 

TURSPTN = 24*0.0 
$END 
$FANIN 

AFAN = 16.98436 , 

FANRPM = 4898.000 , 

NB1 = 27, 

UT = 1249.580, 

I SHAPE = 0 , 

IDOOR = 0 

$END 
$COREIN 
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AB 

BL 

ICR 

NFN 

PLD 

P3 

P4 

P7 

T3 

T4 

T7 

W3 

WFP 4 
$END 
$TURIN 
CSS 
MDOT 
NBLADE 
TS 

TURCAL 

TURRPM 

VTR 


1111.27 , 

0.75 

2 

210 

0.15 

193.9 

180.4 

34.41 

1189.8 , 

2672.5 , 

1638.1 , 

249.1 

3 . 9999999E-02 


151.42 
354.8 , 

88 , 

1580.61 , 

0.0 

4898.0 

1311.559 


$END 

$ JN8C4 IN 
APRI 
TPRI 
VPRI 
WPRI 
ATHP 
PLUGD 
VSEC 
WSEC 
ASEC 
TSEC 
SPOKE 
PEN 

D I VANG = 

EJL 

EJD 

XMAR 

EJLIN 

FLIN 

EJASP 

CFG 

PSIO 

DELPE 

DELPI 

DELMIX = 

$END 


8.476 , 

1350.24 

1921.09 
541.20 
7.47 , 
0.0 

643.98 
475.90 
12.084 
544.19 
20 . , 

0.925 , 

0.0 , 
11.0 
5.15 

0.979 , 

7.09 

2500.0 , 

1.42 

0.95 

0.0 

0.0 , 

0.0 , 

0.0 
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5.2 Full System Input Files 


The full system input files include ground reflections, EGA and jet shielding. Spectral 
corrections for all noise sources are included. 


5.2.1 Full System Approach 


Approach test case for the HSRNOISE code 
Full system predictions with Deltas 
$GEOM 

HMIC =4.0 , 

SL =0.0 

$END 

$WEATHER 

PSIAMB = 14.7 , 

RH = 70.0 , 

TFAMB = 77.0 , 

SIGMA = 1.1E6 

$END 

$FLIPATH 

ALTEVO = 394.0 , 

FPA = -3.0 , 

PAE = 7.1 , 

VAIR = 267.7168 

$END 
$ SOURCE 

I AFRAM = T , 

IATM = 1 , 

IEGA = F, 

IGRND = T, 

I FAN = T, 

ICORE = T, 

I JET = T, 

ITUR = T, 

NENG = 4 

$END 
$DELTAS 

FANDIBB = 17*-8 . 5 , 

FANSPBB = 24*1.0 , 

FANBPF = -8.5 , 

FAN2BPF = -8.5 , 

FAN3BPF = -8.5 , 

FAN4BPF = -8.5 , 

CORDIBB = 17*0.0 , 

CORSPBB = 24*0.0 , 

TURDIBB = 17*0.0 , 

TURSPBB = 24*0.0 , 

TURDITN = 17*0.0 , 

TURSPTN = 24*0.0 , 

AIRFDICW = 8*1.4, 1.0, 0.7, 0.5, 0.3, 0.0, -0.2, -0.5, -0.8, - 

1.0, 
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AIRFSPCW 


= -12., - 11 ., -10., -9., 2 *-8 . , -9., 

-9.5, 3*-9 . , 3*-10 . , 9*-9 . 0 , -8.0, 

AIRFDILE = 8*1.4, 1.0, 0.7, 0.5, 0.3, 0.0, -0.2, -0.5, -0.8, - 

1.0, 

AIRFSPLE = -12., -11., -10., -9., 2*-8., -9., 

-9.5, 3*-9 . , 3*-10 . , 9 *-9 . 0 , -8.0, 

AIRFDITEI = 8*1.4, 1.0, 0.7, 0.5, 0.3, 0.0, -0.2, -0.5, -0.8, - 

1.0, 

AIRFSPTEI = -12., -11., -10., -9., 2*-8., -9., 

-9.5, 3*-9 . , 3*-10 . , 9 *-9 . 0 , -8.0, 

AIRFDITEO = 8*1.4, 1.0, 0.7, 0.5, 0.3, 0.0, -0.2, -0.5, -0.8, - 

1.0, 

AIRFSPTEO = -12., -11., -10., -9., 2*-8., -9., 

-9.5, 3*-9 . , 3*-10 . , 9 *-9 . 0 , -8.0, 

AIRFDIHT = 8*1.4, 1.0, 0.7, 0.5, 0.3, 0.0, -0.2, -0.5, -0.8, - 

1.0, 

AIRFSPHT = -12., -11., -10., -9., 2*-8., -9., 

-9.5, 3 *-9 . , 3*-10 . , 9 *-9 . 0 , -8.0, 

AIRFDIVT = 8*1.4, 1.0, 0.7, 0.5, 0.3, 0.0, -0.2, -0.5, -0.8, - 

1.0, 

AIRFSPVT = -12., -11., -10., -9., 2*-8., -9., 

-9.5, 3 *-9 . , 3*-10 . , 9*-9 . 0 , -8.0, 

AIRFDIMG = 17*0.0, 

AIRFSPMG = 24*0.0, 

AIRFDING = 17*0.0, 

AIRFSPNG = 24*0.0 

$END 
$FANIN 


AFAN 

= 

16.98436 

/ 

FANRPM 

= 

3609.0 

/ 

NB1 

= 

27 

/ 

UT 

= 

923 . 13 

/ 

I SHAPE 

= 

1 

/ 

IDOOR 

= 

0 


$END 




$COREIN 




AB 

= 

1111.27 

r 

BL 

= 

0.75 

r 

ICR 

= 

2 

r 

NFN 

= 

90 

r 

PLD 

= 

0.15 

r 

P3 

= 

92.0 

r 

P4 

= 

85.2 

r 

P7 

= 

20.28 

/ 

T3 

= 

940.7 

r 

T4 

= 

1880.5 

r 

T7 

= 

1201.5 

r 

W3 

= 

141.2 

r 

WFP 4 

= 

3.9999999E 

-02 

$END 




$TURIN 




CSS 

= 

151.42 

/ 

MDOT 

= 

200.20 

/ 

NBLADE 

= 

88 

f 

TS 

= 

1172.34 

r 

TURCAL 

= 

10.0 

r 

TURRPM 

= 

3609.0 

/ 
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VTR 

= 

937.8076 

$END 

$ JN8C4 IN 

APRI 

= 

8.476 

TPRI 

= 

989.84 

VPRI 

= 

1048.88 

WPRI 

= 

327.87 

ATHP 

= 

7.474 

PLUGD 

= 

0.0 

VSEC 

= 

369.42 

WSEC 

= 

316.23 

ASEC 

= 

12.352 

TSEC 

= 

514.51 

SPOKE 

= 

20. , 

PEN 

= 

0.925 , 

D I VANG 

= 

0.0 

EJL 

= 

11.0 

EJD 

= 

5.15 

XMAR 

= 

0.98 

EJLIN 

= 

7.09 

FLIN 

= 

2500.0 , 

EJASP 

= 

1.42 

CFG 

= 

0.95 

PSI0 

= 

0.0 

DELPE 

= 

0.2 , 

DELPI 

= 

-1.0 , 

DELMIX 

= 

0.1 

$END 

$AIRFIN 

ICWING 

= 2, 



IVTAIL = 1, 
IHTAIL = 1, 
ILEDGE = 1, 
ITEIN = 1, 
ITEOUT = 1, 
INGEAR = 1, 
IMGEAR = 1, 

AW = 8878.0, 

BW = 154.3, 

AV = 435.0, 

BV = 16.8, 

AH = 726.0, 

BH = 71.0, 

AFIN = 100.0, 
BFIN = 10.0, 
FLAIN = 10.0, 
NFSIN = 1, 
AFOUT = 200.0, 
BFOUT = 20.0, 
FLAOUT = 10.0, 
NFSOUT = 1, 

NNG = 1, 

NWNG = 2, 

DNG = 3.3, 

RLNG = 13.7, 
NMG = 2, 

NWMG = 8, 



DMG = 3.7, 
RLMG = 16.0, 
$END 


195 





5.2.2 Full System Sideline 


Sideline test case for the HSRNOISE code 
Full system predictions with Deltas 
$GEOM 

HMIC = 4.0, 

SL = 1476.0 

$END 

$WEATHER 

PSIAMB = 14.7 , 

RH = 70.0 , 

TFAMB = 77.0 , 

SIGMA = 1.1E6 
$END 

$FLIPATH 

ALTEVO = 700.0000 , 

FPA = 8.100000 , 

PAE = 17.00000 , 

VAIR = 374.5672 

$END 
$ SOURCE 

IAFRAM = F , 

IATM = 1 , 

IEGA = T, 

IGRND = T, 

I FAN = T, 

ICORE = T, 

I JET = T, 

I JETSHLD = T, 

ITUR = T, 

NENG = 4 

$END 
$ DELTAS 

FANDIBB = 17*-4 . 5 , 

FANSPBB = 24*1.1 , 

FANBPF = -4.5, 

FAN2BPF = -4.5, 

FAN3BPF = -4.5, 

FAN4BPF = -4.5, 

CORDIBB = 4*-0.1,-0.4,-0.7,-1.0,4*-l.l,-1.0, 
-0.9, -0.8, -0.6, -0.5, -0.3, 

CORSPBB = 24*0.0 , 

TURDIBB = 4*-0.1,-0.4,-0.7,-1.0,4*-l.l,-1.0, 
-0.9, -0.8, -0.6, -0.5, -0.3, 

TURSPBB = 24*0.0 , 

TURDITN = 4*-0.1,-0.4,-0.7,-1.0,4*-l.l,-1.0, 
-0.9, -0.8, -0.6, -0.5, -0.3, 

TURSPTN = 24*0.0 
$END 
$FANIN 

AFAN = 16.98436 , 

FANRPM = 5618.0 , 

NB1 = 27 , 
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UT 

= 

1431.22 

r 

I SHAPE 

= 

1 

r 

IDOOR 

SND 

:orein 


1 


AB 

= 

1111.27 

r 

BL 

= 

0.75 

r 

ICR 

= 

2 

r 

NFN 

= 

330 

r 

PLD 

= 

0.15 

r 

P3 

= 

294.5 

r 

P4 

= 

274.6 

r 

P7 

= 

50.15 

r 

T3 

= 

1362 . 1 

/ 

T4 

= 

3239.4 

r 

T7 

= 

1992.8 

r 

W3 

= 

340.8 

r 

WFP 4 

= 

3.9999999E 

-02 


$END 

$TURIN 

CSS 

MDOT 

NBLADE 

TS 

TURCAL 

TURRPM 

VTR 

$END 

$ JN8C4IN 
APRI 
TPRI 
VPRI 
WPRI 
ATHP 
PLUGD 
VSEC 
WSEC 
ASEC 
TSEC 
SPOKE 
PEN 

DIVANG 

EJL 

EJD 

XMAR 

EJLIN 

FLIN 

EJASP 

CFG 

PSIO 

DELPE 

DELPI 

DELMIX 

$END 


151.42 

487.50 

88 

1919.03 

10.0 

5618.0 

1508.941 


8.48 
1668.84 
2467 . 14 

698.95 
7.47 
0.0 

564.42 
440.94 

11.95 
546.03 
20 . 

0.925 , 

0.0 
11.0 
5.15 

0.919 , 

7.09 

2500.0 , 

1.42 
0.95 
0.0 

- 1.8 , 
-1.3 , 
- 0.1 
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5.2.3 Full System Centerline 


Centerline test case for the HSRNOISE code 


Full system predictions with Deltas 
$GEOM 


HMIC 

= 

4.0 

SL 

= 

0.0 

$END 



$WEATHER 



PSIAMB 

= 

14.7 , 

RH 

= 

70.0 

TFAMB 

= 

77.0 

SIGMA 

= 

1 . 1E6 

$END 



$FLIPATH 



ALTEVO 

= 

1273.0 , 

FPA 

= 

2.30 

PAE 

= 

11.7 

VAIR 

= 

377.7744 

$END 



$ SOURCE 



IAFRAM 

= 

F, 

IATM 

= 

1, 

IEGA 

= 

F, 

IGRND 

= 

T, 

I FAN 

= 

T, 

ICORE 

= 

T, 

I JET 

= 

T, 

ITUR 

= 

T, 

NENG 

= 

4 

$END 



$ DELTAS 



F AND IBB 

= 

17*-6 . 7 , 

FANSPBB 

= 

24*1.1 , 

FANBPF 

= 

-6.7, 

FAN2BPF 

= 

-6.7, 

FAN3BPF 

= 

-6.7, 

FAN4BPF 

= 

-6.7, 

CORDIBB 

= 

4 *-0.1, -0.4, -0.7, -1.0, 4 *-1.1, -1.0, 



-0.9, -0.8, -0.6, -0.5, -0.3, 

CORSPBB 

= 

24*0.0 , 

TURDIBB 

= 

4 *-0.1, -0.4, -0.7, -1.0, 4 *-1.1, -1.0, 



-0.9, -0.8, -0.6, -0.5, -0.3, 

TURSPBB 

= 

24*0.0 , 

TURDITN 

= 

4 *-0.1, -0.4, -0.7, -1.0, 4 *-1.1, -1.0, 



-0.9, -0.8, -0.6, -0.5, -0.3, 

TURSPTN 

= 

24*0.0 

$END 



$FANIN 



AFAN 

= 

16.98436 

FANRPM 

= 

4898.000 

NB1 

= 

27, 

UT 

= 

1249.580, 
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= 1 , 
= 0 


I SHAPE 
IDOOR 
$END 
$COREIN 
AB 
BL 
ICR 
NFN 
PLD 
P3 
P4 
P7 
T3 
T4 
T7 
W3 

WFP 4 
$END 
$TURIN 
CSS 
MDOT 
NBLADE = 
TS 

TURCAL = 
TURRPM = 
VTR 
$END 

$ JN8C4 IN 
APRI 
TPRI 
VPRI 
WPRI 
ATHP 
PLUGD 
VSEC 
WSEC 
ASEC 
TSEC 
SPOKE 
PEN 

DIVANG = 
EJL 
EJD 
XMAR 
EJLIN 
FLIN 
EJASP 
CFG 
PSIO 
DELPE 
DELPI 
DELMIX = 
$END 


1111.27 , 

0.75 

2 

210 

0.15 

193.9 

180.4 

34.41 

1189.8 , 

2672.5 , 

1638.1 , 

249.1 

3 . 9999999E-02 


151.42 

354.8 

88 , 

1580.61 , 

10.0 
4898.0 
1311.559 


8.476 

1350.24 

1921.09 
541.20 
7.47 
0.0 

643.98 
475.90 
12.084 
544.19 
20 . , 
0.925 , 

0.0 
11.0 
5.15 

0.979 , 

7.09 

2500.0 , 
1.42 

0.95 , 

0.0 

- 1.8 , 
- 1.1 , 

0.1 
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